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Abstract 



This thesis examines the topics of disorder and electron-electron interactions 
in three distinct quantum systems. Firstly, the Anderson transition is studied for 
the body centred cubic and face centred cubic lattices. We obtain high precision 
results for the critical disorder at the band centre and the critical exponent using 
the transfer-matrix method and finite size scaling. Comparing the critical disorder 
between the simple cubic, body centred cubic and face centred cubic lattices, an 
increase in the critical disorder is observed as a function of the coordination number 
of the lattice. The critical exponent is found to be ~ 1.5 in agreement with the 
value for the simple cubic lattice. Energy-disorder phase diagrams are plotted for 
both lattice types. 

Next, we consider the Aharonov-Bohm effect for an exciton in a ID ring 
geometry. The aim is to determine how the addition of a constant electric field 
in the plane of the ring effects the Aharonov-Bohm oscillations, which occur as 
a function of the magnetic flux threading the ring. We develop a self consistent 
equation for the ground state energy, which is then solved numerically. Oscillations 
in the ground state energy have an increasing amplitude as a function of electric 
field strength until a critical electric field value. At this point, oscillations in the 
oscillator strength become inverted, with the oscillation minimum reaching zero at 
half a magnetic flux quantum. This suggests a possible process for controlling the 
formation and recombination of excitons through tuning the applied fields. 

The final and largest section of the thesis is concerned with collective excita- 
tions of graphene in a strong perpendicular magnetic field. The excitations, which 
are most strongly mixed arc identified and used as a basis to diagonalise the Hamil- 
tonian, which includes the Coulomb interaction between electrons and holes. In this 
way the oscillator strengths and energies of collective excitations are obtained. The 
good quantum numbers for collective excitations are identified. In particular, we 
study those arising from the SU(4) symmetry, which is due to two spin and two 
valley pseudospin projections. This enables us to determine the multiplet structure 
of the states. In addition to neutral collective excitations or excitons, we investigate 
the possible formation of charged collective excitations or trions from nearly full or 
nearly empty Landau levels. The localisation of neutral collective excitations upon 
a single Coulomb or (^-function impurity is also examined. 
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Chapter 1 

Introduction 



The behaviour of electrons in quantum systems depends on the way they interact 
with their surroundings, their interactions with other electrons and the presence of 
external fields. In this work, we shall explore these themes in three different systems. 
Firstly we examine Anderson localisation [I] in three-dimensional crystalline solids. 
This neglects electron-electron (e-e) interactions, but takes into account the effect 
of lattice imperfections, commonly referred to as disorder, which occur in any real 
crystalline material. There exists a critical disorder above which the electronic wave 
function becomes localised, due to quantum interference and the material displays 
insulating behaviour. Secondly, we consider the Aharonov-Bohm (AB) effect [2], 
which is another quantum interference phenomenon, exhibited by systems including 
an external magnetic field in a certain geometry. It provides a pleasing illustration 
of how the magnetic vector potential is not merely a mathematical convenience and 
can actually influence a system's physical properties. The excitonic AB effect is of 
particular interest, due to competing effects of the magnetic field and electron-hole 
(e-h) interactions. Thirdly we study graphene [3], a two-dimensional carbon system, 
which has been the subject of an immense research effort in the past few years. The 
work presented here, which comprises the bulk of the thesis, focuses on the collective 
behaviour of its electrons in a strong magnetic field in both the presence and absence 
of impurities. Much of the theory developed for the two-dimensional electron gas 
can be used to describe graphene, although there are some important differences 
between the systems. Each of these fields is introduced more fully in the following. 
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(a) (b) (c) 

Figure 1.1: (Images courtesy of A. Rodriguez, R.A. Romer, K. Slevin and L.J. A. 
Vasquez). Three-dimensional electronic wave function plots for a simple cubic lattice 
at the band centre for (a) an extended state (W = 1), (b) a state at the critical 
disorder (W = Wc = 16.5) and (c) a localised state {W = 20). In (a) and (c), the 
system has 100^ lattice sites and in (b) 240'^ lattice sites. At a lattice site with index 
j, the wave function amplitude is denoted by ipj and the size of the box is directly 
proportional to The colour of boxes indicates which slice of the system the 

lattice site is in i.e. it clarifies the third dimension. 

1.1 Anderson localisation 
1.1.1 Background and basic ideas 

Understanding electronic transport in crystalline solids is one of the main goals of 
condensed matter theory. Paul Drude was first to seriously tackle the problem in the 
early 1900s by using classical kinetic theory. Following this came various quantum 
mechanical models such as the free electron model, Bloch's theorem for electrons in 
a periodic potential and the tight binding method. One very fundamental question 
to ask is: what makes some materials metallic whilst others are electrically insu- 
lating? Electronic band structure calculations go some way towards answering this 
question. For example, a material with a completely filled valence band separated 
from the conduction band by an energy gap large enough to make excitation be- 
tween the two bands improbable, would be expected to exhibit insulating behaviour. 
However, such band theories do not include e-e interactions and consequently pre- 
dict some materials, known as Mott insulators, to be conductors when they are in 
fact insulators [4||6]. Another important aspect affecting the itinerancy of electrons 
is disorder. Here I refer not to the structural disorder of amorphous lattices, rather 
to the compositional disorder present in crystalline solids. This may be caused by 
dislocated or missing atoms (vacancies) and impurity atoms. 
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Density of states 




Localised 



Figure 1.2: Rough sketch of the density of states in the presence of disorder. The 
shaded region represents localised states and the unshaded region extended states; 
these are separated by the mobility edge. 



The framework for such a theory was first laid down by P.W. Anderson 
in 1958 |ljj. In this seminal paper he models electronic transport in a disordered 
lattice via quantum jumps between localised states. This model is only reasonable 
at temperatures below about lOK, where all scattering is elastic. Both electron- 
phonon and e-e interactions are ignored. Electrons are treated as spinless particles. 
The disorder is introduced via a random potential at the lattice sites with the 
disorder parameter, W, characterising the width of the probability distribution of 
onsite potential energies. In 3D there is a critical disorder at which there is a phase 
transition between extended and exponentially localised states of the electron - 
the disorder-induced metal-insulator transition (MIT) or Anderson transition (AT) 
[7|[8]. Figure [lT] shows some examples of squared wave function amplitude plots 
for extended, critical and localised wave functions at the band centre for different 
disorder values. I use Wc to denote the critical disorder at the band centre. 

Let us discuss this AT in more detail. In the absence of disorder, in a 
perfectly periodic potential due to the ionic lattice, all electronic states are Bloch 
states and extend throughout the system. Hence the system has a finite conductivity 
as the absolute zero of temperature is approached and by this definition behaves 
as a metal. Upon the introduction of disorder, localised states form at the band 
edges. A schematic diagram of the density of states for a fixed disorder value below 



Wc is shown in Fig. 1.2 from which it can be seen that localised and extended 
states exist in different regions of the band. The lines of constant energy separating 
them are called mobility edges Thus we see that for a fixed disorder below 
Wc, there is a localisation-delocalisation transition as a function of energy. As the 
degree of disorder is increased, more states become localised and the mobility edges 
move towards the band centre. They meet at the band centre for critical disorder. 
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Wc, at which point the whole band consists of locaHsed states. The transport 
properties of the system depend on those states close to the Fermi energy. If the 
Fermi energy is well within the extended region, the system will behave metallically. 
Conversely if it is within the localised region, the system acts as an insulator and 
the conductivity will vanish in the limit of absolute zero temperature. Keeping 
the Fermi energy constant and increasing the disorder from zero will result in a 
transition from metallic to insulating behaviour at the critical disorder. 

The mechanism for the localisation of an electron in the presence of disorder 
is quantum interference of the electronic wave function. Classically one would expect 
an electron to be itinerant if its energy exceeds the average lattice site potential and 
localised if it possesses a lower energy. In quantum mechanics the situation is much 
more complicated. The electron acts as a wave packet, so that upon encountering 
a potential barrier, it splits into a transmitted and reflected wave. In this way the 
electronic wave function may interfere with itself and this results in an increased 
probability of backscattering. This has been studied particularly in the limit of weak 
disorder [10] , where constructive interference between clockwise and anticlockwise 
paths results in a positive correction to the resistivity: the weak localisation correc- 
tion. Such behaviour is the precusor to the strong Anderson localisation discussed 
here, when the interference is so strong that the electronic propagation is prohibited 
altogether. 



1.1.2 Scaling theory 

After Anderson's original paper, the problem was re-expressed in the language of 
renormalisation group theory and it was shown that the AT can be described as a 
second order phase transition |11[|12| . In 1979 Abrahams et al. developed the one- 
parameter scaling theory (OPST) of localisation [13| , which examines the scaling 
of the conductance 5 of a hypercube of volume L*^ with the system size L . The 
key assumption of the theory is that the quantity /3 = ^f^||y depends solely on 
the conductance g (L) and not separately on L. The qualitative dependence of j3 



upon g is shown in Fig. 1.3 and estimated by interpolating between the known 



forms in the limits of large and small g 14 . In the limit of small g, the electronic 



wave function is localised, so g = goe . Here A is the localisation length, which 
is the exponential decay length of the envelope of the wave function. Hence for 



small 5, /3 (5) = ln(g(/(7o) (see Fig. 1.3). For large 5, there is metallic behaviour, 
so g = aL'^^'^, where a is the conductivity, yielding f3 (g) = d — 2. An MIT can 
take place at the critical disorder when /3 (g) = 0. Thus the main result of OPST 
is revealed: a transition is only possible in dimensions d > 2 and all states are 
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Figure 1.3: A qualitative plot of the scaling function f3 (g) against In (g) for different 
dimensions d. 



localised for any finite amount of disorder in lower dimensions. Despite this, many 
experiments reaching temperatures as low as lO^^K claim to have observed an MIT 
in effectively 2D systems such as silicon MOSFETs (metal oxide semiconductor field 
effect transisitors), where e-e interactions play a role |15| . The interaction strength 
may drive a transition and can be tuned by altering the electronic density. At higher 



densities the e-e interactions are weaker, which can result in metallic behaviour 16 



1.1.3 Characterising the Anderson transition 

Previously, the temperature dependence of the conductivity, a, was alluded to as 
a way of distinguishing between metallic and insulating behaviour. An alternative 
way, used most commonly in theoretical studies, is to examine the form of the wave 
function, in particular the localisation length, A. It tends to infinity in the thermo- 
dynamic limit as the critical disorder is approached from the insulating region. In 
practice, localisation lengths can only be calculated for finite system sizes; a mate- 
rial is metallic if its localisation length as a fraction of the system size increases with 
increasing system size and insulating if it decreases with system size. Fortunately 
both the conductivity [17] and the localisation length [T] are self-averaging quan- 
tities; the values obtained by averaging over different impurity configurations at a 
fixed disorder, correspond to the value that will be observed with highest probability 
in the thermodynamic limit. 

Studies of the AT seek to determine the critical disorder Wc and the critical 
exponents for the transition, s and ly, such that close to the transition (toc|VF — Wc\^ 
and ^oc|iy — Wc| where ^ is the correlation length for an infinite system. Values 
of these critical exponents determine the universality class of the system. One may 
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also consider a transition as a function of energy, E, at constant disorder, as should 
be clear from the explanation of mobility edges. In this case, the above relations 
describing the behaviour of a and ^ still hold (with the same critical exponents) but 
with W and Wc replaced by E and Ec respectively. It should be noted that the 
quantities Ec and Wc are not universal, depending, for example, on the distribution 
of onsite potential energies |18| and, as shown in the work presented in Chapter 2, 
on the lattice type. 

Sadly we are far from achieving a generally accepted analytical theory for 
Anderson localisation in d dimensions. However, for d = 2 + e, where e « 1, 
it has been shown using perturbation theory and field theoretic techniques that 



s = {d — 2) u 11 , so that s = in 3D. Early attempts to extend this approach to 
d dimensions solved a self consistent equation for the diffusion coefficient, resulting 
in u = 1/ {d — 2) |19| , which is equal to 1 in 3D. This is in agreement with many 
experimental studies, which reported u k l.O (20] - |24| . The most recent analytical 
study ^sj finds u = ^ + ^j^, agreeing more with numerical studies, which find 
K 1.5 for d = 3. 

A common numerical approach is to obtain localisation lengths via the transfer- 
matrix method (TMM) |7|; this is explained in Chapter 2. Another technique is to 



diagonalise the Hamiltonian directly and examine energy level statistics 26-29 
The distribution of nearest neighbour level spacings is exponential for localised 
states, whereas metallic states obey GOE (Gaussian orthogonal ensemble) statis- 
tics. The 3D eigenstate at criticality is multifractal, so that multifractal analysis 



techniques are also commonly employed 30 . Finite size scaling (FSS) is used to 



correct for the use of finite system sizes in an attempt to model physics in the ther- 



modynamic limit 18 . The exact values of critical parameters are still not agreed 



upon. Recent TMM studies have found z/ = 1.57 ± 0.02, Wc = 16.54 + 0.02 M 



and u = 1.515 + 0.033, Wc = 16.50 + 0.05 31 for the case of a box distribution 



of onsite potentials. See Ref. 32 for a review. Energy level statistics 33 yield 
= 1.45 + 0.08, Wc = 16.35. Multifractal analysis |34] finds i' = 1.58 + 0.03, 
Wc = 16.57 + 0.03. Clearly an accurate determination of the critical parameters at 
the transition is still an active area of research. 

1.1.4 Experiment 

There has been a lot of experimental work in which an MIT has been observed 



20 , 23 . However, a true AT may not have been experimentally observed in crys- 
talline solids. This is because e-e interactions cannot be sufficiently reduced and 
phonons are always present, so that a transition between extended and localised 



6 



o 

CO 



-•««•« I I I 



A/J-^.8 



c 
o 

o 

(£ 
O 



^ m ^ ^ m m m m 



c^j .::,3; ^ ^ 40- 



10 100 180 340 440 540 640 
Time (ms) 



A/J = 4.2 



Figure 1.4: (Taken from Ref. 35 ). Absorption images of a BEC in a quasiperiodic 
optical lattice at different times and disorder values. The disorder, A, is measured 
in units of the hopping parameter, J/h = 153Hz. 



behaviour cannot unambiguously be attributed to disorder alone. However, it has 
been seen in a broad range of other physical systems. Some time ago, light local- 
isation was observed due to scattering from suspended polystyrene balls [36] and 
more recently due to scattering from semiconductor powders [37]. The localisa- 
tion of ultrasound waves was also seen and the measurement of localisation lengths 
possible [38| . These systems have the advantage that interparticle interactions are 
automatically eliminated. 

For a long time it has been the goal to directly see Anderson localisation of 
matter waves. Imaging the local density of states provides a direct way to probe 
the electronic probability density. This has been done using scanning tunnelling 



microscopy in thin semiconductor samples 39 and also by scanning tunnelling 



spectroscopy in Quantum Hall systems 40 . However, the e-e interactions have 



been seen to be strong in these experiments. The last ten years has seen incred- 



ible advances in the ability to control gases of ultracold atoms 41 ; such systems 
have potential for use as quantum simulators to experimentally mimic theoretical 
Hamiltonians. They are particularly advantageous, since the disorder is tunable 
and interactions between atoms can be controlled and made negligible. Using these 
techniques has allowed the very recent observation of exponential spatial localisation 



of matter waves in ID 35,42 . Quantum kicked rotors, in which ultracold atoms are 
exposed to a pulsed laser, have been seen to be analogous to spatially disordered 
systems, opening another way to experimentally probe the AT [43]|44| . In Ref. [35| 
a transition was seen for a Bose-Einstein condensate (BEC) in a quasi-periodic ID 
optical lattice in an experimental realisation of the Aubry- Andre model [45|. Fig. 



1.4 shows absorption images of the BEC as a function of time after the potential 
initially confining the atoms is removed. The disorder increases from top to bottom. 
At large times, one can see that the condensate has spread out for low disorders, 
but remains confined for higher disorders, clearly indicating localisation. 



1.2 The Aharonov-Bohm effect 

The AB effect |2j is a purely quantum mechanical phenomenon, which cannot be 
understood using classical physics. I shall begin by illustrating it for a particular 



example. Let us consider the situation (shown in Fig. 1.5 ) where a beam of electrons 
is split in two, such that some of the electrons pass above a region which contains a 
magnetic field, while some pass below it and they recombine on the other side. It is 
important that the magnetic field is confined to a region from which the electrons 
are excluded. This could be achieved by using a tightly wound current carrying 
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Figure 1.5: Schematic diagram of the experimental set up required to view the 
AB effect. The magnetic field is confined to a region from which the electrons are 
excluded. 



solenoid. Each electron is subject to zero Lorentz force from the magnetic field, so 
that classically one might expect it to behave in exactly the same way as it would if 
there were no solenoid. However, each electron experiences a finite vector potential, 
A(r), where r is the position vector, since for any path C enclosing the solenoid, 

S)dr.A{r) = f dS. (V x A) = f dS.B = (1.1) 

J Js Js 

c 

where the flux $ is finite. Let V'o(r) be the electron's wave function and E its energy 
when there is no solenoid. In the presence of a solenoid, the Hamiltonian is 

H=^ (1.2) 



where e is equal to negative the charge on an electron and nie is the electron's mass. 

eA 

he ' 



One can show that for ip{r) = e*"('"Vo(r)) with Va — 



Hi^ = e'"("-) (^V^^ V'o(r) = E^;. (1.3) 

Hence the addition of a solenoid alters the wave function by a phase a(r). Let 
be the wave function for a beam that goes on one side of the solenoid and -02 that 
for a beam on the other. When the two beams recombine, the intensity of the total 
beam depends on |^i + V'2p5 which in turn depends on the phase difference 

Aa = ^ (Cdr.A(r) = = 27r^. (1.4) 
hcj he $0 



Here <I>o = is the universal flux quantum. Thus any measurement of the in- 
terference pattern will depend upon the flux even though the electrons never 
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experience the magnetic field. 

This very striking effect was first predicted by W. Ehrenberg and R. E. Siday 
in 1949 46 . However, the idea remained largely unnoticed until it was rediscovered 
ten years later by Y. Aharonov and D. Bohm [2]. The latter paper also describes 
an electric AB effect, whereby two electron beams experience different time varying 
electrical potentials without being in an electric field and their interference pattern 
depends on the field. It was later found by M. V. Berry that the phase given by 



Eq. (1.4) is a particular example of a more general geometrical phase, which a par- 
ticle's wave function acquires, upon completion of an adiabatic loop in parameter 
space |47|. Indeed the inclusion of such phases is currently revolutionising semiclas- 
sical transport theory [48| . 

The AB electron interference pattern has been seen in multiple experiments 
[49| - [52] , which use a variety of techniques to confine the magnetic field including 
whiskers, solenoids, toroidal magnets and the Meissner effect. The AB effect has 
also been found to affect the transport properties in materials with non simply con- 
nected geometries. Oscillations in the magnetoresistance of metal rings with period 



$0 have been predicted and experimentally observed 53 . The rings must be small 
enough that the electron's phase is not randomised by inelastic scattering events 
as the electron moves about the arm of the ring. However, it was argued that the 
thermodynamic properties may still have $0 periodicity for the case when the elastic 
mean free path is small compared to the size of the ring |54|. It is interesting to 
consider the AB effect in superconducting quantum rings. It can be shown using 
Ginzburg-Landau theory that the free energy oscillates as a function of ^ with min- 
ima at integer multiples of <I>o/2 55 . One can attribute the extra factor of 2 to the 
formation of Cooper pairs of electrons. Hence the flux threading a superconducting 
quantum ring is quantised in units of ^o/2 [52,56 ; this is achieved by supercurrents 
which generate an additional magnetic field to that applied externally. Prom this 
it follows that all thermodynamic quantities have this flux periodicity 57 . For ex- 
ample, the critical temperature for a normal-superconducting phase transition has 



been seen to be #o/2 periodic 58 ,59 



1.3 Two-dimensional electron gas in a magnetic field 

In the two-dimensional electron gas (2DEG) model, electron motion is confined to 
a plane, which is embedded in three dimensions, meaning that any electromagnetic 
fields present are not restricted to the plane. During the latter half of the previous 
century it was a popular research topic, largely motivated by its application to 
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heterostructures, where the electrons are trapped in a layer less than Inm thick 60 . 
In particular, electrons in the inversion layer of silicon MOSFETs, one of the main 
electronic components used in computers, behave in accordance with the 2DEG 
model. In addition, effects of fundamental physical interest such as the integer and 
fractional quantum Hall effects (6l] and Wigner crystallisation [62] , have been found 
to exist in two-dimensional electron systems. The 2DEG model can also be used 
to describe electrons in thin films and to some extent in graphene. Hence it is now 
appropriate to examine the single particle problem for the 2DEG in a perpendicular 
magnetic field. Also, in Chapter |4j when the many body graphene problem in a 
magnetic field is addressed, the techniques used will be inspired by the interacting 
2DEG model. We consider the 2DEG in a portion of the x — y plane of area S, 
with a perpendicular magnetic field, B = Be^, where is a unit vector in the 
z direction. For most of the work presented here it will be convenient to use the 
symmetric gauge, Ag (r) = x r, but for some problems use of the Landau gauge, 
Al (r) = BxBy, will make for a simpler solution. We review below the single particle 



problem in both gauges 63 



1.3.1 Single particle problem in the symmetric gauge 

The symmetric gauge is the appropriate choice for systems with a disk geometry. 
The single particle Hamiltonian for an electron can be shown to be 

^|g = ^c''''(«^«+0, (1-5) 
where and a are given (up to a phase) by: 

Here z = x + iy, £b = \ r!^ the magnetic length and wJ^^ = is the cyclotron 
resonance frequency. The eigenkets of are 

|nm>=-i=(at)"(6t)'"|00>, (1.7) 
-s/nlm] 



where 



1/2 ^„ d \ .+ 1 f z* 



and n, m are non-negative integers. The state, |00), satisfies a|00) = 6|00) = 0, from 
which its coordinate form can be deduced. In the coordinate representation, the 
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eigenstates for general n, m ar^ 



<r|nm> = Vnm (r) 



' \ n—m / 2 



1 


/ ^2 


exp ( 


."44 











where is a generalised Laguerre polynomial. The eigenenergies corresponding to 
tpnm (r) are 

En = fuvf^ (n + (1.10) 

From this we see that electrons occupy discrete energy states equally sep- 
arated by the cyclotron resonance energy (CRE), huj^^ . These are Landau levels 
(LLs) and are labelled by the LL index n. The LLs are degenerate with respect 
to the m quantum number, which is called the oscillator quantum number. The 
filling factor v is defined by = where N is the total number of electrons in 
the system and Nq = = ^ is the LL degeneracy associated with m. From Eq. 



(1.7) it is clear that (a) raises (lowers) an electron by one LL and (b) raises 
(lowers) an electron's m quantum number within the same LL. The operators a, 
and b, b^ thus form two independent sets of ladder operators. Each set obeys 
bosonic commutation relations i.e. [a, a^] = [i, = 1. 

It is instructive to consider the physical meaning of the n and m quantum 



numbers. Fig. 1.6 indicates how the electron coordinate r may be decomposed 
according to r = R + X, where X is the guiding centre and R is the relative 
coordinate. One can show 

<R2> = (2n + l)4, <X2> = (2m + l)4, (r^) = (n + m + 1) 24. (1.11) 

Hence if we consider the electrons classically as particles undergoing circular motion 
in a magnetic field, the quantum number n determines the radius of the circular 
motion, whereas the quantum number m determines the distance of the guiding 
centre, about which the circular motion occurs, from the origin. The averages in 



Eq. (1.11) are taken with respect to the state given in Eq. (1.7). By setting n = in 
Eq. (1.11), we can understand the physical meaning of the magnetic length is] it is 
the root mean square radius of the cyclotron motion of an electron in the degenerate 
ground state. 



^This may be proved using Eq. ( 1.7 1 and induction. 
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Figure 1.6: Diagram to explain the decomposition of the electron coordinate into a 
guiding centre, X, and a relative coordinate, R. The shaded region indicates the 
area where an electron in a given LL with a given value of the oscillator quantum 
number, m, will exist. 



The orbital angular momentum projection is well defined for the states |nm). 
It is easily shown that the operator of angular momentum projection is 

4 = r X -ihV = h {a' a - b^) , (1.12) 



which clearly commutes with the Hamiltonian in Eq. (1.5) and yields an eigenvalue 
"iz = {n — m) h for the state |nm). 

In the equivalent problem for a hole, the single particle Hamiltonian is the 



same as that for an electron, given in Eq. (1.5), but with a,a^ replaced by b,b^. 
Hence it is b^ (b) which raises (lowers) the LL index of a hole and a,a\ which are 
responsible for the LL degeneracy. It follows that a hole with LL index , oscillator 
quantum number ruh and coordinates r/j, has wave function ^m-hnh (r/i) and orbital 
angular momentum projection, (m/j — nh)h. From Eq. ( |1.11[ ), it is also clear that 
the LL index for a hole is associated with the guiding centre of cyclotron motion 
and the oscillator quantum number for a hole is associated with the radius of the 
cyclotron orbit. This is the opposite way round compared to an electron. 
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1.3.2 Single particle problem in the Landau gauge 

The Landau gauge is often used when modelhng rectangular systems of area, 
S = LxLy. The non- interacting Hamiltonian for the electron is 

^EG = -^^^ - ^^T^ - 2i/ieSx— + e^B'^x^ . (1.13) 



2me V dx"^ dy"^ dy 

Clearly \H^Q,Py\ = 0, where Py is the y-component of the canonical momentum, 
p = —iKV. Since py = hky is a good quantum number, we try a solution of the 
form, e^^y^g{x), which gives 



Hkce^'^ygix) = e^'^yy fu^^'' {f^f x + \)9{x). (1.14) 

Here 

are Bose ladder operators, satisfying {fx,f\i\ = ^xx', where X = —kyi'^^. Hence 
the electronic eigenstate is 

^^"^ - — exp (-i^) Cn{x - X), (1.16) 



'Ly V 



where 



Cn{x) = (r\nX = 0}= ^<r| |00> = exp (-^) (pi 

(1.17) 

and Hn is a Hermite polynomial. The eigenstate (p^j^ corresponds to the energy 
hoj^'^ {n + which is (reassuringly!) the same as that obtained using the sym- 
metric gauge. The wave function for the hole can be found similarly and is 

tfi (r) = ;^exp (i^^ Ux - X). (1.18) 



Clearly Eqs. (1.16) and (1.18) describe an extended plane wave in the y-direction 
and a Gaussian peak in the x-direction centred about x = X. Since ky is an integer 
multiple of I'KjLy, X can take the values X = — —^^5 • • • , —^"^-^^j^^^^Lx- The n 
quantum number denotes the LL index as for the symmetric gauge. The angular mo- 
mentum projection, niz, is not a good quantum number in the Landau gauge. The 
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symmetric and Landau gauges are related according to As (r) = Al (r) — VA (r), 
where A (r) = ^Bxy. The unitary operator, which transforms a state in the sym- 
metric gauge into a state in the Landau gauge is Q = exp (— ^A (r)). 



1.4 Graphene 
1.4.1 Background 

Graphene is the name given to a single layer of graphite [sj. Hence it is an atom- 
ically two-dimensional structure consisting of carbon atoms, which are sp^-bonded 
together to form a honeycomb lattice. Historically there has been a lot of interest in 
carbon crystalline lattices, particularly in those which are comprised predominantly 
of benzene rings (with the hydrogen atoms removed), such as graphite, fullerenes |64| 
and carbon nanotubes (65|. Graphene can be viewed as the building block of all 



these structures and as such has been studied theoretically for a long time [66 67 
Many scientists were of the opinion that it could not be synthesised, as it would be 
unstable with respect to scrolling up. However, our current understanding is that 
graphene is most likely created every time someone writes with a graphite pencil. 
Various early works reported the production and observation of extremely thin layers 



of graphite 168 -71 . However, physicists at the University of Manchester, UK, lead 



by K.S. Novoselov and A.K. Geim were the first to isolate true monolayer graphene 



in 2004 and to clearly demonstrate the number of layers in their samples 72 . They 
realised the potential of this result and are largely responsible for making graphene 
one of the most active research areas in condensed matter physics today. For their 
contribution they have recently been awarded the 2010 Nobel prize in physics [73] . 

There are currently several techniques for the fabrication of monolayer graphene. 



The first samples 72 were created using mechanical exfoliation, whereby graphene 



layers are removed from graphite flakes using sticky tape. Other commonly used 



techniques 74 ,75] include epitaxial growth, thermal decomposition of the surface of 
SiC crystals, chemical vapour decomposition, the reduction of graphene oxide and 
unzipping carbon nanotubes. Monolayer graphene can be identified under a care- 
fully tuned optical microscope, when placed upon an appropriate substrate of the 
correct thickness to provide a noticeable contrast. Graphene is the first atomically 
two-dimensional material to have been synthesised and has inspired a search for 
further examples of atomic membranes ]76]. 



According to OPST described in Section 1.1.2 , the electronic states in graphene 
should be localised, provided there is some disorder present, which is always true 
for a real sample. Hence one would expect graphene to be an insulator. Examining 
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(a) 





Figure 1.7: Diagram of the graphene lattice structure (a) and its Brillouin zone 
(b). In (a) the red and blue circles are used to indicate carbon atoms on the A 
and B sublattices. The dashed rhombus indicates a unit cell. In (b) the yellow 
shaded area represents the Brillouin zone with the high symmetry points K ,K' ,M ,T 
indicated. The circles and squares denote positions equivalent to the K and K' 
points respectively. 



the problem from a different perspective, the chemist Linus Pauling, predicted in 
the 1960s that graphene, like benzene, has a resonant valence bond structure and 



hence should be a Mott insulator 77 . More recent theoretical studies of graphite 



find a gap for high enough interaction strengths |78l|79|. A study based on Monte 



Carlo simulations of graphene 80 ,81 , indicates that a gap should be present when 
the surrounding environment has a high enough dielectric constant. However, the 
size of the gap is not exactly determined and may not be large enough to induce 
insulating behaviour. Indeed, numerous experiments have found that graphene is 
in fact an excellent conductor [82| . The cause of the apparent discrepancy between 
theory and experiment is unclear. It may be due to the intermediate strength of e-e 
interactions in graphene, which are large enough to delocalise Anderson electrons, 
but not strong enough to open a gap. In addition, localisation by extended defects 



will not occur due to Klein tunnelling 83 , 84 



1.4.2 Massless Dirac Fermions 

The structure of the graphene lattice is indicated in Fig. |1.7[ a). It is bipartite 
and composed of two triangular sublattices, commonly labelled A and B, so that 
there are two atoms per unit cell. The distance between nearest neighbour atoms is 
a ~ I.42A. The reciprocal lattice of a triangular lattice is also triangular, leading to 
a hexagonal Brillouin zone (BZ), as shown in Fig. fLTjh). Only two of the six vertices 
of the BZ are inequivalent i.e. cannot be joined by a reciprocal lattice vector. These 
are typically labelled the K and K' valleys and behave like a pseudospin degree of 
freedom, so I shall sometimes use the alternative notation fj" and J].. 
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The dispersion relation for graphene was derived by P.R. Wallace in 1947 |66| 
assuming no disorder and no e-e interactions, as outlined in the following. Graphene 
has only one conduction electron per atom, which is in the 2pz state (see the following 
section). Let 0(r) be the 2pz orbital wave function for an electron in an isolated 
atom. Then in the tight binding approximation, the wave function for an electron 
in the bulk material is given by a linear combination of and <^ b ■, with 

= J]exp(27rik.rA)</)(r-r^) , (1.19) 

where the sum over ta is over all points in the A sublattice; is defined anal- 
ogously. Neglecting the overlap of 2pz wave functions centred on different lattice 
sites, it can be shown that the dispersion is E (k) = Haa + |-ffyl_B|; where 



dT^^\H^B (1.20) 



and Haa is given by replacing B with A in Eq. ( 1.20 ). Here N is the number of unit 
cells and H the Hamiltonian. By substituting the expressions for $^ and $b into 
Haa and Hab and considering only nearest neighbour and next nearest neighbour 
hopping, with energies t and t' respectively, it can be shown that the dispersion 
relation measured relative to the energy of a 2pz state in an isolated atom is 

ii;(k) = ±i(3 + /(k))^-t7(k). (1.21) 

The function / (k) is given by 

/ (k) = 2cos ^Vsfcyoj + 4cos (^k^a^ cos (^^^^kyO^ . (1.22) 



Fig. 1.8 shows a plot of the dispersion relation for particular values of t and t' . The 
upper band is called the tt* band and the lower band the vr band. Evaluating / (k) 
at any corner of the BZ, such as K = |^ ^1, and K' = |^ ^1,— yields 
/ (k) = —3. Hence if we include only nearest neighbour hopping {t' = 0), which I 
shall assume from now on, E (k) = at each corner of the BZ. There are two atoms 
and hence two electrons per unit cell, so the vr band is completely full and the tt* 
band completely empty. The bands meet at = 0, so the Fermi energy is Ep = 
and the Fermi surface consists of the six points, which are the BZ corners. The 
touching of conduction and valence bands means graphene may be thought of as a 
zero-gap semiconductor. 
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Figure 1.8: (Produced using source code from the Wolfram Demonstrations 
Project). Plot of the electronic dispersion relation, -E(k), for the case t = 3eV 
and t' = i.e. no next nearest neighbour hopping. 



18 



An approximate form of the dispersion relation at low energies, for example 



close to the K point, may be obtained by substituting k = K + q into Eq. (1.21) 
and assuming qa = |q|a « 1. This gives 

3 

E (q) K +-atq = +hvFq, (1.23) 
where vp = ^ is the Fermi velocity; using t ~ 2.8eV, vp ~ lO^ms^^. The linear 



dispersion relation given by Eq. (1.23) is that of massless relativistic Dirac particles, 
but with the speed of light, c, replaced by vp. Notice that the energy depends only 
on the magnitude of q, so that the dispersion relation at low energies looks like two 



cones with their apexes touching at the K-points, as shown in Fig. 1.8 These are 
refered to as Dirac cones and the K and K' points as Dirac points. 

One normally assumes that electrons are not scattered between the K and 
K' points, due to the large change in momentum that such a scattering would 
require. However, it is possible in the presence of (5-function impurities as we shall 
see in Chapter 7. The single particle Hamiltonian for electrons in the vicinity of 
the K-point was shown using the continuum approximation to be the 2D Dirac 
Hamiltonian [67] : 

Hj^ = vpa ■ n, (1.24) 

where 11 = —iKV — ^ A is the kinematic momentum with q the particle's charge and 
the vector cr contains the Pauli matrices. In the absence of electromagnetic fields, 
the eigenstates with energies +hvpq ar^ 




where S is the area of the graphene sheet, tan^ (q) = qy/qx and H denotes the K 
valley. The appropriate Hamiltonian for electrons in the vicinity of the K'-point can 
be found by changing dy —dy in Eq. (1.24), which results in eigenstates ^'q" (r). 



which can be obtained from ^'q~ (r) by switching the vector components. In both 
^'q~ (r) and ^'q~ (r), the first and second vector components are the amplitudes for 
finding the electron in the A and B sublattices respectively. The conical disper- 
sion relation has been seen explicitly in angle-resolved photoemission spectroscopy 



(ARPES) measurements 85 . In addition, other relativistic behaviour such as Klein 



tunnelling has been observed experimentally for electrons in graphene 183 , 84 



Let us now examine the form of the density of states, p{E). It is shown in 



^Please note the difFerence between q the wave vector, its modulus q = |q| and q the particle's 
charge. 
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Figure 1.9: Density of states, p{E), where the energy is measured in units of the 
nearest neighbour hopping parameter, t. Next nearest neighbour hopping is ne- 
glected {t' = 0). 



Fig. 1.9 for the case t' = 0, where there is nearest neighbour hopping only. An 
analytical solution has been derived for this case |86| . Close to the Dirac points, it 
is easy to show 

P(E) = (1.20) 

The effect of t' is to distort the density of states and move the minimum away 
from ii^ = js] . The effect of an isolated impurity on the density of states has also 



been studied 87 



1.4.3 Graphene chemistry 



A carbon atom has six electrons with orbital configuration Is 2s 2p . The 2s, 2px 
and 2py orbitals hybridise to form three sp"^ in-plane orbitals each containing one 
of the four outer shell electrons. Covalent cr-bonds form between sp'^ orbitals in 
neighbouring carbon atoms and are responsible for graphene's honeycomb lattice 
structure. The remaining electron exists in the 2pz orbital, which is perpendicular 
to the lattice plane. As mentioned in the previous section, the 2pz electrons may 
hop from one atom to another and hence are delocalised. In graphite it is these out 
of plane orbitals which result in the weak coupling between graphene layers. 

The strength of the cr-bonds between the carbon atoms results in graphene 
having some remarkable properties. Its intrinsic strength was measured via nanoin- 
dentation of suspended samples using an atomic force microscope tip and it was 
found to be one of the strongest materials ever tested with a breaking force of 1.8//N 
for a tip of radius 16.5nm |88|. In addition, such strong bonds make it difficult for 
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impurity atoms to replace carbon atoms in the honeycomb lattice, so that graphene 
is naturally a very clean material. The purity or quality of a sample is characterised 
experimentally by the mobility, which is the ratio of the drift speed to the applied 
electric field. Graphene samples have been found to have mobilities 10-100 times 
higher than semiconductors such as silicon. Suspended samples, where there are no 
substrate impurities present, have the highest mobilities. The low impurity density 
allows electrons to travel through the sample in a ballistic manner with a mean free 
path of up to Ifim k lO^a [sj. These impressive electronic and mechanical properties 
make graphene an excellent candidate for a wide range of applications (see Section 



1.4.5). 



1.4.4 Graphene in a magnetic field 

Since its initial isolation, there has been a lot of research on the properties of 
graphene in a perpendicular magnetic field. Interestingly it was predicted that in 
the absence of an applied field, a pseudo-magnetic field may be induced, by straining 



the sample in a particular way, giving rise to gauge fields B ~ lOT [89 90 . Ex- 
perimentally, pseudo-magnetic fields greater than SOOT have been seen in graphene 
nanobubbles |91|. Early on, the anomalous integer quantum Hall effect was ob- 



served, with plateaus detected at filling factors v = +2{\n\ + 1), where n is the 



Landau level index 92,93 . Subsequently, in higher fields, plateaus were seen at 



0, +1, +4, signifying a lift in the fourfold (spin and valley) degeneracy, most 



likely due to the many-body effects 94 . Most recently, fabrication of high quality 
samples has allowed the fractional quantum Hall effect to be observed at u = |, 



indicating again the presence of a strongly correlated electron state [95 96 . Clearly 
many body effects strongly impact the physics of graphene in a high magnetic field. 
However, before turning to the complex many electron problem addressed in this 
thesis, we must examine the single particle problem in a magnetic field. 

Single particle problem for graphene in a magnetic field 

As in the zero field case, we assume the Hamiltonian is block diagonal in the valley 
index: 

( t I ) ■ 

We examine the problem in the symmetric and Landau gauges separately. For 
the symmetric gauge it is convenient to use the points K = ^ (^~^^ ^.nd 

K' = — 1^ ^1, The block for an electron in the K valley can be found by 
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replacing H^. — H^. in Eq. (1.24). This yields 




(1.28) 



where a, a' are defined in Eq. ( 1.6 ) and hwc = h ^^^ is the CRE (between the n = 



and n = 1 LLs) for graphene. The single particle eigenstates of Hd for an electron 
in the K valley are 



<r|n sm) = ^^^^^^(r) = ^'„mff(r)xs = 



/ V'|n|m(r) ^ 

5nV'|n|-l»n(r) 



V 



(1.29) 



where ipnm are the 2DEG wave functions defined in Eq. (1.9), Xs is the spin part 



(s =t, i), dn = 22(''"'0 ^) the normalisation constant and 5^ = sign(n) with Sq = 0. 
The single particle eigenstates, ^'nj|sm(r), for an electron in the K' valley are given 



by swapping the sublattice spinor components in Eq. (|1.29|). The eigenenergies, 
e„ = sign{n)hujc 



|n|, are degenerate with respect to spin, valley and oscillator quan- 
tum number. However, the quadruple degeneracy due to spin and valley pseudospin 
may be lifted in the presence of high magnetic fie^ 



it is clear that Iz , Hd 



7^ 0. Instead we define |97 

h 



ds. From Eqs. ( |1.12[ ) and ( |1.28 ) 

the operator 



Jz 



(1.30) 



where II2 is the 2x2 identity matrix and (x) denotes the Kronecker product. Then 
jz is a generalised angular momentum projection operator, with the second term 
introducing a sublattice component and satisfies j^jz, Hj^ 
wave function ^'m-sm(r), where r =11, J|, jz = \n\ 



m ■ 



0. For an electron with 
2, takes half integer values. 



Using Eq. (1.11) one can show that for graphene 



(nTsm\r'^\nTsm} = (2(|n| + m) + 1 + 5n,o) ^s, 



(1.31) 



so as for the 2DEG, large \n\,m values mean the particle is far from the origin. The 
equivalent problem for a hole can be followed through analogously, beginning with 
the single particle Hamiltonian in the K valley 




(1.32) 
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When solving the problem in the Landau gauge, we use the conventional 
choices of valley, K = |^ (^1, and K' = |^ (^1,—^^. As for the analogous 



problem in the 2DEG, we try a wave function of the form, e* 
electron in the K valley: 



9a{x) 
9b{x) 



for an 



gA(x) 
9b{x) 



fx 




JX 



9a{x) 
9b{x) 



(1.33) 



Hence the wave function for an electron in the K valley is 



(r|n i\ sX') 



lexp 



.Xy 



( SnC\n\-l{x - X) \ 
C\n\{x - X) 






Xs 



(1.34) 



where Cn is defined in Eq. (1.17). It has the same eigenenergy, e^, which was 
calculated for the symmetric gauge. The wave function for a hole in the K valley is 



(r|n i\ sXy 



(1.35) 



/ C\n\{x-X) \ 



l^y \ 



V 



■SnC\n\-lix - X) 






Xs 



I 



The corresponding wave functions in the K' valley, ^'^nisX^'^) ^^"^ *^r4sx('^)' given 
by swapping the sublattice spinor components in Eqs. (1.34) and ( 1.35[ ) respectively. 
It is instructive to reflect upon the differences between the single electron problem 



in the 2DEG and graphene, as emphasised in Fig. 1.10 First note that in graphene, 



unlike for the 2DEG, the LL index n may take negative integer values and these 
correspond to negative energies or hole states in the neutral system. In the 2DEG the 
ground state energy, is positive, whereas in graphene there exists a state with 

energy zero. This state has the unique property that the electronic wave function 
is conflned entirely to one sublattice, which is opposite in the different valleys. In 
the 2DEG, the LLs are equally spaced, whereas in graphene, the spacing between 
neighbouring levels with indices n and n + 1 depends on \/n + 1 — -^/n, which is a 
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decreasing function of n. Thus in graphene, high LLs may become mixed particularly 
in the presence of disorder broadening. The LL energies for graphene ~ s/B (with 
the exception of the n = LL which is independent of B), so that increasing the 
magnetic field strength can help define neighbouring LLs better, assuming that the 
disorder broadening does not grow faster with B. This is one of several reasons for 
requiring a strong magnetic field for the results presented in this thesis. 



1.4.5 Applications 

A large proportion of current graphene research is device motivated. Its remarkable 
electronic and mechanical properties make it an ideal candidate for a huge and 
varied range of applications. I shall summarise some of the most important ones 
here; please see Refs. (OSHTOOl for a more comprehensive overview. 

One aspect, which has received a lot of media coverage, is the use of graphene 
within the electronics industry. Currently, integrated circuits are made with sil- 
icon using complementary metal-oxide-semiconductor technology. However, it is 
well known that in approximately ten years, the reasonable limits of silicon will 
be reached in terms of minituarising devices, increasing their speed and reducing 
the cost. Indeed, there is a physical limit imposed by the fact that silicon loses 
its crystalline structure beyond 10 nm due to increasing thermal fluctuations of the 
atoms. Thus a new material is required for progress to be made and both car- 
bon nanotubes and graphene are being seriously investigated for use in MOSFETs. 
The properties recommending graphene for such a purpose are its high mobility 
electronic transport, two-dimensional nature, strength and high melting tempera- 
ture. Consequently there has been a sustained effort towards the fabrication of high 
performance graphene FETs. Most recently a group at IBM achieved cut-off fre- 
quencies of 100 GHz [101 , which is markedly better than the highest performing 
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silicon MOSFETs. These have a high chance of being used in analog transistors 
for communications technology and imaging. However, the fact that graphene has 
no intrinsic band gap means that it is difficult to turn off^ the electrical current, so 
on-off current ratios are much too low 102] , ruling out their use in digital circuits 
for the moment. Since digital transistors make up the bulk of applications, over- 
coming this problem is now one of the main goals of graphene research. Numerous 
ways of opening a gap have been suggested. There already exists a gap in bilayer 
graphene, which is tunable by a gate voltage 103 . Graphene strongly interacting 
with certain substrates may open a gap, particularly if this induces an asymmetry 



between the two sublattices [104 . A gap can be opened by constricting the path 
of electrons, for example in nanoribbons 105 , quantum dots (106 or nanoconstric- 



tions 107 . Chemically modifying graphene by doping it, with hydrogen 108 109 



or 



fluorine 110 , for example, will also induce a transition to a gapped insulating state. 
We have already seen that an excitonic gap should exist in the regime of strong e-e 
interactions. However, there are problems associated with each of these approaches 
and it is still not possible to control the band gap of graphene to the degree required 
for competitive digital transistors. Indeed, the full extent of graphene's future role 
in the electronics industry remains unclear. 

Graphene is also used in creating nanoelectromechanical systems (NEMs), 
due to its strength and stiffness [111 . For example, drum resonators have been made 
from graphene oxide with high quality factors [112] . It is hoped that graphene will 
be used to realise the ultimate goal of chemical/biological sensors, namely single 
molecule gas detection 113 . The basic idea is that each adsorbed molecule acts 
as an acceptor or donor, causing step-like changes in the resistance. The difficulty 
is that such small signals are easily drowned out by the noise. This is reduced in 
graphene, since it naturally has a low number of impurities. Other useful properties 
in this context are that it remains a conductor even at low carrier densities |92] and 
its 2D nature means its entire volume is available to receive adsorbates. 

Another biological application, which has been recently proposed [114| , is 
using graphene to sequence DNA. How to determine the sequence of nucleotide 
bases on a single molecule of DNA at a reasonable cost and speed, is currently one 



of the most important and challenging problems in biotechnology research 115 



A previous successful technique, known as shotgun sequencing, involves cutting 
up the molecule 116 . In another method, which does not modify the molecule, 
the DNA is drawn through a biological nanopore by an electric field, in an ionic 



solution [117 . The basic idea is that different bases result in different ionic currents. 
Researchers are now trying to use holes in graphene instead of biological pores. 
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Graphene certainly seems to offer several advantages. Its thinness means that only 
one base will be in the hole at any given time allowing single base resolution. It can 
be used both as a membrane and electrodes. It is also stable and the size and shape 
of the hole may be chosen quite precisely. 

Graphene is also being investigated for flexible touch screen technology. The 
material currently used, indium tin oxide, is in short supply. A graphene sheet 



supported by a polymer is used to create the screen's transparent electrodes 118 
The final aim is to use this in ultrathin computing devices, which can be folded 
away and easily transported. 

With much of the basic physics underlying graphene now understood, there 
is a strong emphasis on growing larger, cleaner sheets to be used in technological 
applications. There has been a lot of progress [75] since the original fabrication of 
/im size samples via mechanical exfoliation. The current record is a 30 inch sample, 
which was grown on a flexible copper substrate via chemical vapour decomposition 
using roll-to-roll production, that allows it to be transfered to a target substrate 



of choice 119] . There is also a lot of interest in the nature of the impurities in 



graphene, which limit the mobility. This question is important for improving device 
performance, but also pertains to the fundamental properties of the material. It 
is also relevant to the work presented in Chapter [TJ where we examine collective 
excitations of graphene localised on impurities of different types. 

1.5 Layout of thesis 

The first two Chapters are each self contained. In Chapter [2j I explain the numerical 
techniques required to study the AT in BCC and FCC lattices and give the resulting 
calculated values of the critical parameters. In Chapter [3j a simple model is used 
to study the excitonic AB effect in a quantum ring and determine how it is affected 
by the presence of an in-plane electric field. Chapters |4]j7] are all concerned with 
collective excitations of graphene in a strong perpendicular magnetic field. The 



single particle problem has already been treated in Section 1.4.4 A large amount 
of the theoretical framework and background knowledge required to understand the 
problem is laid out in Chapter [4} The focus is on neutral collective excitations 
in pristine graphene. In Chapter [5j the dynamical symmetries associated with e- 
h complexes, which have both spin and valley pseudospin quantum numbers, are 
investigated. Chapters [6] and [7] are mainly concerned with showing and discussing 
the results. Chapter [6] examines charged collective excitations comprised of three 
particles, which may form at certain filling factors. Chapter [7] deals with the effects 
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of disorder on neutral collective excitations. A short summary of the most important 
findings and possible further avenues of research is given in Chapter [8] 
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Chapter 2 

Critical parameters for the 
Anderson transition in BCC 
and FCC lattices 



The vast majority of Anderson transition (AT) studies use the most elementary 
crystalhne structure, namely the simple cubic (SC) lattice. However, the SC lattice 
does not commonly exist in nature; the only element known to adopt it is the alpha 
phase of polonium [120 . Indeed most metals have body centred (BCC) or face cen- 
tred (FCC) cubic lattices. The aim of this study is to obtain the critical parameters 
and energy-disorder phase diagrams for the AT in BCC and FCC lattices and to 
compare the results with the SC lattice. No change in the critical exponent, v ~ 1.5, 
is expected, since all the systems belong to the same orthogonal universality class. 
We use the transfer-matrix method (TMM) to calculate the electronic localisation 
lengths and finite size scaling (FSS) to analyse the data for different system sizes 
and acquire the critical parameters, defined in the thermodynamic limit. 

Although the localisation properties for an FCC lattice have been studied 



recently for a vibrational problem [121 , to the best of my knowledge no critical 
parameters were reported for an electronic AT in BCC and FCC lattices prior to 
this work [122 . More recently, there has been further interest in modelling more 
physically relevant lattices 123,124 . In Ref. [124 a CsCl-type lattice was studied 
with site selective disorder and they found v k 1.5 — 1.6. This models mixed valence 
semiconductors in which dopants favour the atoms on one of the SC sublattices over 
the other, so that the disorder is not evenly distributed over the lattice as a whole. 
The BCC lattice limit was also taken and agreement with our results obtained. 
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2.1 Numerical approach 



2.1.1 The transfer-matrix method 

The Anderson Hamiltonian for an electron in a lattice with compositional disorder 
is 

Ha = X! £i4ci - 2 *^^^l^r (2-1) 

Here cj(cj) are the creation (annihilation) operators for creating (annihilating) an 
electron in a state localised on an atom at site i. The wave function for such a 
localised state, (r|cj|0), is a Wannier function 
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the Wannier functions form 



an orthonormal set. The two terms in Eq. (2.1) describe competing effects. The 
first term gives the disordered onsite potential energy, which has a localising effect 
on the electron. The onsite energies, Si e [—W/2,W/2], are randomly distributed 
according to a uniform distribution with height W^^. The second term is the 
familiar tight binding term giving the kinetic contribution, which delocalises the 
electron. The tij are the energies required for the electron to hop from site i to site 
j. We include nearest neighbour hopping only and do not include disorder in the 
hopping term. For a ID system, for example, this gives tij = i. All energies 

are then measured in units of the hopping energy, t. 

In the tight binding approximation, the ket describing an electron's state, 
can be expressed as a linear combination of the kets for an electron localised 
on a particular lattice site: 

m = J]Mm- (2.2) 

i 

Substituting lip} into the Schrodinger equation, Hyi|'0) = E\ipy, gives in the ID case, 
the lattice Schrodinger equation for the wave function amplitudes: 

- i'i-i - ipi+i = Etpi, (2.3) 
where I have set t = 1. Rearranging this yields the matrix equation 



(2.4) 



where 

Tj is termed the transfer-matrix since it transfers between the amplitudes ipi^ipi-i 
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and V'j+ijV'i (7, 126| . Eq. (2.4) is solved recursively after choosing initial conditions 
ipojipi- The amplitudes tp]\f,'^N+i are given by 





= T^Tiv-i . . . Ti =T7V . (2.6) 

A correct solution is achieved when the iteration is begun from both ends and 



the wave function and its derivative can be matched at some site in the bulk 126 
However, the amplitudes themselves are not needed to determine the localisation 
length. The expression for the localisation length is found by applying Oseledet's 
multiplicative ergodic theorem [127 128| to the matrices T^y, r^r. 



The salient points of the theorem (for a rigorous mathematical treatment see 
Refs. pTlpSp are: 



i. The following matrix exists: 



hm (r^rT)^-. (2.7) 

Af— >oo 



ii. The eigenvalues of T have the form e'^" < 00. The are termed the Lyapunov 
exponents (LEs) of tat. Note that the matrix r^r, which is a product of the transfer 
matrices defined in Eq. (2.5) is symplecticQ Hence the corresponding eigenvalues of 



r are e-^. 

iii. For each vector, x, in the eigenspace of F, the following limit exists: 

7(x) = lim ^ln||Tivx||, (2.8) 

where 1 1 • 1 1 is the standard Euclidean norm, 
iv. 

7(x) = max({7i}), (2.9) 

where {71} is a subset of the LEs such that the energies e^' have eigenvectors, which 
the vector x has a non-zero projection on. 

Now consider the propagation of a wave through the system. For a ID atomic 



chain it was shown 129 



InT 

7 = -lim— , (2.10) 

Af— >oo iv 

where T is the transmission amplitude at the end of the system. From this it follows 
T/v ~ e^^^ , where T/v is the transmission amplitude at the A^*^ site. This expression 
only makes physical sense for 7 > 0, since more of the wave is reflected and less 



matrix, A, is symplectic if it satisfies AJ = J, wfiere J = ^ ^ 0^ )• 
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transmitted as it propagates through the system. The locahsation length is then 
defined by 

A = -. (2.11) 

7 

It is the exponential decay length of the transmission and thus also the tails of the 
localised wave function. 

In the TMM for a 3D system, the lattice is modelled as a quasi- ID bar 
consisting of N 2D layers, each containing M x M lattice sites, where N » M. 
The localisation length depends on M, the energy E and disorder W. In 3D it 
is scaled according to the system size (width), defining the reduced localisation 
length, Am = X/M. The length of the system depends on the number of transfer 
matrices required to achieve the desired convergence. The transfer-matrix equation, 



Eq. (2.4), still holds, but instead of being scalars, the ipi are vectors of length M x M 
containing the amplitudes for the basis of localised states, ct|0), for lattice sites in 
the i^^ layer. The transfer-matrix in 3D is 

rr _ ( -Cj+\(-^lAf2 - Hi) 



1 



M2 




(2.12) 



Here Om'^ ^^^d ^A/^ are the x zero and identity matrices respectively, Hj is 
the 2D Hamiltonian for the i^^ layer, containing the diagonal disorder and hopping 
between sites in the same layer and Cj is the connectivity matrix describing the 
connections between the i^^ and (i — 1)**^ layers. 

Element Cjk of the connectivity matrix equals 1 if the site j in one slice is 
connected to the site k in the other; otherwise Cj^ = 0. In the case of the SC lattice, 
each site has only one connection to the succeeding (preceding) layer; therefore all 
Cj are unit matrices. For BCC and FCC lattices, the connectivity matrices take 
a more complicated form, but with purely diagonal disorder, they are constant so 
that the inverse C^^ needs to be calculated only once at the beginning of the TMM 
calculations for a given size M. Nevertheless, the additional need to multiply all 
states at each step of the TMM with a dense matrix C^^ reduces the speed of the 
calculation and hence restricts the attainable system sizes. The construction of the 
C^^ matrices is not necessarily always possible for a given lattice along all possible 
lattice vectors. Rather, only selected directions, boundary conditions and M values 
will lead to non-singular Cj matrices. The identification of permissible directions 
for the application of the above TMM requires some care. 

The numerical technique for obtaining the LEs and thus the localisation 
lengths is perhaps different from what one might expect from Oseledet's theorem. 
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Importantly, one does not diagonalise the matrix F. Instead parts iii. and iv. of 
the theorem are utihsed. The procedure is to begin with a set orthonormal vectors, 
{uj}, and arrange them as the columns of a matrix - the M'^ x identity matrix 
for simplicity. The matrix Ti = ri then acts upon 




(2.13) 



where 1^^^ and 2 are the x identity and zero matrices respectively. Choos- 
ing = 0j^2 imposes hard wall boundary conditions at that end of the system. 
Such an action yields a new set of vectors 



(2.14) 



This process is repeated with the vectors from the previous step being left-multiplied 
by Ttv to get the new set of vectors for that set. The idea is that provides the 




largest LE in the sense of Eq. (2.8), U2 provides the next largest LE etc. with 



u^2 giving the smallest positive LE. The negative LEs are not calculated; they are 
simply the negative counterparts to the positive LEs. It is worth noting that the 
eigenvalues still come in pairs, e-'^', in this case, even though the matrix in Eq. 



(2.12) is not generally symplectic. The smallest LE will give the largest localisation 



length, so we define 

A = — (2.15) 

Tmin 

for the 3D case. 

Another important part of the numerical procedure is the reorthonormalisa- 
tion of the vectors, {u^}, at regular intervals. It is the most time consuming part of 
the code, but is required for two reasons. Firstly, after a few transfer-matrix multi- 
plications, the vectors tend to lose their orthogonality. Orthogonality is essential to 
ensure that the vectors all yield different LEs, so that the smallest may definitely 
be found. Secondly, the size of the vectors increases exponentially with the number 
of multiplications, so renormalising them to one periodically, is needed to avoid the 
problem of numerical overflow. The reorthonormalisation is executed after roughly 
every 10 steps. However, the frequency of performing it is varied according to how 
much the norm of the vector leading to the smallest LE changes. The reorthogo- 
nalisation is achieved using the Gram-Schmidt process (see, for example Ref. [130[ ). 
We have {uf} as the set of vectors, which need orthogonalising. Let {vf} be the set 
of orthogonal vectors constructed from the {uf^}, which generate the same vector 
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space as the {uf}- The Gram-Schmidt method sets and for i > 1 

vf = uf- proj^^_^ {uf) =uf- g fS,vf, (2.16) 

where Vi is the vector space generated by [v^Y-Z^^ and proj^ the operator which 
projects a vector onto the vector space V . It is clear that the above definition yields 
orthogonal vectors, which may then be normalised to produce an orthonormal set. 

After each orthogonalisation step (before normalisation), an approximate 
localisation length and its relative error is calculated. If this error is less than some 
set value, the program terminates. After each orthogonalisation step, the vectors 
are placed in descending order of their norms, ensuring that it is the final vector, 
^^^^2, which approximates the eigenvector corresponding to 7min- Suppose after A'' 
transfer-matrix multiplications, there have been no orthonormalisation steps and 
that, for the sake of simplicity, these steps are equally spaced. The estimate for 7min 
is then 

1 "0 

7^,n = ^SHKf/'^°l|. (2.17) 



The relative error in the estimate of 7min, is the absolute error, (57^jjj, which is the 

.if/"«||/Ar}, divided by 7^1,: 



standard deviation of the {ln||u'^^'''"''||/A^}, divided by 7^ 



'mm 
' mm 



]n||,/^/"0| 



„i^/"-o| 
hp I 



2 

--1. (2.18) 



The relative error in the estimate of 7min is equal in magnitude to the relative error 



in the estimate of A and Am, due to the relation Eq. (2.15) 



2.1.2 The lattice structures 

The structure of the BCC lattice is displayed in Fig. |2.1[ a). The construction of the 
TMM quasi- ID bar proceeds along a (100) vector. In this case each site within the 
slice is connected to four sites in the preceding slice and to four sites in the succeeding 
one. There are no connections between sites within the slice, which means that the 
Hamiltonian matrix Hj is diagonal. We use periodic boundary conditions in both 
transversal directions, which results in the connectivity matrix for a slice of M x M 
sites being singular for all even M, thus restricting the system sizes we can use. 



Using a helical boundary condition 131 in one or two directions provides the same 
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Figure 2.1: (a) Three layers of the 3D BCC lattice along a (100) lattice vector. 
Light grey spheres mark the 1®* and 3"^^ layer, dark grey ones indicate the central 
layer. Lines between layers denote the connections between lattice sites. The con- 
nections to the upper-left sphere in the central layer are emphasised by broad lines, 
illustrating its 8 neighbours. The 4 thick light grey lines connect the 1^* and the 
central layer, the black ones go from the central to the 3'^'* layer, (b) Structure of 3 
layers of the 3D FCC lattice along a (111) lattice vector. The broken lines mark the 
cubic unit cell of the lattice. Sites in the I''* and 3'''^ layer are dark grey; sites in the 
central layer are light grey. The thick light grey lines represent connections between 
lattice sites in the same layer for one particular site. The thick black lines represent 
connections between lattice sites in neighbouring layers for another site. The thin 
lines indicate connections to other sites. Some sites in the upper-right corner are 
removed for clarity. 
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singularities and hence offers no advantage. 

Fig. |2.1[ b) shows the structure of the FCC lattice. It proved convenient to 
construct the TMM bar along a (111) vector, so the subsequent layers of the bar are 
close packed. Within the layer, each site has six connections to nearest neighbours. 
In addition there are three connections to the preceding and three connections to 
the succeeding layer. The resulting connectivity matrix can be inverted for each 
size of the M x M TMM slice, but only when we use a mix of periodic boundary 
conditions in one direction and helical boundary conditions in the other. We note 
that it has been shown previously that critical exponents and transition points are 
independent of the boundary conditions 7, 132 . See Appendix [A| for examples of 



the connectivity matrices for system size M = 3. 



2.1.3 Finite size scaling 



As we saw in Section 1.1. 3[ the AT is characterised by a divergent correlation length 



for the infinite system, ^, so that at fixed energy E, S^{W)'x\W — Wc\^'^ and at 
fixed disorder W, ^{E)(x\E — Ec\^'^ , where is the critical exponent and Wc, E^ 
are the critical disorder and energy, respectively, at which the AT occurs |7|. In the 
following discussion, let's assume the case of fixed energy and varying disorder; the 
converse case of fixed disorder and varying energy proceeds analogously. 

In order to extract the critical parameters from the calculated values of 



Am(VF), one applies the FSS procedure outlined in Ref. 133 . The correlation 
length, ^, may be obtained from the localisation lengths for finite system sizes 
A.m{W) by using the one-parameter scaling law Km = f{M/^) [134j . The FSS 
can be performed numerically by minimising the deviations of the data from the 
common scaling curve; this has two branches for extended and localised states. The 
critical parameters are then obtained by fitting the ^ values as obtained from FSS. 
Better numerical accuracy for the FSS procedure can be achieved by fitting directly 
the raw data from TMM calculations using the method applied previously to the 



TMM data for the 3D SC lattice 118 135 . We introduce a set of fit functions which 



include two kinds of corrections to scaling, (i) nonlinearities of the W dependence 
of the scaling variables and (ii) an irrelevant scaling variable which accounts for a 
shift of the point at which the Am{W) curves cross. We use |135| 

Am = fiXrMy^XiMy), (2.19) 

where Xr and Xi are the relevant and irrelevant scaling variables respectively. Note 
that y < 0, so that the contribution from the irrelevant scaling variable is reduced 
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for larger system sizes. The function Am{W) is then Taylor expanded 



Am = X! xrM"^/n(XrMV-), (2.20) 

fn = J]a^kx'^M'^/r (2.21) 

Nonlinearities are taken into account by expanding Xr and Xi in terms ofw = {Wc — W) /Wc 
up to order nir and mi, respectively, 

TTtr m; 

Xr{w) = ^-^"^ XiH = X! ^rnW"'^ (2-22) 

with 61 = Co = 1. Note that this FSS procedure assures the divergence of ^ and 
hence it is not the divergence itself but rather the quality of how the model fits 
the computed reduced localisation lengths Am which determines the validity of the 
scaling hypothesis. 

2.1.4 Nonlinear fitting 

The nonlinear fit was carried out using the built-in Mathematica function 



"NonLinear Regress" . It uses the Levenberg-Marquardt method 135 1361 to find the 



parameters {aij, bi, Ci, y, i^, Wc), which provide the best fit of the data ({M, W, Am}) 



to the model / given in Eqs. ( 2.19 )-( 2.22) with fixed values of rij-, ni, rrij-, mi. The 



data points are weighted according to ,^ n , m ^■? > which is the inverse variance of 
Am- Importantly, it also provides errors for Wc and v and statistical measurements 
of the goodness of fit, which indicates the validity of the model. This information 
enables us to choose the most appropriate values of rir, n\, mr, m\. The goodness 
of fit should be maximised and the errors minimised whilst simultaneously keeping 
the number of parameters low. The number of parameters (not including Wc and 
v) is {n\ + 1) (nr + 1) + m-i + mr + 2 for the case ni > and nr + m-r + 1 when ni = 0. 



One of the outputs of "NonLinearRegress" is x 137 . It is related to the 
issue of how to quantitatively judge how well a set of parameters fit the data for 
a particular model. This is done using the concept of maximum likelihood estima- 
tors. One calculates the probability that the data set could have occured, given a 
particular set of parameters. Suppose the data is {{2:1, yi}, {x2, 2/2}, • • • , {^^n, yn}} 
and the goodness of fit of the model y{x) is to be tested. For data which may 
take continuous values, such as our localisation lengths, the aforementioned prob- 
ability is zero. Instead we calculate the probability that a data set contained 
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in {{xi, yi + Ay}, {x2, 2/2 + Ay}, . . . , {xn, Hn + Ay}} could have occured, given the 
model y{x), where Ay « 1 is fixed. This is 



Pec Y\ exp 



Vi-y {xi) 



Ay. 



(2.23) 



The expression for P assumes that each data point, yj, is normally distributed 
about y{xi) with standard deviation ai and that the errors are independent. The 
assumption of a normal distribution should hold due to the central limit theorem 
[T38J. However, one should bare in mind that this is not always the case in practice. 



Maximising Eq. (2.23) is equivalent to minimising 

y{xi) 



X 



Vi 



(2.24) 



The method of minimising to find the best-fit parameters, is refered to as a least- 
squares fit. The Levenberg-Marquardt method 135 ,136 minimises iteratively. 



An initial guess for the parameters is used to calculate an initial value for x^- The 
parameters are then altered to further lower x^- This is done using either the inverse 
Hessian method or method of steepest descent, depending on how well the function 
X^(a) (a is a vector containing the parameters) approximates the behaviour of 
close to its minimum value. 

To examine the probability distribution of different values of we must 
define the incomplete gamma function: 



where 



Tp{a,x) 



T{a) 



1 

— dtt-\ 
r(a) Jo 



(2.25) 



f 

Jo 



dt t 



(2.26) 



is the ordinary Gamma function. Let /i be the number of degrees of freedom: 



// = #data points — ^parameters. 



(2.27) 



Then rp(^, gives the probability that the observed chi-squared value for a correct 
model is less than x^. It is usually the complement of Tp, i.e. Fq = 1 



Fp, which 



is considered. The closer rq(^, ^y) is to 1, the better the model. However, values 



rq(^, ^) K 1 are suspicious and may mean that the errors have been overestimated. 
A reasonable fit is thought to have x^ ~ //. 
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Figure 2.2: Phase diagram for the BCC lattice. The cyan region represents the 
approximate location of the phase boundary. Its edges (the solid black lines) were 
determined by comparing localisation lengths with errors ^ 10% for system sizes 
M = 7 and M = 9 in the {E, W) plane. The solid red squares (■) are points calcu- 
lated by performing high-precision FSS on localisation data with an error ^ 0.1%. 
The hollow squares (□) are reflections of the solid squares in the E = axis. The 
diamonds (♦) denote the band edges at W = 0. They have been joined to the phase 
boundary edges calculated for higher disorders as a guide to the eye. The dashed 
lines are the theoretical band edges +{Z + W/2), where Z is the coordination num- 
ber. The horizontal dotted line is the BCC estimate 21.13 for Wc of Ref. 139 . The 



light grey, shaded area in the centre contains extended states; states outside the 
phase boundary are localised. Error bars are within symbol size for ■, □. 

2.2 Calculations and results 



2.2.1 Phase diagrams 



Fig. |2.2| and Fig. 2.3 show the phase diagrams for the BCC and FCC lattices, 
respectively. Originally a grid of W versus E values was created with separation 
AE = AW = 0.5t. At each point the nature of the electronic wave function was 
determined by comparing the reduced localisation lengths. Am, calculated for system 
sizes M = 7 and M = 9 with error ^ 10%. If Ag > (<)A7 at the same values of 
E and W then the point {E, W) in the phase diagram is identified as extended 
(localised) . The edges of the phase boundary were obtained by averaging separately 
over the three extended and localised points {E, W) nearest to the boundary and 
then connecting such averages using a spline fit. Data points are not obtained for 
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Figure 2.3: Phase diagram for the FCC lattice. Symbols, lines and shaded areas 



have the same meaning as in Fig. 2.2 with the diamonds (♦) representing the band 
edges —12 and 4 at zero disorder [120 and the dotted line corresponding to the 
self-consistent estimate Wc = 33.08 139 . Error bars are within symbol size. 



lower disorder values, as the fluctuations in the LEs, due to the small system sizes, 
become too big; higher values of disorder smooth out these fluctuations. 

A striking difference between the phase diagrams is that for the BCC lattice 
the phase boundary is symmetric about the line E = 0, whereas for the FCC lattice 
it is not. This is due to the bipartiteness of the BCC lattice which consists of two 
SC sublattices, one displaced half the distance along a body diagonal of the other. 
Hence for any site in one sublattice, its nearest neighbours are in the other sublattice. 
Such connections result in states coupled by a bipartite symmetry transformation — 
which is exact for the case of no diagonal disorder — with eigenenergies of the same 
magnitude but opposite sign having approximately the same localisation lengths; 
this produces a symmetric phase diagram. The FCC lattice is non-bipartite, so 
such a symmetry in its phase diagram is not observed. 



In Figs. 2.2 and|2.3[ we also indicate via horizontal lines previous results for 



critical disorder strengths |139|, which were based on the self-consistent theory of 
localisation and obtained for a momentum cut-off of 2pF, where pp is the Fermi 
momentum. These estimates agree well with our results for both BCC and FCC 
lattices. The SC result Wc = 13.91 of Ref. 139 , however, deviates more strongly 



from recent Wc = 16.54 + 0.02 [l8||32}[T35] estimates. 
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Figure 2.4: Reduced localisation lengths Km versus disorder W for BCC lattice. 
System sizes M are 3( ), 5(n), . . . , 15( ). Error bars are within symbol size. Lines 
are fits to the data given by Eqs. (2.19) - (2.22) with = 3, ni = 2, mr = 3, mi = 1. 



2.2.2 Critical parameters at = 

The TMM calculations were performed for system sizes up to M = 15. In order 
to examine the localisation properties at the band centre for the BCC lattice and 
the barycentre (average energy eigenvalue) for the FCC lattice, we set S = in 



Eq. (2.12). A value of the critical disorder Wc was approximated using the phase 



diagrams described above and then localisation lengths A were calculated for a range 
of W close to this approximate value with accuracy ranging from 0.1% for small 
system sizes M to about 0.14% for the largest. Note that the term critical disorder 
is used to indicate that there are no further extended states at E = for disorders 
W > Wc, extended states may still exist for W > Wc at other energies E, as shown 



in Fig. 2.3 



The reduced localisation lengths for the BCC lattice are displayed in Fig. 
2.4, Since Am increases (decreases) as M increases for extended (localised) states, 



the Am curves for different system sizes should cross at Wc- Note how the crossing 



point of the curves in Fig. 2.4 shifts with changing M. In most cases this indicates 
the need for an irrelevant scaling variable introduced via non zero values of rii and 



TTii in Eqs. (2.20) and (2.22). Fig. 2.5 shows the results of the scaling procedure 
for Uj- = 3, ni = 2,mi. = 3, mi = 1. The scaling curve exhibits localised (bottom) 
and extended (top) branches as expected for the MIT. Divergence of the scaling 
parameter at W ~ 20.75 indicates the critical value of the disorder. Table [2rT| a) 
gives some examples of models providing the best fits and the resulting critical 
parameters. The values of the critical disorder and critical exponent are obtained 
by averaging over all the best fit models. 
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Table 2.1: Critical parameters for the MIT in the BCC lattice. All errors quoted 
are standard errors, (a) 3 examples of FSS results with varying Tir, nj, rrir, mj at 
fixed energy E = 0. We use 91 data points, equally spaced in the indicated intervals 
(cp. Fig. 2.4), for each set of nr,ni,mr,mi. Varying rir, ni, mr, mi, we obtain 41 
best fit models in order to produce the indicated averages, (b) Similar FSS results 
obtained for 3 out of 14 best fit models from 82 non-equally spaced data points 
at fixed W = 15 for the indicated energy intervals. Some of the best fit models 
use irrelevant scaling, although the examples explicitly detailed do not, so that 
rrii = Hi = and the parameter y is not used, (c) Results at W = 17.5 (cp. Fig. 2.8) 
for 3 out of 8 best fit models with 108 non-equally spaced data points used in each 
FSS procedure. All the best fit models use no irrelevant scaling. The numerical 
fitting procedure continued in all cases until convergence was reached or (a,b) 5000, 
(c) 1000 iterations had been completed. When averaging, non-converged results 
were neglected. 



(a) 



AM E AW riy- n\ rrij- nii Wc v \y\ 

3 - 15 20.3 - 21.5 2 1 20.95(1) 1.67(5) ~- 

3 - 15 20.3 - 21.5 3 1 1 4 20.92(2) 1.51(9) 1.7(5) 

3 - 15 20.3 - 21.5 3 2 3 1 20.75(3) 1.70(9) 3.0(5) 



averages: 20.81(1) 1.60(2) 



(b) 



AM 


AE 


W 




mr 




V 


9 - 13 


9.9 - 10.9 


15 


2 


1 


10.38(1) 


1.32(5) 


9 - 13 


9.9 - 10.9 


15 


2 


2 


10.38(1) 


1.22(5) 


9 - 13 


9.9 - 10.9 


15 


3 


4 


10.40(1) 


1.03(3) 


averages: 










10.39(1) 


1.21(2) 


(c) 


AM 


AE 


W 


rir 


rrir 




V 


7- 15 


10.5 - 11.5 


17.5 


2 


1 


10.98(1) 


1.55(6) 


7- 15 


10.5 - 11.5 


17.5 


3 


2 


10.99(1) 


1.48(6) 


7- 15 


10.5 - 11.5 


17.5 


3 


4 


10.99(1) 


1.36(7) 


averages: 










10.99(1) 


1.45(3) 
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Figure 2.5: Scaling function (solid line) and scaled data points for the BCC lattice 
and nr = 3, ni = 2, mr = 3, mi = 1. Symbols denote the same values of M as in Fig. 



2.4 Inset: Dependence of the scaling parameter ^ on the disorder strength W for 



the 13 W values shown in Fig. 2.4 The dashed vertical line indicates the divergence 
of ^ at = 20.75. In all cases error bars are within symbol size. 



Let me emphasise that given the precision of our data, simply presenting the 
data with minimal value for a given set of rir, rrir, n\,m\, systematically underesti- 
mates the true error. The reason is twofold. Firstly, there is no a priori justification 
for which set of n^, m^, ni, rrii values to use. Also the variation in results for different 
n-f:,m^,n\^m\ is typically larger than the error bars suggested for each individual 
set. Secondly, one can change the range of system sizes and disorder /energy values 
for the localisation data to be included in the fit. This again leads to changes in 
the critical parameters which are usually beyond the error bars generated in each 
individual nj., mr, ni, mi fit. Hence, in the absence of a clear criterion for which of 
these fits to choose, the strategy is to (i) delete all obviously erroneous fits, i.e. those 
which do not converge or which converge to unphysical values; (ii) to average over 
the remaining results with a proper estimation of accumulated error, based on the 
individual errors for each n^,m^,n\,m\ choice. Last, all this should be done while 
keeping the number of parameters — as determined by n^, rrij-, ni, mj — as small as 
possible. The accumulated standard error is given by 



(2.28) 



where n is the number of models being averaged over and ai is the standard deviation 
for the i*^ model. Note that in ultrahigh precision studies 34 , the issues mentioned 



above are less of a problem, so such an averaging procedure is not used. 



42 



26 



26.5 27 
W 



27.5 



Figure 2.6: Reduced localisation lengths Km versus disorder W for FCC lattice. 
Symbol sizes M are 3( ),4( ),...,15( ). Error bars are within symbol size. Lines 
are fits to the data given by Eqs. (2.19) - (2.22) with = 2, mr = 2. Lines for even 
M have been removed for clarity. 



Table 2.2: Critical parameters for the MIT in the FCC lattice at = (cp. Fig. 



2.6). We use 91 data points for each FSS and the obtained 31 best fit models 



average as indicated. The fitting procedure was continued until convergence was 
reached or until 5000 iterations had been completed, although only models for which 
convergence was reached were included in the averaging process. No irrelevant 
scaling was necessary, so = m\ = 0. 



AM 


E 






mj. 






3 - 15 





26 - 27.5 


1 


2 


26.73(1) 


1.58(2) 


3 - 15 





26 - 27.5 


2 


2 


26.73(1) 


1.58(3) 


3 - 15 





26 - 27.5 


3 


3 


26.73(1) 


1.67(5) 


averages: 










26.73(1) 


1.60(1) 



Results of the TMM calculations for the FCC lattice are shown in Fig. 2.6 
In this case the lines for constant M cross at the same point — at least within the 
accuracy of the calculated Am — indicating that the use of the irrelevant variables 
in Eqs. (2.20) and (2.22) is not necessary in most cases and n\ = m 



of the fit for ?7-r = 2, mr 



0. Results 

2 are displayed in Fig. 2.7 The transition at ~ 26.73 

' ' as 



is clearly indicated. More examples of best fit models can be found in Table 2.2 
well as the average values of the critical parameters. 



2.2.3 Critical parameters away from the band centre 

We also perform calculations where we fix the disorder and allow the energy to vary 
across a critical value for the transition. It is known that such investigations 
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Figure 2.7: Scaling function (solid line) and scaled data points for the FCC lattice 
and TT-r = 2,mr = 2. Symbols denote the same values of M as in Fig. 2.6 Inset: 



Dependence of the scaling parameter ^ on the disorder strength W for the 13 W 



values shown in Fig. 2.6 In all cases error bars are within symbol size. 



are numerically more difficult due to the influence of density of states effects 140 
Results of the TMM and FSS calculations for the BCC lattice with W = 17.5 can 



be seen in Fig. 2.8 The poorer quality of the fit compared to the calculations where 
the energy was fixed at zero is evident. This is not due to using lower accuracy 
data, as the maximum raw-data error remained at 0.1%. Hence we attribute it to 
complications arising from a varying density of states close to Ec at the attainable 



values of M. Results for the critical parameters are shown in Table 2.1 (b) and (c) 
for 1^ = 15 and W = 17.5 respectively. The low value of for the case W = 15 can 
be attributed to the use of fewer data points in the FSS and only using three values 
of M. We note that this is consistent with the lower values of u obtained in the 



diagonalisation studies 133 as mentioned in the Section 1.1.3 It appears that the 



FSS procedure systematically reduces the values of the critical exponent for data 
from smaller systems or of lower accuracy. 

Results for the TMM and FSS calculations for the FCC lattice with W = 18 



can be seen in Fig. 2.9 Table 2.3 gives examples of the best fit models and shows the 



resulting average critical parameters. Note that both estimates of v in Table 2.2 and 



Table [2^ are consistent with the result 1.57(2) for the SC lattice [18] . The TMM 
data for this system has the highest quality out of all the data presented in this 
study. Thus for this case only, we include the values and Fq (///2, x^/2) values to 
demonstrate the goodness of fit. For all the models ~ which indicates a good 



fit. We remark that the results presented in Tables 2.1 , 2.2 have lower precision, but 
this is clear since the accuracy of the individual data points is less, due to smaller 
system sizes and fewer available data points very close to the transition. 
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2 4 6 

log|g(^/iW) E 

Figure 2.8: Localisation data for the BCC lattice with W = 17.5. System sizes M 



are 7(a), 9( ),..., 15( ) as in Fig. 2.4 Error bars are within symbol size. Left: 
Scaling function (solid line) and scaled data points using = 2,mr = 1. Right: 
Reduced localisation lengths Am versus energy E. Lines are fits to the data given 



by Eqs. (2.19) - (2.22) with = 2,mr = 1. 



- 0.65 




Figure 2.9: 
9(+),ll(0), 



Data for the FCC lattice with W = 18. System sizes M are 
..,15(- ) as in Fig. 2.7 Error bars are within symbol size. Left: 
Scaling function (solid line) and scaled data points using = l,mr = 2. Right: 
Reduced localisation lengths Am versus energy E. Lines are fits to the data given 
by Eqs. ^IM - (IZ^ with = 1, mr = 2. 
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Table 2.3: Critical parameters for the MIT in the FCC lattice at TV = 18 (cp. 
Fig. 2.9). We use 83 data points to perform the FSS and use 8 best fit models to 
obtain the averages. None of the best fit models use irrelevant scaling. The fitting 
procedure was continued until convergence was reached or until 1000 iterations had 
been completed; only models for which the method converged were used to obtain 
the averages, n is used to denote the number of degrees of freedom. 



AM 


AE 


W 






Ec 


V 


x' 






9 - 15 


8.52 - 8.88 


18 


1 


2 


8.683(3) 


1.63(5) 


79.0 


0.44 


78 


9 - 15 


8.52 - 8.88 


18 


2 


2 


8.687(4) 


1.65(5) 


78.0 


0.44 


77 


9 - 15 


8.52 - 8.88 


18 


3 


1 


8.685(3) 


1.62(5) 


77.6 


0.46 


77 


avera^ 


;es: 








8.684(2) 


1.63(2) 









2.3 Conclusions 

The transfer-matrix method and FSS were used to determine the critical parameters 
of the Anderson transition for the BCC and FCC lattices. The values of the critical 
exponent v are in good agreement with the results obtained previously for other 
systems belonging to the Gaussian orthogonal universality class. The increase of 
the critical disorder from 16.54 for the SC lattice to 20.81 for BCC and 26.73 for 
FCC lattice may be attributed to an increasing number of nearest neighbours which 
for the above structures equals 6, 8 and 12, respectively. More nearest neighbours 
connected to a given site provide more paths for electronic transport, so stronger 
disorder is needed to localise eigenstates of the system. The universal localisation 
properties of a 3D system and the presence of an MIT are however not affected 
in accordance to the scaling theory of localisation |13j and are in agreement with 
results [141 showing that they depend only on the dimensionality of the system, 



but not on the number of nearest neighbours in the lattice. 

Our results and their interpretation are consistent with investigations of clas- 



sical bond and site percolation models on SC, BCC and FCC lattices. In Ref. 142) 
it was found that the percolation thresholds for these lattices decrease with increas- 
ing number of nearest neighbours; more neighbours allow for easier formation of a 
percolating cluster, or, as in our case, the formation of extended states. 
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Chapter 3 

Excitonic Aharonov-Bohm 
effect with external fields 



3.1 Introduction 



The manipulation of light and excitons is an area which has sparked much recent 
interest 143 -147 . The speed of light has been slowed down to an amazing 17ms^^ 



in an ultracold gas of sodium atoms 147 by the use of electromagnetically induced 



transparency techniques 1481. Bose- Einstein condensation of an excitonic gas is 



a phenomenon considered theoretically a long time ago 149 . However, it is only 



recently that the ability to grow high quality layered semiconductor structures has 
allowed for a successful experimental investigation of excitonic superfluids. Spatial 
patterns in photoluminescence measurements of an exciton gas in GaAs/AlGaAs 



coupled quantum wells have been observed 1501 and are thought to be a signature 



of quantum degeneracy [151 . Zero Hall voltages measured in bilayer quantum 
well systems for particular layer separations and magnitudes of a magnetic field 
applied perpendicular to the layers, are also understood to indicate the presence of 
an excitonic condensate [152 . There has also been experimental evidence for the 



formation of a polariton BEC 1153 and a room temperature polariton laser 154 



Such discoveries are very important for quantum computing applications 
155 . In current computers, electrons are used for information processing and pho- 



tons for communication. Conversion between the two media places limitations upon 



the machine's efficiency. The creation of an exciton-based integrated circuit 156 



could provide an exciting solution to this problem. The possibility of using excitons 
for data storage has also been investigated using indirect excitons in self-assembled 
quantum dots [144| and coupled quantum wells 143 . Here the electron and hole 
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are separated using a gate voltage, significantly prolonging their lifetimes. 

In this Chapter we study the Aharonov-Bohm effect (ABE) for an exciton in 
a nanoring with an in-plane electric field. Oscillations in the physical properties of 
the exciton as a function of magnetic field are signatures of the quantum mechan- 
ical interference between the electron and hole wave functions. A mechanism for 
controlling the excitonic lifetime by tuning the electric and magnetic fields is discov- 
ered. The excitonic Aharonov-Bohm effect (XABE) is particularly interesting; as 
a neutral particle, one would not expect an exciton to couple to a magnetic vector 
potential. However, an exciton's finite size and the internal motion of the electron 
and hole mean that the XABE is possible. AB oscillations have been predicted for 
optically created excitons in the ring geometry by a variety of different models. Ini- 
tial treatments used a simple ID continuous model that could be solved analytically 
with a perpendicular solenoidal magnetic field through the centre of the ring and 



a contact electron-hole interaction potential [157 -159 . AB oscillations were found 
both in the excitonic ground state (GS) and in the corresponding oscillator strength 
(OS) for rings with diameters smaller than or comparable to the excitonic Bohr 
radius, qb = h^e/ixe^ (e is dielectric constant and /x reduced mass) 



158 . 



Other models tried to capture the essence of what can currently be achieved 
experimentally. An important property that has been widely treated is the finite ring 
width [160 -167 . In a 2D Hubbard model approach with a short range electron-hole 
interaction, excitonic AB oscillations in the GS energy were seen to survive for finite 



ring widths [163 . In contrast, a continuous approach introducing the ring geometry 
by a confining potential and using a screened Coulombic interaction found no GS AB 
oscillations for any finite ring widths [160 168 . However, oscillations in other low 



lying energy states were detected 160[ for rings of small enough diameter and width. 
Eccentric ring geometries have also been examined for different confining potentials 



167 . Although the eccentricity reduces the amplitude of oscillations, it also causes 



previously dark states to become optically active, so a slightly eccentric ring could be 
a good candidate for the experimental observation of the XABE. Another approach 
thought to enhance the effect is the polarisation of the exciton [164 165, 169 . This 
occurs naturally for particular choices of semiconductor materials. For example, in 
self assembled InAs rings on a GaAs substrate, the electron is most likely to be 
found towards the centre of the ring, whereas the hole is largely confined to the 
outer radius of the ring [170 . 



The role of the electron- hole interaction has also been investigated [159 ,164 



the oscillation amplitude was found to decrease as the Coulombic interaction be- 



comes less shielded, in agreement with intuition. The effects of impurities 164 ,171] 
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and more general disorder [172 173 have also been studied. The introduction of an 
isolated impurity potential breaks the rotational symmetry of the ring, resulting in 
the coupling of excitonic states with different angular momenta and anticrossings 
appearing in the energy levels. The absorption coefficient is also modulated, since 
scattering by impurities causes optical emission from previously optically inactive 
states. The introduction of disorder to the onsite potential energies has been seen 
to suppress the ABE for certain interaction strengths, when many realisations of 
disorder are averaged over [173 . However, the shift in linear absorption peak as a 
function of magnetic flux is dependent on the realisation of disorder. For certain 
realisations and interaction strengths, the disorder actually enhances the effect. In 
view of studies conducted on attractive two particle systems in a disordered ring, 
which show a non monotonic dependence of the localisation length upon the inter- 
action strength [174 , it is probable that disorder enhances the ABE for particular 
interaction strengths even when averaging over disorder realisations is undertaken. 



It is only very recently (after the publication [175] of the work described 
in this Chapter) that the XABE has been observed experimentally for the ring 
geometry [iTGpTT] . InRef. (l76] the energy peak in the photoluminescence spectrum 
of an ensemble of type-I InAs/GaAs quantum rings was seen to oscillate as a function 
of magnetic field. In Ref. [177] oscillations in the photoluminescence peak of a single 
InAs/GaAs quantum ring were observed. It should be noted that such experiments 
do not strictly meet the criteria set out in Section 1.2 for the observation of the ABE, 
since the magnetic field penetrates the ring and so the particles are not confined to 
a magnetic field free region. However, the observed effect can be modelled as if the 
magnetic field was concentrated in the middle of the ring, so is considered to be an 
ABE. Prior to this, AB oscillations were observed for negatively charged excitons 



or trions (a bound structure of two electrons and a hole) [178 . AB oscillations were 
also seen for excitons confined to type-II quantum dots [179 180 . In this case, only 
one of the particles (either the electron or the hole depending on the material) exists 
in a ring geometry created by the circumference of the quantum dot. The other is 
localised in the centre of the quantum dot. Hence these AB oscillations are really 
due to the motion of a single particle. 

With quite a large body of theoretical work already existing on the XABE, 
the aim of the project described below was to see how the application of an in- 
plane electric field might affect the amplitude of AB oscillations in the exciton's GS 
energy and OS. It was previously suggested that the electric field would enhance 
the XABE 181 ; our results suggest this is only partially true. Interestingly they 
indicate the vanishing of the OS for particular values of electric and magnetic field. 
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Figure 3.1: The geometry of a ring in the presence of electric and magnetic fields. 



pointing to the possibility of trapping light. This result may lead to applications in 
data storage. 

3.2 Theoretical approach 
3.2.1 Single particle problem 

First consider the case for a single charged particle, since it is the wave functions of 
the electron and hole, which provide a basis for the desired excitonic wave function. 
The standard Hamiltonian for a particle with charge q and mass mo in the presence 
of an electric potential (p and a magnetic vector potential A is 

H = ^ ^ + q<p. (3.1) 



The ring geometry is illustrated in Fig. 3.1 The ring is modelled as a ID circle 



of radius p. It is threaded by a magnetic flux perpendicular to the plane of the 
ring and is in the presence of a constant in-plane electric field in the x— direction of 
magnitude Uq. Let Be and ip^*^^ denote the azimuthal coordinate and wave function 



of the electron respectively. Then using the Hamiltonian in Eq. (3.1) with q = — e, 
the Schrodinger equation for the electron may be written in terms of dimensionless 
quantities as 

( - 2i$ J- - $2 + uocos {Be) + A(^)) ^p('^ (9^) = 0, (3.2) 

where eo = 2m~p^' ^ ~ ~^Ic (iiumber of flux quanta), uq = and A^^^^ = 
where is the energy of the electron. In the absence of the electric field (mq = 0), 
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Eq. (3.2) has translational symmetry and the solutions are plane waves [158 : 



2tt 



= {N- 



(3.3) 



For the hole: 



27r 



(3.4) 



where A^, N' are integers and I have assumed the electron and hole masses are the 
same. 

The electric field breaks this symmetry and we use a Fourier Ansatz for the 
electronic wave function 



00 

^('^)(^,) = -= 2 c(^)e'^'^ 



(3.5) 



By substituting this into Eq. (3.2) and using the orthogonality of the complex ex- 



ponential functions, one may show that for any integer m, 



X{e)c(e) = _ (e) _ ^ (e) _ ^ (e) 



(3.6) 



Eq. (3.6) can be expressed as an eigenvalue problem: 



(3.7) 



where 



M 



UQ 

2 



2 



-1 - $)^ 



UQ 

2 



up 
2 



up 
2 



'op 
2 



Up 

2 



V 

(3.8) 

and c^*^) = ^. . . , c^^, . . . , c^j, Cq^'', c^"^-*, . . . , Cm'', ■ ■ Truncating the summation in 
Eq. (3.5) at some finite integer Mmax allows Eq. (3.7) to be solved numerically to ob- 



tain the dimensionless electronic energies, and the coefficients c^*^). These can be 
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Figure 3.2: The two lowest single particle energy states as a function of magnetic 
flux <I> at various values of electric field no with Mmax = 30. 



substituted into Eq. (3.5 ) to approximate the electronic wave functions (^e)- The 
hole wave functions 'ip^ (Oh) and energies A^,-' are obtained by solving an eigenequa- 
tion similar to Eq. (3.7) with <I> and uq replaced by — <I> and —uq respectively. 



Fig. 3.2 shows the single particle energies plotted as a function of magnetic 
flux <I> for the case N, N' = 0, 1. The AB oscillations are clearly visible. The energy 
values are the same for the electron and hole. This can be seen by studying the 
); it is clearly similar to a matrix with $ and uo replaced by — $ 
-Uq. For the case uq = we obtain quadratic behaviour in agreement with 



matrix in Eq. (3. 
and 



Eqs. (3.3) and (3.4). 



Fig. 3.3 shows the single particle energies plotted as a function of electric field 
Uo for various values of ^ and N, N' with ^ N,N' ^3. The N, N' value associated 
with a particular curve can be deduced by noting its value at <I> = 0, uo = and 



comparing this with Eqs. (3.3) and (3.4). All the curves are symmetric about the 
line no = due to the symmetry of the problem. This symmetry can also be seen 



in the matrix in Eq. (3.8), which is similar to a matrix equal to it, aside from a 



sign reversal in the off-diagonal entries. For the case $ = 0, Eq. (3.2) reduces to 



the Mathieu equation fl82| and the single particle energies are proportional to the 
characteristic Mathieu values. 



3.2.2 Self-consistent solution for the exciton 

We are now in a position to solve the problem of an exciton in a nanoring with 
external fields. The excitonic wave function can be expanded in a basis of electron 
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Figure 3.3: Single particle energies as a function of external field uq at various 
values of magnetic flux $ with Mmax = 30. 



and hole wave functions: 



(3.9) 



N,N' 



The Hamiltonian is given hy H = He + + V, where He is the Hamiltonian for 
the electron, Hh is the Hamiltonian for the hole and V the interaction term, which 
depends on the distance R between the electron and hole. We assume a short range 
interaction 



v[R{9e-eh)] 



eo 



2-1^05 {Be 



(3.10) 



where vq is the average interaction strength in units of eg. Values for vq are chosen as 
integer multiples of — l/vr^, to ensure that the inverse decay length for the excitonic 
wave function, 7 = 7r|wo|eo-, 
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, is comparable to the inverse of the ring 
circumference. The Schrodinger equation for the exciton may now be written as 

Ann' {\%^ + - a) V'^^ {Oe) V-Jf.^ {Oh) + 2t:vo5 {Oe - Oh) ^ (^e, Oh) = 0, 

N,N' 

(3.11) 

where A is the excitonic energy in units of eo. Multiplying by j^^^'' {Be) ip^^ (Oh) 
for particular N,N' e 7^ and integrating over 9e,9h ^ [0,27r] allows us to obtain an 



expression for the coefficients in Eq. (3.9): 

2ttvo 



A 



NN' 



+ a!v^^ - A 



Jo 



(3.12) 
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Using Eq. (3.9) with Of, = 9h = 0, we may write 



JO ^ jy; Jo 



Let 



and 



/*27r A A 

Jo 



(3.13) 
(3.14) 



Ittvq 



NN'No% 



A^^) + A^^) - A Jo 



deV'J^) (e) {6) (0) (0) . (3.15) 



Then Eq. (3.13) can be expressed as 



G 



2 Gnn'Pnn'No% 



(A) 



(3.16) 



N,N' 



In order to find approximate solutions to (3.16), we cut-off the sums at a maximum 
value A^max for N, N' , Nq, Nq. We have chosen A^max = 15 for calculation of the 
excitonic energies and A'^max = 11 for the oscillator strengths; we have tested that 
the results do not change appreciably for the range of <I>, uq and vq considered 
here. Mapping N,N' ^ K and No,%^ Kq according to = (iV - 1) iY^ax + N' , 



allows Eq. (3.16) to be reformulated in 2D: 



Gko=Y^GkPkkAA) 



(3.17) 



K 



The excitonic energies may now be calculated numerically by determining the val- 
ues of A which result in the matrix P^Kq having an eigenvalue equal to 1. Note 
that it is the transpose of PrKq, which must be diagonalised to obtain the correct 



eigenvectors, since Eq. (3.17) is an eigenequation from the left. For fixed values of 



Uq and we test a range of A values with separation 0.01 until we find a pair such 
that for one, the eigenvalue of PrKq is Isss than 1 and for the other, it is greater 
than 1. Previous analytical results at uq = indicate at which energies to look for 



the ground state 158 . A crude approximation for the excitonic energy is found 



by linearly interpolating between these two points. The process is repeated closer 
to this approximate value using a smaller separation of 0.002 to produce a set of 
coordinates {A, PrKq (A) eigenvalue}. We numerically find a function that interpo- 
lates these points and obtain a final solution for the excitonic energy by setting this 
equal to 1 and solving for A. The coefficients ^AfAf' may then be found by using 
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Figure 3.4: Excitonic energy as a function of magnetic flux $ at various values 
of electric field uq for interaction strength vq = — 2/7r^. Only every 4**^ data point 
is shown for clarity in each curve. Inset top: The energies of the first 3 states at 
^ = 0.45 as a function of uq. Inset bottom: The first 3 states at = 0.5 as a 
function of uq. Results displayed in insets were provided by Vivaldo L. Campo, Jr.. 



the obtained values of A and Gnn' in Eq- (3.12). 



3.3 Numerical results for the energy of the exciton 



Fig. 3.4 shows the resulting excitonic ground state energies plotted as a function 



of magnetic flux for different electric field strengths. For small enough electric 
fields, the excitonic energy oscillates as a function of the magnetic field as seen 
previously for uq = 157 . This is due to the electron and hole, which were created 



simultaneously at the same position, having a finite probability to tunnel in opposite 
directions and meet each other on the opposite side of the ring. The dependence 
of the oscillation amplitude, A(<1> = 0.5) — A(<I> = 0.0), upon electric field strength 
is shown in Fig. |3.5| for different values of the interaction strength. In all cases, 
there is initially a small increase in amplitude, as speculated in [181 and then it 
decreases strongly as a function of electric field strength. This change happens when 
the polarisation energy, (no|cos(0)|) = 2uo/7r, becomes comparable to the binding 
energy of an exciton. Using the binding energy for an exciton on a straight line, 
7r^fQ/2 158 , yields an estimate for the critical uq value, Uc = tt^Vq/A. As seen in 



the inset of Fig. |3.5[ a better fit to computed Uc values is given by the function 
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Figure 3.5: Amplitude of the Aharonov-Bohm oscillations as a function of external 
field uq at various values of the interaction strength vq. Only every 6*^^ data point 
is shown in each curve for clarity. Inset: critical external field function of 

vq. The dashed curve, Uc = 7r|t;o|(7r|t;o|/4 + 2/5), fits the numerical points quite well. 
The vertical dashed lines in the main plot indicate udvo). 



Uc = tt\vo\{tt\vo\/4: + 2/5). More precisely, the change of behaviour corresponds to 
a level crossing between the ground and first excited state, which only occurs at 



^ = 0.5, as shown in the insets of Fig. 3.4 The level crossing is associated with 
an exact symmetry at $ = 0.5 upon exchanging h ^ ~(^eh- It is also evident 



in Fig. 3.6, where the excitonic energy is plotted as a function of the electric field 
for various values of magnetic flux. The quadratic behaviour for small enough 
magnetic fluxes and the negative curvature of A are as expected from second order 
perturbation theory in uq. For magnetic flux values close to <I> = 0.5, the excitonic 
energy remains almost constant as a function of electric field until the critical uq 
value, when the level crossing occurs. 



3.4 Vanishing oscillator strength 



The OS is the strength of the transition from the excitonic ground state, ^'(^e, (^h), 
to the vacuum via photoemission. It can be shown to be directly proportional to 
the quantity 



^-271 

d9^ {6, 9) 

JO 



(3.18) 
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Figure 3.6: Excitonic energy as a function of uq for interaction strength vq = —2/tt'^ 
and various values of Only every 6*^^ data point is shown. The dotted lines show 
the results of second order perturbation theory performed on the ground state. The 
vertical dashed lines indicate Ur. 



where the wave function, ^{9e,9h), is correctly normalised 183| . Figure 3.7 shows 
F plotted as a function of magnetic flux for different electric field strengths; the 
AB oscillations are clearly visible. Interestingly, the oscillations become inverted for 
large enough values of uq. Most importantly, for these uq values, the OS vanishes 
at $ = 0.5. This means physically that an exciton created at some <I>, uq with 
F is prevented from recombining once uq are tuned to values for which 
F = 0. In addition, no excitons can be created when F = 0. Suppression of 
ground state exciton emission has already been seen for polarised excitons in the 



ring geometry for particular magnetic field strengths 1169 . However, this is different 



to the phenomenon we observe, since in Ref. [169], electron- hole recombination is 
prevented by confining the particles to different nanostructures. 

The transition from optically active to dark coincides with the change in 



slope of the amplitude of the excitonic AB oscillation as in Fig. 3.5 and also the 
reinstatement of the parabolic behavior in Fig. |3.6[ In order to explain the mech- 
anism for this behaviour, the probability density function for the exciton is shown 



in Fig. 3.8 at different values of -uq and vq. In Fig. |3.8[ a) the magnetic flux and 
electric field have been set to zero. The dark grey diagonal region indicates a peak 
in probability density where the electron and hole coordinates are close together, so 
the exciton is very small. In Fig. |3.8[b) the flux has been increased from its zero 
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Figure 3.7: Oscillator strength as a function of <I> for electron- hole interaction 
strength vq = — 2/7r^ and various values of the electric field uq. The horizontal 
dashed line indicates -F = 0. Only every 6**^ data point is shown for clarity in each 
curve. 



value, whilst keeping the interaction strength and electric field constant. There is a 
slight narrowing of the peaked region and the contours are closer together, showing 
that the exciton has become smaller. In Fig. |3.8[ c) the magnetic flux and interaction 
strength are the same as in Fig. |3.8K a), but the electric field has been increased. We 
see the dark grey region has begun to split in two, suggesting the exciton is begin- 
ning to separate. In Fig. |3.8[ d) the magnetic flux has been increased from its zero 
value in Fig. |3.8[ c). The result is two dark grey peaks indicating that the exciton 
has been broken |144] into an electron and hole, which are located on opposite sides 
of the ring. 

Figure |3.7| suggests a method for controlling the excitonic lifetime. The 
procedure is (i) begin with <I> = and uq = (ii) increase the flux adiabatically 
to <I> = 0.5 so an exciton can easily be created from the vacuum (iii) decrease $ 
(iv) increase uq until it appreciably exceeds the critical value (v) increase flux to 
^ = 0.5. The oscillator strength is now zero, so the exciton is trapped and unable 
to decay until the external fields are further tuned. This has important implications 
for excitonic data storage. 
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Figure 3.8: Excitonic probability density \^{9e,9h)\'^ at vq = —2/tt'^ and (a) uq = 0, 
$ = 0, (b) no = 0, = 0.4999, (c) uq = 0.5, $ = 0, and (d) uq = 0.5, $ = 0.4999. 
The solid and dashed lines are contour lines separated in height by l/lOvr^. 



3.5 Experimental work 



After publication of the work presented here, the XABE was observed for the first 
time in type-I structures in which both the electron and hole move in a ring ge- 
ometry [176 177 . In both of these experiments, the effect was enhanced using an 
electric field, raising it to a measureable level. Experimentally, AB oscillations in 
the excitonic energy appear on top of a quadratic dependence on the magnetic field, 



termed the diamagnetic shift 179 ,184 . One might also expect to see the linear 



Zeeman splitting, although this is only observed when single rings are studied, not 
for ensembles of rings [177 . 



In Ref. 176 , an array of self assembled InAs/GaAs rings were studied us- 
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Figure 3.9: (Taken from Ref. |177|). Density plot for the photoluminescence inten- 
sity of light emitted by an exciton in a InGaAs/GaAs quantum ring in the presence 
of external magnetic and electric fields, as a function of the field strengths. 



ing magneto-photoluminescence techniques. The emission spectra display a weak 
shoulder at lower energies indicating the presence of two emission bands. Clear 
oscillations were seen in the integrated intensity as a function of magnetic field 
for each band, although they have different periods and their maxima and minima 
occur in the opposite order. The different periods are simply due to different ring 
sizes. The inversion of the oscillation pattern is caused by an inbuilt electric field for 
the larger quantum rings, due to piezoelectricity from uniaxial strain fields formed 
in the growth process. It is known that such strains are particularly strong for 



InAs/GaAs rings 185 . It was also found using a theoretical model that inversion of 
the maxima/minima pattern in oscillations of the oscillator strength can be caused 
by increasing the temperature, so that excitons are excited out of the bright ground 
state into dark excited states. However, this was not the correct explanation for the 
observed behaviour. 



In Ref. [177], InGaAs quantum dots are grown using molecular beam epitaxy 
and capped with GaAs. Ring structures are then etched from these and the magneto- 
photoluminescence spectra studied for a single quantum ring. A vertical electrical 
field is applied, which can be used to alter the electron-hole separation and thus 
the period of AB oscillations. As suggested in our study, it is found that by tuning 
the electric and magnetic fields, the photoluminescence intensity is almost quenched. 



This is shown in Fig. 3.9, which was taken from Ref. 1177 . Subsequently, the XABE 



was studied theoretically in the presence of a perpendicular electric field 186 . The 
electric field was used to tune the ground state energy, so that weak AB oscillations 
were found. The amplitude and period of these oscillations can be tuned by varying 
the strength of the electric field. 
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Chapter 4 

Neutral collective excitations in 
pristine graphene 



4.1 The importance of interactions in graphene 

How to model e-e interactions is one of the greatest challenges of condensed matter 
physics. This is because at any given instant each electron moves in a potential, 
which depends upon the configuration of all the other electrons. Hence problems 
involving a few electrons are difficult to solve even using numerical approaches and 
those treating a macroscopic system are intractable via direct consideration of the 
forces at a microscopic level. To make progress, it is necessary to "zoom out" from 



these microscopic details. In the mean field approach 1120 , for example, the inter- 



action of an electron with other electrons in the system is modelled, by considering 
the remaining electrons as a smooth charge density representative of the average 
electronic configuration. In fact in many systems, the addition of e-e interactions 
to the non-interacting problem only results in quantitative changes. Such systems 



are well described by Landau's Fermi liquid theory 187] in which the problem of 
interacting electrons is replaced by one of non-interacting quasiparticles, which are 
elementary excitations of the system. This results in the renormalisation of their 
non-interacting properties. Systems in which e-e interactions induce a qualitative 
departure from the non-interacting behaviour are termed strongly correlated and 
are fascinating to study. 



The role of e-e interactions in graphene [188 is not fully understood. Some 



calculations suggest it is a 2D Fermi liquid 189 with a Fermi velocity renormalised 



by interactions, whereas others indicate strong correlations 190). An important 



parameter is the ratio of the characteristic energy of Coulomb interactions, Ei^t, to 
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^2 ^2 .2 

the kinetic energy, i^kin- For the 2DEG, E'kin = ~ 2mF' ^'^^ Fermi energy, 

2 

where i is the average electron separation and E'int = f^, where e is the dielectric 
constant which depends on the surrounding environment. Hence the ratio is 



-E'kin / 2DEG ' 

where oq is the Bohr radius. Thus we see that for the 2DEG, e-e interactions 
become more important relative to the kinetic energy as the electron separation is 
increased. Indeed below a critical electron density, the 2DEG will undergo Wigner 
crystallisation [62], forming a triangular lattice in order to minimise the potential 
energy of interactions. For the case of graphene, E'km ~ hvphp, yielding 

^ ~ — , (4.2) 
E'kin ehvF ' 

2 

which is independent of i. In fact = is called the graphene fine structure 
constant. Depending on the dielectric constant, ~ 10^^, which is large relative 
to the fine structure constant in QED, a = 1/137. The fine structure constant 
also describes the coupling between matter and light and a large value for graphene 
compared to non-relativistic materials means that it has a high opacity, considering 



it is only one layer thick, which has been confirmed experimentally 191 



Another important consideration is to what extent the e-e interactions are 



screened in graphene 192 . The Thomas Fermi wavelength, Atf, gives the length 
beyond which the potential from a given electron is screened out. In 2D, Atf — 
so for graphene Atf ~ '^^rT^ where rig is the carrier density. For neutral 
graphene Ue 0, as the Fermi energy is approached, since the conduction and 
valence bands touch at this point. Hence we conclude that for graphene with little 
or no doping, the e-e interaction potential can be treated as the long-range Coulomb 

2 

interaction, U {r) = j^, where the only source of screening is from e. 

In the presence of a magnetic field, e-e interactions play an important role. 



due to the single particle LL energies (calculated in Section 1.4.4) being highly 
degenerate. For the 2DEG 

Eint\^*^ . /e^me^c 1 ^^^^ 



-Ekin/2DEG ^ Be^n? -/B 



Hence Ei^t/E]^^^ depends on the magnetic field strength for the 2DEG and interac- 
tions can be made less important relative to the kinetic energy, by going to higher 
field strengths. In contrast, for graphene, -E'int/-E'kin = ck^, as in the absence of a 
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magnetic field. The average electronic separation cancels out as before, although for 
B it scales as £b- The size of ii^int/^kin in graphene can be considered medium. 
It is small enough to minimise spontaneous inter-LL excitations or mixing of LLs 
in the absence of an external energy source, but large enough not to be ignored 
altogether. In the absence of LL mixing, in a partially filled LL, the kinetic energy 
is effectively quenched and the interaction energy is the only remaining energy scale. 

In high magnetic fields, the LLs may be split into four sublevels due to Zee- 
man splitting (hujs) and an additional valley pseudospin splitting (tkOv)- This LL 
splitting is not yet well understood; it is even unclear under what circumstances 
spin or valley splitting is larger. Lifting of the four-fold degeneracy has been seen 
in quantum Hall experiments in high magnetic fields {B > 20 T) for the n = 
LL, with evidence for the n = +1 LLs observed, but somewhat weaker j94, 193 



2 

Large splittings on the order of the characteristic Coulomb energy, have been 
predicted |194| and observed |94|. Infrared spectroscopy experiments have found 
non-monotonic behaviour of the CRE as the filling factor of the zeroth LL is in- 



creased, which was interpreted as an indication of LL splitting 1951. Theoretical 



studies 196-201 state that it is a many body effect, due to some spontaneous 
symmetry breaking between the two sublattices. 

Another way in which the LL energies are affected by e-e interactions is via 
the exchange self energy correction, which is calculated below (Section 4.7). In 
graphene, energies for different LLs are renormalised differently, so that the energies 
for the optical excitation of an electron between LLs are altered and this occurs 
differently for different resonances; such effects have been seen experimentally. In 
contrast, in the 2DEG, the unique CRE remains unchanged by e-e interactions. 



This is known as Kohn's theorem 202 and is discussed in detail in this Chapter, in 
particular how it breaks down in graphene. Within the electric dipole approxima- 
tion, optically induced excitations carry no spin. The case of spin fiip transitions 
within the same LL, often referred to as spin waves, is dealt with for the 2DEG by 
Larmor's theorem [203 . It states that the energy of a spin wave excitation at wave 
vector k = is unchanged by interactions and remains equal to the bare Zeeman 
splitting, E (k = 0) = gjisB. Larmor's theorem also holds in graphene [204 . 



The central theme of the work on graphene presented here is collective exci- 
tations (CEs) of the monolayer system at zero temperature and in the presence of a 
strong perpendicular magnetic field. We focus on a system close to neutrality, where 
the conduction electrons behave as massless Dirac particles. Thus all calculations 
use ground states, with integer filling factors v = —2, —1,0, 1, 2, corresponding 
to the possible sublevel fillings of the n = LL. Considering only LLs close to neu- 
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trality, which have maximum separation for neighbouring levels, should minimise 
the effects of LL broadening present in real samples due to disorder and phonons; 
this work requires different LLs to be distinguishable. It should be noted however 
that the broadening of LLs in graphene is poorly understood at the moment, with 
qualitative disagreement between theory [205 , 206 and experiment [207 , 208] . This 
is hardly surprising, since an accurate theory requires knowledge of the type of dis- 
order causing the broadening, which is difficult to determine and may be different 
in different experiments. A neutral CE of the ground state occurs when an electron 
in one of the uppermost filled LLs within the Dirac sea is excited to an empty state 
in a higher lying LL leaving behind a hole; the bound state of an electron and hole 
is an exciton. The formation of charged CEs is also possible and will be discussed 
below. The term "collective" is used because of the LL degeneracy and also the fact 
that different excitations are mixed by the Coulomb interaction. 

This Chapter focuses on the clean system and presents some of the theory 
required for the calculations. This includes identifying the good quantum numbers, 
which allows classification of CEs according to the geometrical symmetries of the 
system. Classification of CEs according to their dynamical symmetries is discussed 
in Chapter [5j Optical properties of the CEs are of particular interest, so the selection 
rules are derived here for circularly polarised light. In addition, I contextualise the 
work in terms of pre-existing research. 

The majority of work on CEs of graphene in a strong perpendicular mag- 
netic field follows a similar approach to that used to study the analogous problem 
for the 2DEG 209-211 . In the 2DEG, bound electron-hole (e-/i) complexes in a 
magnetic field have a well defined linear momentum of magnitude k. A classical 
way to visualise this is that at a certain e-h separation, the attractive and repul- 
sive forces exerted on them, due to their Coulomb interaction and the Lorentz force 
respectively, cancel. Hence the excitonic dispersion relation is well defined with 

— 2 



the e-e interactions causing a deviation from the bare CRE of order £"0 = ^2 tIb 
In Ref. [210| the dispersion relation for both filled and half-filled (spin polarised) 
LLs is calculated using a Hartree-Fock approach to first order in the parameter 
£^int/-E'kin- For small values of k, the dispersion relation looks like that of magneto- 
plasmons (charge density waves) or spin waves. Corresponding studies of graphene 
have found excitation spectra qualitatively similar to those for the 2DEG, but with 
stronger many body corrections [212 213 . Chapter [7| addresses the addition of 
impurities and the resulting localised neutral CEs, which appear as discrete states 
below the band of extended CEs. The identification of these localised states re- 
quires knowledge of the width of the band and so for completeness we review in 
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Section 4.8 the derivation of the excitonic dispersion relation for graphene. In addi- 



tion to spin and valley pseudospin preserving CEs, spin flips, spin waves 204 and 



skyrmions 194 214] (coherent charged topological CEs) have also been studied in 



the context of graphene in a perpendicular magnetic field. 

4.2 Neutral and charged collective excitations 



Recall from Section 1.4.4 that in the symmetric gauge, single electron states in 
graphene are labelled by the quantum numbers n, r, s, m. Let Cnrsm be the fermionic 
creation operator for an electron in graphene, so that ^'A/'m(r) = (r|c]y^j^|0), where 
J\f = {tits} is a collective index and |0) represents the vacuum. I shall use the 
hole representation, cj^j-m —>■ d)j\fm ^Mm ~^ dj\fm for all filled levels. The usual 
anticommutation relations for fermions are obeyed: 

{cMmi : ^ = ^m^m^ , {'^Mma > ^jv^mi, ^ = ^Af2Af^ ^m^m'^ ■ (4-4) 

Consider a neutral CE, corresponding to an electron being promoted from a 
state M2 to an empty state A/i in a higher lying LL. This excitation can be created 
from the ground state, {v}, by the creation operator 

Qm,N2 = S ^^1^2 (mi , ) c]^^^^ S^^^^ . (4.5) 

m,i,m2=0 

Suppose now that the filling factor \s u = + e with e « 1 and /i = 0, +1, +2. 
Then an integer number of sublevels of the n = LL are filled and there are a few 
electrons in the next sublevel (see Fig. |4.1[ a)). In this case it is possible for an 
neutral CE to bind to an additional electron forming a negatively charged CE or 
trion, . This excitation can be created by acting on by the creation operator 

00 

^MMATs = S ^MA/-2A/-3 i^i , m2 , mg) ^3^3 S^^^^ . (4.6) 

mi, m2, 013=0 

Similarly if the filling factor is = /i — e, so that the uppermost filled sublevel is 
almost completely filled, a positively charged collective excitation, X"'", may form 



(see Fig. 4.1 ^b)). The creation operator for the X'^ can be formed in an analogous 
way to that for the X^ . As we shall see in Chapter [6} the charged CEs appear 
as discrete states below the excitonic continuum. The existence of trions in semi- 



conductors was first proposed in 1958 215| and observed experimentally during the 
1990s [216] . 
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Figure 4.1: Schematic diagram of the formation of (a) X trions at fihing factor 
V = 2 + £ and (b) X'^ trions at fihing factor v = 2 — e. 



4.3 Hamiltonian of collective excitations 

The Hamihonian, H = Hq + U^^ + Uee + Uhh-, for coUective excitations is 

M ,m M ,m 

Ml ,Af2 '"^2 

A/';^,A/';^'"2.'"3 

It is important to note that the range over which given M quantum numbers are 
summed, depends on whether they index an electron or hole operator. The in are 
the single particle energies, which have been renormalised by self energy corrections 



(see Section 4.7 below). I assume that the previously mentioned spin and valley 
splitting is negligible compared to the LL energy: hu^ < hug « ei and hence set 
hujy = fujjs = in all calculations. The only effect of such splittings is to fix the spin 
and pseudospin quantum numbers for the different ground states. 

The last three terms give the pairwise particle interactions. We assume 
that the direct interaction with the Dirac sea is cancelled by a neutral positive 
background charge from the carbon atoms. The lA are the matrix elements of the 
Coulomb interaction, U {\vi — r2|) = e^/e|ri — r2|, in the basis of single electron 
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wave functions (see Eq. ( 1.29| )) defined by 



J J dr^dr^^-;^, (ri)®^')^, (r2)t/ (|ri - ral) ^-TVi (ri)®^iV2(r2), 

(4.8) 

M' m! M' m' 

where = {nmr] is a collective index. Note that Z^a/^ \r^^^oC(5„, „/ 5„„o' (5^,-r"5T-ox' ; 
the long-range Coulomb potential cannot provide enough momentum to scatter be- 
tween the valleys [217|. Eq. (4.8) can be written in terms of the 2DEG matrix 
elements, 



Unll^lni^i = J Jt^r?*2V'4m',(ri)V'4m'2(r2)f^ (^1 - r2 1 ) Vnama (r2) V'nimi (ri) , 



according to 



(4.9) 



^NimiN2m2 



^sxs'^ 5s2s'2 ^txt[ ^T2T^ CLni ^n'j 



Tj\K\"^'l \"-2\"T-2 
\ni\mi \n2\m2 



"■2 "2 |ni|mi |n2| — lm,2 



J-*? rrhi|-l»7i'i |n'2|-lm'2 

-hO„iO„20„'_^0„/^Ly, |n2|-lm2 



(4.10) 



The first of the two-body terms in Eq. (4.7) gives the e-h interaction, tlehi where 



-^A/-{m'iA/->'2 
''^Nim-iN2m2 



, ,N'[m'^M2'm2 j,M2m2N'[m\ 



(4.11) 



The two parts correspond to the direct and exchange interactions; the latter com- 
prises the random phase approximation. The last two terms, Uee and Uhh^ give the 
e-e and h-h interactions respectively. These terms are only needed when considering 



charged CEs. Techniques for calculating the matrix elements in Eq. (4.9) are given 
in Ref. |218| and Appendix |B] 



4.4 Geometrical symmetries of the interacting Hamil- 
tonian 



We wish to diagonalise the Hamiltonian, H, in order to obtain the eigenenergies 
and amplitudes of CEs. This complex problem can be made easier by determining 
the maximal set of mutually commuting operators containing H, since these op- 



erators will have simultaneous eigenstates. According to Noether's theorem 219 



corresponding to each symmetry of a physical system is a conserved quantity and 
by Heisenberg's equation of motion, the operators of conserved quantities commute 
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with the Hamiltonian. Thus to find a set of mutually commuting operators, it 
makes sense to consider the symmetries of the system. We examine at this point 
only the geometrical symmetries; the dynamical SU(4) symmetry arising from spin 
and pseudospin degrees of freedom is covered in Chapter [5] 

Let us consider a clean monolayer of graphene in the symmetric gauge. 
Then there is rotational symmetry and the resulting conserved quantity is the total 
generalised angular momentum projection, J^. This corresponds to the operator 
'^z = where {jz)i is the generalised angular momentum projection for the 



i**^ electron defined in Eq. (1.30). One can show explicitly [J^,-^] = 0, confirming 
that Jz is a well defined quantity for CEs. For a neutral CE consisting of an elec- 
tron in an A/i state and a hole in an N2 state, Jz = \ni \ — mi — \n2 \ + m2 and the 



sublattice parts cancel. Thus one of the summations in Eq. (4.5) can be removed 



and the operator for neutral CEs can be re-written more simply as 

00 

mi=0 

where m2(mi) = Jz — \ni\ + mi + \n2\. For the state, Jz = \ni\ — mi — \n2\ + 



^2 + \n^\ —ms — h, so that one of the summations in Eq. (4.6) may also be removed. 



The quantum number Jz has a direct geometrical meaning for neutral CEs 
determining the average positions of the electron and hole relative to the origin. 



Using Eq. (1.31), one may show 



(MiM2JzK - 4WMz) = {2Jz + 4(|n2| - |ni|) + 5„,,o - 5„,,o)4, (4-13) 

where \M1N2Jz) = Q^MiN'2Jz\^'^' Therefore, the sign and magnitude of Jz (together 
with the LL indices ni and 712), determine whether the electron or hole is closer 
to the origin. This has physical relevance for the systems studied in Chapter [7| 
which have a single impurity located at the origin. For example, for the excitations 
involving the initial and final states in the n = LLs (ni = and n2 = 0), we 
get (rl - tDo^o = 2Jzi ^, so that for Jz > 0, for example, the excitation is a 
low-energy (high-energy) state on the donor (acceptor) impurity. For transitions 
between the zero and the first LL's in Graphene {ni = 1 and n2 = 0) we have 
(r^ — rg)o^i = {2Jz — 3)^^, while for transitions between the LL's n2 = —1 and 
ni = 0, (rl-rlU^o = {2Jz + S)£l. 

Due to the homogeneity of space, the operator of magnetic translations is 
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it k 

k = k=] k = 2 

Figure 4.2: Schematic diagram illustrating how the operator simultaneously 
raises the value of k and lowers the value of Jz- 



conserved. For a single particle it is 



-Ti X B. 



(4.14) 



The total operator summed over all electrons is then k = kj. Noting that (in the 
symmetric gauge only) k(B) = n(— B) simplifies calculations. Conservation of k is 



easily seen from Eq. (4.14), since 



dt 



dt 



Qi dvj 
c dt 



X B = 



(4.15) 



by the Lorentz force law; in Eq. (4.15) and below t is used to represent time. This 
together with the fact that k does not depend explicitly on time, yields [H, k] = 
by the Heisenberg equation of motion. However, [k^, ky] for charged complexes 
and [J^jk] 0. Specifically 



ifikf 



(4.16) 



Instead one can show [ j^, k^] = and hence H, Jz and k^ are mutually commuting. 
The eigenvalues of P are {2k + 1)(^/£b)^ with A; = 0, 1, . . . [22oj. k is called the 
oscillator quantum number and indeed /c = m for a single particle. Analogously to 
the single particle case, k determines the distance from the origin of the guiding 
centre of a few particle complex. For a uniform system this leads to the Landau 
degeneracy in k. Few particle complexes occur in families starting with a seed state, 
which has A; = and Jz = J^^', the latter cannot be guessed a priori. The off- 
1,2,... are generated by successive action of the r aisin g ladder operator 
iky\, since [k^,A;_] 



spring k 



2{h/lBYk^ (see Fig. 



creasing eigenvalues of Jz equal to ji'^-' — k because of the relation \Jz,k 



4.2|). They have de- 
= —hk^. 



which is clear from Eq. (4.16). In addition to having the same energies as each 



other, states in a given family also exhibit the same optical properties. 

In order to use the k quantum number to remove a degree of freedom, a Bo- 
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goliubov transformation of the basis must be made. This is a challenging approach, 
which unfortunately we did not have time to investigate. Translational symmetry 
is broken in any numerical calculation, since they can only treat finite system sizes; 
the system size is controlled by the value at which summations over the m quan- 
tum number are cut off. Thus even in our calculations for charged CEs in pristine 
graphene, k is not strictly speaking an exact quantum number. However, for large 
enough system sizes the eigenstates have a large amplitude for a single k state and 
in Chapter |6] we are able to clearly see the degeneracies described above in the dis- 
crete states which appear below the continuum. All the original work presented here 
on neutral CEs includes the effects of a single impurity thus breaking translational 
symmetry, so k is not a good quantum number. If we were to remove the impurity, 
k would in principle be defined, but then the states would form a continuum hiding 



the families. In Section 4.8 below, previous work on neutral CEs of clean graphene is 
examined, but in this case the Landau gauge is used resulting in different quantum 
numbers. 



4.5 Optical selection rules 

We begin by determining the perturbation to the single particle Hamiltonian given 



in Eq. (1.27) by an incoming circularly polarised beam of light. We work in the 
electric dipole approximation, ignoring the magnetic component of the light wave. 
The electric field for right (left) circularly polarised light, which I shall refer to as 
the (c ) polarisation, is the real part of the complex field, E = £e^^'^^{ex + ie^), 
where £^oc|E|, a; is the angular frequency and e^,. and are unit vectors in the x and 
y directions respectively. Using Weyl's gauge in which the scalar potential vanishes, 
we obtain E = —d{5A)/dt, which gives 5 A = —i'E/io, where 6 A is the correction by 
light to the vector potential A due to the constant magnetic field. It follows that 
the correction to the single particle Hamiltonian due to incoming a- light is 

6H+ ■■ 



(4.17) 



where a+ = ax i: iay. Here we have assumed that any linear momentum transfered 
to the system by a photon is negligible, so the excited electron cannot change valley, 
which would require a large momentum boost. In addition, the interaction with 
light does not change the spin of an electron, at least within the electric dipole 
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approximation. Using (un)primed symbols for the (initial) final state of the electron, 
these two conservation rules are expressed as t = t' , s = s' . 



Using Eqs. (1.29) and (4.17) the single particle optical selection rules can be 



easily found by calculating the dipole matrix elements, since the oscillator strength 
is directly proportional to the modulus squared of these matrix elements. The 
selection rules for the orbital quantum numbers for an electron, with r = r' =fl~ and 
illuminated by light for example, can be found by considering the dipole matrix 
element 

(n il sm' \ 6H+ \ n il sm) oc {(\ n'\,m' \ Sn'(\n'\ - l,m' \) ^ (4.18) 

^ f I \n\,m} \ 
\ Sn I \n\ -l,m) J 
oc 5„'(|n'| — l,m' \ |n|,m) 

Smm'- (4.19) 

From this and similar calculations, the single particle optical selection rules are 
found to be 



m 



m 



n 



n 



+1 



(4.20) 



for the a- polarisations, irrespective of the valley index (so long as it is conserved). 



Eq. (4.20) is the starting point for calculating the optical selection rules for 



CEs. The correction to the many body Hamiltonian in Eq. (4.7) due to incoming 
a- light is 

6n+ = J] (Af'm' I 6H+ I AAm>4„,4^. (4.21) 

Af,m 
Af',m' 

It can easily be shown that (A/iA/'2 J2|(5?^+|i^)oc5j^^+i, so that the selection rules for 
a neutral CE are = +1 and Sz = Tz = 0, which are the spin and pseudospin 



projections of the CE. Eq. (4.21) will be used to calculate the oscillator strengths 
for CEs in later Chapters. 



4.6 Mixing of excitations 



Excitations created by Eq. (4.12) with the same value of Jz, but different values 



of ni,n2 are mixed by the Ueh interaction term in the Hamiltonian H. The true 
neutral CE is thus given by a superposition of the different excitations, |A/iA/'2 Jz). 
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This is best illustrated by an example. Consider illuminating a graphene sample 
with filling factor v = — 1, . . . , 2 with fj"*" light of frequency Coc (the tilde denotes 
renormalisation of the bare frequency due to e-e interactions). Then electrons will 
be excited from the n = LL to the n = 1 LL (0 1). These excitations will be 
mixed to others with the same Jz, for example —1^0 excitations (for = — 1, 0, 1). 
The latter are in resonance with the 0^1 excitations and thus significantly mixed. 
However, excitations that are not in resonance, —l—>2 for example, are mixed with 
amplitudes which depend on ~ 10^^. We neglect such virtual transitions and thus 
only ever consider the resonant 0^1 and —1^0 excitations; for the 2DEG this is 
called the high magnetic field approximation. Even for the case ~ 1, the mixing 



in the 2DEG between non-resonant excitations has been seen to be small 1221 



In graphene there are 16 resonant excitations, shown in Fig. 4.3 for the 



particular case u = —1. The Kronecker's deltas in Eq. (4.10) together with Eq. 



(4.11) restrict the possible spin and pseudospin states of mixed transitions. Two 
excitations, \M1M2Jz} and |A/\'A/'2 J^), are mixed by (i) the direct interaction if si = 
s\, Ti = t[, S2 = s'2 and T2 = and (ii) the exchange interaction if s\ = S2, 
Ti = T2, s'l = s'2 and t[ = The energies and eigenstates (the numerical values 
of ^A/'iA/'2Jz) foi^ CEs are computed numerically by cutting off the summation over 



the m quantum numbers in Eq. (4.7) at a finite value and diagonalising the matrix 



composed of the matrix elements of the Hamiltonian 

(M[rn[j\f2rn'2\H\MirniJ\f2rn2} = '5m2,J.-|ni|+mi + |n2|'5m'2,J,-|ni|'+m'i + |n'2|(4.22) 



'■^'^A/'imiA/'2ni2 )' 

This truncation of the basis corresponds physically to choosing a finite system size. 



4.7 Self energy corrections to single particle energies 
4.7.1 Kohn's theorem 

To see how Kohn's theorem breaks down in graphene, let us first examine the reasons 



that it holds for the 2DEG. Kohn's theorem 202 states that for the interacting 



electron gas in the presence of a magnetic field, the CRE is independent of the 
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Figure 4.3: Diagram of lowest energy inter-LL transitions for v = —1. Shading 
indicates the Dirac sea, open circles (o) holes and closed circles (•) electrons. Black 
transitions are optically dark, whereas coloured transitions are optically active. Red 
(blue) indicates transitions which are active in the (<7 ) polarisation provided 
they have = 1 (Jz = — !)• Excitations in the dashed box are mixed by the 
exchange Coulomb interaction, whereas excitations connected by a dashed line are 
mixed by the direct interaction. 




interaction strength. The Hamiltonian can be written as 



2mf 



+ U, 



(4.23) 



where is the number of electrons and U = '^iju{ri — r^) is the interaction 
potential. If the momentum of the whole system is defined as II = "^itli, then it 
can be shown that 



dn 

~dt 



rrieC 



n X B, 



(4.24) 



which is the many electron version of the Lorentz force law. The derivation of Eq. 
(4.24) uses nj = 0, which holds 3ecause II acts only on the centre of mass 
coordinate, which U is independent of. Defining II-i- = II^: + illy and using Eq. 



(4.24) gives 



(4.25) 



Therefore if ^'o is an eigenstate of H with energy E^,^i = n+^'o satisfies 



EG^ 



^-1 



(4.26) 
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and so is also an eigenstate with energy ^1 = E'^ + hLof'^. The perturbation 5H^^ to 
the Hamiltonian H due to incoming circularly polarised light in the a- polarisation 
can be derived within the electric dipole approximation in a similar way to that 



for graphene given by Eq. (4.17). Assuming light of a long enough wavelength 



(wavelength in infra red regime or higher) 



ie£ 



nipUjc 



-n+e 



—iwt 



(4.27) 



Eq. (4.27) indicates that an electron in state \I'o, which is excited by a photon can 

jc- This 



only be promoted to ^i, so that there is a sharp absorption peak dl uj = Uc 
shows that the CRE is independent of the interaction strength, proving Kohn's 
theorem for the 2DEG. 

To summarise, the CRE for the 2DEG is unaffected by interactions due to 



Eqs. (4.26) and (4.27). For the case of graphene, we see from Eq. (4.17), that 



the Hamiltonian, describing the coupling between light and electrons, is di- 

rectly proportional to a+ or cj- , depending on the valley and polarisation direction. 
However, [lfD,o"+] is not directly proportional to cj+, where Hd is the graphene 



Hamiltonian defined in Eq. (1.24), so there is no graphene analogy to Eq. (4.25) 



Thus we conclude that optical excitation energies in graphene can be renormalised 
by e-e interactions and proceed to calculate these self energy corrections. This topic 



has been addressed in several theoretical papers 212,213,222,223 . 



4.7.2 Renormalisation of single particle energies 

Consider the ground state of graphene in a magnetic field with an integer filling 
factor. An electron in a given LL, n, will interact via the exchange interaction with 
each other electron in the system with the same spin and pseudospin, irrespective 
of what LL they are in. The two-electron state remains unchanged by such an 
interaction. The overall effect is that the energy of this electron is lowered below 
the bare single particle energy, e-n- Let's derive this self energy correction using the 
symmetric gauge. We have seen that the exchange interaction energy for an electron 
in state Mimi and a hole in state ^2^,2 is given by the matrix element ^^^mfjv^^m2 ■ 
Analogously, the interaction energy due to exchange between two electrons in states, 
Mm and M'm! is given by ~^vW/^n^' '^lisre the sign change is due to switching 
between a positively charged hole and a negatively charged electron. Indeed, self 
energy corrections are easier to visualise by thinking in the electron representation 
as opposed to the usual e-h representation. Since 

Kmj^^ > 0' the exchange inter- 
action between two electrons is attractive. One may assign the energy — ^^jVm^VW 
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to each electron. From Eq. (Q, U^'^^^ocd^s'Srr' and W^^^r'n™ depends only 
on n, m, n' , m' , not on the specific value of r, s. Let E^e {n, n') denote the correction 
to the energy of an electron in the n}"^ LL due to exchange with electrons in the n'*^ 
LL, so that 

^SE '^O = ~ 2 2 ^nrsm nVW • (4-28) 
m' 

Clearly E^E{n,n') = EQ,E{n\n). As usual, the graphene matrix elements can be 
expressed in terms of 2DEG matrix elements according to Eq. (4.10). Details of 



how to calculate E'sE {n,n') are given in Appendix [B| where it is shown to be inde- 
pendent of m. The complete self energy correction for an electron in the n}^ LL is 
of course given by summing E^e {n,n') over all filled LLs n' . We shall see that the 
bare single particle energies are renormalised differently for different n, leading to a 
renormalisation of the optical excitation energies. 

Let us now focus on the particular example of z/ = —2, ... 2 and calculate the 
self energies for different inter-LL transitions. The possible transitions for i/ = —1 
are shown in Fig. |4.3[ they fall into two different classes with the members of each 
class having the same self energy correction. All self energy corrections for tran- 
sitions occuring at other partial fillings of the n = LL are identical to the self 
energies of one of these classes. To proceed we need the following results, which 
hold for n ^ 1: 

^se(0,0) = -^0, (4.29) 

EsEiO,-n) = -^^-i^, (4.30) 
Za/tt n! 



£se(1,0) = -^, (4.31) 

' ' 8Vvr T{n + 1) 2n - 1 ' ^ ' 



P . 1 X ^0 r(n+^) 8n + 4^-l 

EsE (-1, -n) = ^^7^1—' (4-33) 



where the Gamma function is as defined in Eq. (2.26). Li general, higher lying 
LLs are lowered less by the exchange compared to lower lying LLs, for example 
l-EsE (1)0) I < \Ea,E (0,0) |. This means that although the self energy correction to 
each bare LL energy is negative, the overall self energy correction to the optical 
excitation energy is positive. 

Let's first consider the transitions with no spin or pseudospin flips, enclosed 



by the dashed box in Fig. 4.3, beginning with the 0^1 transition. The total self 
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energy correction to the transition energy is 



tpIts 



5ei -5eo=Y, (^se (1, -n) - Ese (0, -n)) 



(4.34) 



n=0 



where ric is the value of n at which the sum is cut off. This should correspond to 
the region in the zero magnetic field case, where the dispersion relation deviates 



significantly from being linear (see Fig. 1.8). According to the derivation of Eq. 

1 



(1.23), this occurs at the critical wave vector, qc 



of ric can be determined from 



1^ 



Hence an appropriate value 



(4.35) 



which yields 



16300 



(4.36) 



where it is understood here that B is the magnetic field in Tesla. This approach is 



similar to that used in Ref. 222 to estimate n^, although there the authors used 



where VSa corresponds to the lattice parameter of one of the triangular 



graphene sublattices. The estimates for Uc derived here and in Ref. 1222 are both 



lower than that given in Ref. [212 . Using Eqs. ( |4.29[ )- p.32 ) one can show 



Its 
Ors 



^.Eo + ^Y. 4vn 



8V7r 



2n- 1 



(4.37) 



At large values of 



Its 
Ots 



Eq f 



00 

V 



V^r (n - I) 1 



n! 



n 



+ 7 + 



+ Inric 



(4.38) 



where 7 is Euler's constant. Hence -Eg^f diverges logarithmically with Uc and for 
very large values of ric is hardly altered by small changes in ric, which is a pleasing 



result. Reassuringly Eq. (4.38) agrees with the result in Ref. 212 , where it was 
calculated using the Landau gauge. For B = 20T, ric * 810 and £"0^;! k 1.56Eo. 

One can show Eq'^^ = E^^^^, which is due to the particle- hole symmetry 
in graphene. Hence the energies of collective excitations comprised of the four 
excitations without spin or pseudospin flips are only shifted by a constant due to 



self energy corrections. It is clear by studying Fig. 4.3 (see page 73) that all the 
transitions which are not mixed have the same self energy correction, -Eq", as the 
four transitions which are mixed by the exchange interaction. Further calculations 
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show 



oitt ~ offt ~ of[t ~ -ifrt ~ -ifrt ~ -Iff! ~ Ors "I" ^ ) i^^.oyj 

so that all transitions, which are mixed to exactly one other transition by the direct 
interaction, have the same self energy as each other. 

Let's consider the self energy correction for another transition with a different 
bare energy for comparison. For simplicity I choose the filling factor v = —2, so the 
Fermi energy lies between the n = —1 and n = LLs. All sixteen possible —1^0 
transitions will then have the same self energy correction E^^^^. The —1^2 
transition has the lowest bare transition energy, hcodV^ + 1)) following the —1^0 
transition. The self energy correction for the —l—>2 transition is 

E^^i:^\^ = 6e,-6e^i (4.40) 
= Y,(^SE{2,-n)-EsE{-l,-n)) 



32-.hr ' 2 



^^r(n + i) 



n=2 



r(n + 1) 



32(V2 + l)nf - 16n - (36V2 + 48)ni + 3 
(2n - 1) (2n - 3) 



For B = 20r, Tic ~ 810 and E'^^^^^^ k 2.69^0- Note that the self energy corrections 
do not scale in the same way as the single particle energies. 



4.8 Dispersion relation for magnetoexcitons 

Below I present the derivation of the magnetoexcitonic dispersion relation in graphene 



following in part the approach outlined in 2121. It is similar to that used for the 



2DEG (210]. The non-interacting e-h wave function, which depends on momentum, 



is found by manipulating the Hamiltonian through coordinate transformations and 
a gauge transformation until it can be formulated in terms of Bose ladder operators. 
This is used to find the second quantised operator, which creates such an e-h pair 
out of the Dirac sea. The dispersion relation is found by diagonalising the Coulomb 
interaction matrix in the basis of these non-interacting two-particle wave functions. 
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4.8.1 Derivation of non-interacting wave function 

We begin by finding the non-interacting e-h wave function using tlie Landau gauge. 
Let's assume both the electron and hole are in the K valley. As previously remarked, 
the long-range Coulomb interaction does not scatter between the valleys. The valley 
index of either particle can be changed by interchanging the A and B sublattice 
components as usual. The two particle Hamiltonian is then 



Ha 



Hoe 012 + 120 Hi 



Oh 



(4.41) 
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where -ffoe(ri) and ifo/i(i"2) are the single particle electron and hole Hamiltonians 
respectively, as given by Eq. ( |1.24 ). After making the coordinate transformation 



ri + Y2 



y^y 



y^y 



ri - r2 



(4.42) 



and the gauge transformation U = exp ^— ^^p^j , we obtain Hq = U^HqU. Now Ho 

satisfies [Hq, —iKVn] = 0, so that the canonical centre of mass momentum, P, is a 
good quantum number. Since we wish to find simultaneous eigenstates of Hq and 
— iWR,, we replace the operator — i Wr, in Hq by the quantity P. The R dependence 
of eigenstates of Hq is then included by the factor exp(zP.R/^). The dependence of 
Hq on P can be removed by making another coordinate transformation 



il{e, X P)/h. 



(4.43) 



This gives 



where 



Hq 



hujc 



( 

h 

at 
V 



I I Z 



2iB 



6t a \ 

0a 

6t 

at 6 



6t 



i ( ~z* 



(4.44) 



(4.45) 
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and z = Vx + if^y Notice that the operators a^, are the same up to a phase as 
those defined in Eqs. (1.6) and (1.8): 



-at, 6t 



ib^ . Hence 



V;.^(f) = <f|^^(at)'^(6t)"|00> 
\ n\m,\ V / 



(-l)"+'"i"V'nm (f) 



(4.46) 



We can now construct the normalised e-h eigenstates of Hq as 



ifnm (r, R, P) 



a„amexp 



''B 



exp(iP.R/?i) 



/ Sn4'\n\-l,\m\ (?) \ 

•SnSm'4'\n\-l,\m\-l {^) 

^|n|,|m| (r) 

V •Sm^\n\,\m\-1 (?) / 

(4.47) 



These have corresponding eigenenergies (5n-\/|"'l + 'Sm\'\'m'\)hujc, which is the sum 
of the electron and hole single particle energies, as one would expect. Note that m 
here represents the LL index of the hole, not an oscillator quantum number as in the 
single particle case. Thus it is more appropriate to use the notation ni,n2 instead 
of n, m. In addition, the order of the components of the coordinate part of the wave 
function depend on the pseudospins of the electron and hole and of course both 
particles also have spin. Hence we switch to the notation v^a/'iM (riR)P)i which 
is connected to the old notation via ^^niitsmafrsa (r,R-, P) = V'nma (r, R, P) XsiXs2- 



The analogous e-h wave function for the 2DEG is very similar to Eq. (4.47): 



(r,R,P) 



exp 



exp(zP.R/^)'0„^„2 (f) 



(4.48) 



4.8.2 Determining the dispersion relation 



The next step is to express the state (/?a/'iA/'2 (r, R, P) in second quantised form. In 
other words we would like to find the operatoij^ 



QmA/-2(P)= S ^MXi,A/-2X2(P)cViXi^^A/-2X2' 
Xl,X2 



(4.49) 



^Note that Jz is not a good quantum number for these CEs, since we are using the Landau 
gauge. 
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which satisfies 



= v?MM(r,R,P). (4.50) 



The creation operators for electron and hole, c^n^T-^siXi ^^"^ ^ri2T2S2X2 respectively, 
are the Landau gauge analogy of those defined for the symmetric gauge. The 
^A/'i^i,A/'2X2(P) ^'^^ found by using the orthonormality of the c/)^"^^ (r) and 4>^^'^ (r), 
giving 



1 ^ ( .XP, 
2_, exp I I 



X 



) ^AfiX-Pyel/2h"'Af2X+Pyel/2h- 



(4.51) 



This is useful, since we already have the second quantised form of the interacting 



Hamiltonian. It is given in Eq. (4.7) for the symmetric gauge, but the correspond- 
ing operator for the Landau gauge is simply obtained by changing all m ^ X 



in Eqs. (4.7)-(4.11 ), all ^ $ in Eq. (4.8) and ah V 



in Eq. (4.9). Let 



|A/iA/'2P) = Q^^_^j^^{P)\h'}. As discussed in Section 4.6, several states |A/iA/'2P) are 
mixed together by the Coulomb interaction to form the CE. The dispersion relation 
is then given by diagonalising the matrix with elements 



(MiM2P'\H\MiM2P} 



I V— \ i 



Ux-x')p^ 



^0 

X (Cni - e 



^AfiAf[^M2M!,^XX' 



n2 



+ W 



,Af[x'-Pye%/2h AfliX'+Pyi%/2h 



MiX-Py£%/2h M2X+Pyi%/2h 



3pp/ 



(4.52) 



^NxM[0N2M'2\^n^- ^^2)^1^ M^M2 MxM2 



The first term in Eq. (4.52) represents the exchange renormalised single particle 



energies and is independent of P. The latter two terms are not and give the direct 
and exchange contributions respectively. They are defined as 



f.dM[M'2 
" MxM2 



L V Mx-x')p^jj^ix'-Pyel/2hM2X+Pyel/2h , . 

V. 2-1^ ^MiX-Pye%/2hAf^x'+Pye%/2h ^^■^'^^ 



XX' 



and 



" M1M2 



XX' 



i{x-x')p^.,J^2X+Pyf.l/2hM[x'-Pyil/2n 
^MiX-Pye%/2nAf!^x'+Pyel/2h- 



(4.54) 
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The graphene matrix elements, U, can be expressed as a combination of 2DEG 



matrix elements, f7, as we saw for the symmetric gauge in Eq. (4.10). Hence Eqs. 
(4.53) and (4.54) can be written in terms of Ij'^nin^ a-nd lJ^nxn2 ■, which are defined 



by changing M ^ n and lA ^ U vcl Eqs. (4.53) and (4.54) respectively. One can 
show 

C/'ni? = (-1)"^+"'^+^ Jd^rV;:,^, (r) \J (|r + 4(e, x P)//i|) V5„,n. (r) , (4.55) 
which is derived using 



^1%. (r,R,P) 



h 



niX-Pye%/2h 



and Eq. (4.48). The exchange interaction can also be simplified as 



U\7n7 = {-ir+<U[P/K)i,l^^^ (e, X P) (e, x P) . (4.57) 



From Eq. (4.55) one can see that taking the matrix element of the Coulomb interac- 



tion with respect to the two-body wave function, (fj^f-^j^f^ (r, R, P), yields the direct 
interaction only. To obtain both the exchange and direct terms, the many body 
approach described here is necessary. 



4.8.3 Worked examples 
Intra-LL excitations 

We begin with a simple example of the dispersion relation for the case when both the 
electron and hole are in the n = LL. The possible intra-LL neutral CEs are shown 
in Fig. 4.4 (inset) for v = —\] they are not mixed. Since we neglect the energy 
separation between sublevels, there is no single particle contribution. Also each 
transition involves a spin or pseudospin flip, so there is no exchange contribution. 
The direct matrix element is given by 



^ 0TiSi0r2S2 



fjd 00 
^ 00 



-Eoe 4 Jo , ^ 



(4.58) 



where V = PIb/Tt- and /„ ai^e modified Bessel functions of the first kind. Notice that 
the matrix element depends only on the modulus of P and not its direction, which 



is to be expected due to the symmetry of the system. Fig. 4.4 shows the dispersion 



relation. Since the self energy is independent of V, we set Eq^'^^'^ = 0. By taking 



the limits V ^ and P ^ oo in Eq. (4.58), the band width for the case when both 
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V 



Figure 4.4: Dispersion relation for transitions within the n = LL, given by the 
expression If^QQ (see text). Inset: Diagram of the possible intra-LL transitions for 
filling factor v = —1. 



electron and hole are in the n = LL is seen to be Eq. 



Inter-LL excitations 

Here I present the dispersion relation for inter-LL magnetoexcitons for the case 



when an integer number of sublevels of the n = LL is filled (see Fig. 4.3 page 73). 
As described in Section 4.6, we only take the 0^1 and —1 excitations into 



account. We require the following matrix elements: 



d = U 



d Insi 0r2S2 
Irisi 0t2S2 



d Orisi — 1t2S2 
Orisi — 1t2S2 



frd 10 I frd 00 
■-^ 10 + 00 



e 4 



(6 + lo 



(4.59) 



^' _ Insi 0t2S2 
" ~ Onsi -lr2S2 



d 11 
00 



(4.60) 
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(4.61) 



Since the single particle energies are independent of P, they only contribute a con- 
stant term to the dispersion energy. We renormalise the dispersion energies by 
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Figure 4.5: Plots of the various inter-LL dispersion relations described in the text 
corresponding to partial filling of the n = LL. 



setting ei + E^;i = 0. 

We have already seen that there are three types of transitions: those which 
are not mixed to any other transition, those which exist in a pair and those comprised 
of four mixed transitions. I denote the (renormalised) dispersions for these three 
types as E^^\E^'^^ and E^^^ respectively. The resulting dispersion relations are 
shown in Fig. |4.5[ The degeneracy of each dispersion relation depends upon the 
filling factor; for example there are no mixed pairs for the case u = 2. It is clear 
that 

^(1) = d, (4.62) 



which is given explicitly in Eq. (4.59). The mixed pairs of excitations are found by 
diagonalising the matrix 

(4.63) 

where the inclusion of ^ on the diagonal is due to the additional self energy pos- 




sessed by mixed pairs of transitions, as seen in Section 4.7.2 This yields 



(4.64) 



where the different letters signify different branches of the dispersion. For the case 
= — 1, the dispersions corresponding to the four exchange- mixed transitions are 
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found by diagonalising 



( d + ex ex ex 
ex d + ex ex 
ex ex d + ex 



-ex \ 

-ex 
-ex 



\ 



-ex 



-ex 



(4.65) 



-ex d + ex j 



This gives a triply degenerate eigenvalue which equals E^^^ and a non-degenerate 
one: 

(4.66) 



= d, e\^'^ = Aex + d. 



Although the matrices for different sublevel filling factors of the n = LL may look 



different to that in Eq. (4.65), they yield the same eigenvalues and degeneracies. 
In Fig. 4.5, in addition to plotting Ej^'^\ its direct {d = E^^'^) and exchange {Aex) 



contributions are plotted separately. 

We may now determine the excitonic band width for excitations with bare 
energies ei. For all the dispersions calculated above, their lowest value occurs at 
P = and they are bounded above by the P ^ cc limit. Hence, if only transitions 
with no spin or pseudospin flip are considered, the band has width 0.75ii^o- If 
sixteen transitions are considered however, it has width Eq. 



4.8.4 Experimental evidence 

There have been numerous experimental studies investigating the CREs in graphene 
[195 , 207, 224- 226] . Ref. [225] describes a photoconductance study in which modu- 
lations of the samples' conductivity upon radiation are used to determine the CRE. 
An asymmetry between the electron and hole bands was observed, suggesting a 
difference in the electron and hole velocities. This is interesting, since it is widely 
thought that graphene has perfect particle-hole symmetry. To the best of my knowl- 
edge, such an asymmetry has not been reproduced experimentally, though this may 
well be due to the high precision of measurements required to observe it. A pop- 
ular direct probe of the CRE is spectroscopy. Initial far infrared spectroscopy of 
ultrathin graphite [224] reported good agreement with the simple non-interacting 
Dirac picture. However, many experimental papers on monolayer graphene claim 



that their results indicate the presence of many body effects 195 207 2261. The bare 



transition energy and the many body correction both have approximately the same 
dependence on the magnetic field, so this cannot be used to measure the correction. 
The presence of many body effects cannot be demonstrated by examining a single 
resonance energy, since the effect of e-e interactions is just to renormalise vp and 
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B(T) 

1 4 9 16 25 




Figure 4.6: (Taken from Ref. [207]). Cyclotron resonance energies plotted as a func- 
tion of magnetic field, B. The transition Ti corresponds to the —1^0 transition 
(v = —2), which is equal to the 0^1 transition (v = 2). The transition T2 corre- 
sponds to the —1 2 transition (u = +2). The solid blue line is a best s/B fit to 
the Ti transition data, yielding vp = (1.12 + 0.02) x lO^ms^^. The dashed red line 
is a scaling of the solid blue line by a factor of 1 -I- -\/2. 



one cannot turn off the interactions and experimentally determine a bare vp for 
comparison. However, as we have seen in Section [4.7.2 the self energy correction 
for excitations with different bare energies, does not scale according to the ratio 
of these bare energies. Hence by studying more than one resonance it should be 
possible to see a signature of many body effects. 

This approach was used in Ref. [207] , where the energies of the —1^0 
and —1^2 transitions were compared at filling z/ = — 2. The results are shown 
in Fig. 4.6, which is taken from Ref. [207| . A renormalised Fermi velocity vp was 
calculated, by setting the experimentally determined energy of the —1^0 transition 
equal to V2hvp/£B- It was then noted that the experimental data for the —1^2 
transition energy lies slightly above the line ^/2hvp{l + ■\/2)/£b (red dashed line in 
Fig 4.6) and that this shift is outside the error bars, which are contained within 
the symbol size. These results are consistent with the theory discussed in Sections 



4.7 and 4.8, as I shall now explain. We need to calculate the contribution to the 



transition energy from interactions, assuming no spin or pseudospin flips and V = Q 
(by conservation of momentum, since the photon does not supply any momentum). 
Let's denote these interaction energies by -F-10 and -F-12 for the —1^0 and —1 2 
transitions respectively. They have direct and self energy components, since the 
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for both transitions. For —1 0, we have 



exchange energy vanishes at V 
F-io = d{0) + E'^\%^ K OMEo, where d is defined in Eq. (|4.59|). Using 



j~/d 2tiSi — 1t2S2 
2tisi —1t2S2 



rrd 21 , fjd 11 , fjd 20 , frd 10 
21 + 11 + 20 + ^ 10 



e 4 

f28 



(66 + 13p2 + + P^) lo ^ 



p2 

4 



(23 + + -p4^ 



(4.67) 



we have for -1 ^ 2, F_i2 = U''2rlsl-lrTs2 (0)+^"C2 ~ 2.17^0- Hence F_i2/F_io « 
2.7 > 1 + \/2 K 2.4 in quahtative agreement with the experimental data for the 
—1^2 transition lying above the red dashed line in Fig. 4.6 



2riSi 
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Chapter 5 

SU(4) symmetry in graphene 



Symmetry is ubiquitous in physics. The wave function of a set of bosons/fermions is 
symmetric/antisymmetric under the exchange of qiiantum numbers for any pair of 
particles. In this Chapter we focus on the spin and pseudospin degrees of freedom of 
CEs, considering them separately from the other quantum numbers. The two pos- 
sible spin projections and two possible pseudospin projections, yield four equivalent 
possibilities for the spin/pseudospin state of an electron in graphene, which leads to 
the SU(4) symmetry. We first classify the neutral CEs according to their total spin 
and pseudospin and then explore the SU(4) symmetry, by applying techniques used 
in elementary particle physics to classify hadrons. 

5.1 Classification of excitons by spin and pseudospin 

Consider a system of two electrons and suppose each one has spin and no other 
quantum numbers associated with it. Then it is well known that a basis for the 
possible states of this two electron system may be constructed in such a way that 
the basis states are eigenstates of the total spin operator and total spin projection 
operator. For CEs it is clear that the total spin and pseudospin projections, Sz and 
Tz, are always well defined quantum numbers. Explicitly for the CEs QI/j^/^J^}, 
-^]viA/2A/3 Ia*)' '^z = Si — S2 and Tz = ti — T2. We would like to see whether it is also 
possible to assign a total spin, S and total pseudospin, T, to CEs. To do this it 
is necessary to express the operators for these quantities in second quantised form. 
The operator for total spin, S^, can be derived using S"^ = + Sy + and 




(5.1) 
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where i = x,y,z and = {nmr}. We stick to the electron representation to simphfy 
the formulae; it is also the natural choice to determine the total spin of the ground 
state. After rearranging to normal ordered form, we find 

= — ^ (^4cJy|cJ^/|CAr'|CAr| - 2c'^^c'^, ^Cn' [CN] (5.2) 
N,N' 

N 

Note that is symmetric under flipping all spins, as expected and that it gives 
the correct value for the total spin of a single electron. The operator for total 
pseudospin, T^, is found analogously. It can be obtained from the expression for 
S^, by changing f, i^ffi JJ- N ^ N = {nms}. 

The total spin (pseudospin) of a CE is only well defined if the ground state, 
from which it is excited has total spin (pseudospin) zero. The total spin (pseu- 
dospin) of a ground state is found by expressing it in second quantised form and 
then acting on it by (T^). For example, 

i^=2>= fi n n n ^™io>- (s-^) 

Hence one can show that S* = only for u = +2 (or indeed any completely filled 
LL) and T = only for = 0, +2 (or any half or completely filled LL). 

We concentrate on neutral CEs from ground states with S = T = 0. In this 
case, S^Ql^j^Ji^) = [S\QI^j^J\j^) and t^gj^^^ji.) = [t^, Q]^^_^J|z.>. Choosing 
1/ = 2 (without loss of generality), one may show 

[S^<5LtOTt] = -4i^| 0rt("ll,"l2) (5.4) 

ft _ t t ^ 

^A^i '^Imi rt '^^t '^0m2 ri ^NI^lmirl'^N [ C0m2r f 



mi ,m2 



N 



and hence 



S^QLt Ortl^ = 2) = h^ 0rt(mi,m2) (clmirtCOmsrt " Am^r iC0m2T [) \v = 2>. 



mi,m2 



(5.5) 

A similar expression for S^Qirj, OtJ^ ~ ^) '^^^ obtained by flipping all spins in 
Eq. (5.5). Let's set vIitI 0rt("^i)"i2) = ^irj, OTj,("ii; ?7Z2) for all mi, m2 and r, and 
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let Insi T2S2) = Q\^^,^ ora^al'^ = ^^^^ 

ny + ut ^t> + 1 m to + 1 = (5.6) 

and we have shown that the CE, ^ (| 1lt> + Ut J|t> + I Hi Hi) + I Hi Hi» is a 
spin singlet (5 = 0). By similar reasoning it is also a pseudospin singlet (T = 0), so 
we define 



|5 = 0, T = 0; 5, = 0, = 0> = ^ (I Ht Ht> + I Ht Ht> + I Hi Hi> + I Hi Hi» 



(5.7) 



Further calculations find 



|5 = 0, T = 1; 5. = 0, = o> = ^ (I Ht Ht> - I Ht Ht> + I Hi Hi> - I Hi Hi» , 

(5.8) 

|5 = 1, T = 0; = 0, = o> = ^ (I Ht Ht> + I Ht Ht> - I Hi Hi> - I Hi Hi» , 

(5.9) 

|5 = 1, T = 1; 5. = 0, = o> = ^ (I Ht Ht> - I Ht Ht> - I Hi Hi> + I Hi Hi» • 

(5.10) 

Notice that the four states above all have Sz = 0, = 0. The states corresponding 
to different spin/pseudospin projections can be found by acting on these states by 
the spin/pseudospin raising and lowering operators. In second quantised form these 
are 

S+ = Sx + iSy = c^isf^CNi, (5.11) 

N 

S- = Sj;- iSy = cjy^CAri, (5.12) 

N 

with similar definitions for T^,T^. In total there are 16 states; they are listed in 
full in Appendix [C| An alternative way of classifying the states is according to the 
quantum numbers, 5, T, /, I^. The quantity / is the total (spin plus pseudospin) 
angular momentum projection: 

P\S, T, I; h) ={S + f f\S, T, I; h) = I{I + 1) \S, T, I; h) (5.13) 

and Iz = Sz + Tz. 

A similar analysis may be done for trions, where the states can be classified 
as spin/pseudospin doublets or quadruplets. However, this is considerably more 
complicated than for the two-particle case, since it involves 64 states. In addition, it 
doesn't explain the degeneracy in the energies of CEs. Hence instead of pursuing this 
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spin/pseudospin classification, we follow another route, which is able to predict the 
degeneracies of CE energies. To do this, I must first introduce some mathematical 
concepts. 



5.2 Lie algebras 

In this Section, I assume that the reader is familiar with the concept of a dis- 
crete group (see for example Ref. 227] ). An n-parameter continuous group con- 



sists of the elements A{a) = A{ai, a2, ■ ■ ■ , ctn) obtained by varying the parameters 
ai, • • • ! ttra continuously, together with a binary operation. By the closure re- 
quirement for a group, A{a)A{l3) = ^(7). Thus 7 = f{a,(3), where / is a con- 
tinuous function. If / is also analytic (it has a convergent Taylor series expansion 
everywhere in its domain), then the group is a Lie group [228 . Let's define 



B = V„^(a) |„=o, (5.14) 

from which it follows 

A{a) = e^-" = lim f 1 + . (5.15) 

n^co \ n / 



From Eq. (5.15) it is clear that the components of B generate the Lie group. They 
form the basis of a vector space, V, over some field, F. Moreover, one can prove 
that V together with the binary operation, [■,•]: F x y ^ 1/ form a Lie algebra, 
since they satisfy the axioms for all a,b e F and for all x,y, z e V: 

1. Bilinearity 

[ax + by, z] = a[x, z] + b[y, z], [z, ax + by] = a[z, x] + b[z, y] (5.16) 

2. Jacobi Identity 

[[x,y],z] + [[z,x],y] + [[y,z],x] = (5.17) 

3. Antisymmetry 

[x,y] = -[y,x]. (5.18) 

We have seen how a Lie algebra emerges from a Lie group. It is also possible to 
reverse this process and specify a Lie group from a Lie algebra; this is Lie's theorem. 
Let us illustrate these ideas on a simple example. SU(2) is the group of 2 x 2 
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unitary matrices with determinant one. Each element has the form 




Uia,b)=i^^^^ : I, (5.19) 

where a, 6 are complex numbers satisfying |ap + |5p = 1. Hence SU(2) is a Lie 
group with dimension 3. If we consider a spin ^ system, whose state is given in 

/ at \ 

the Pauli two-component formalism by I ) , a physical application of SU(2) is 

J 

/at \ 

that it consists of the matrices which when applied to I yield the new state 

V " / 

resulting from a particular rotation of space. The three independent operators Sx, 
Sy and Sz form the basis of a vector space for the Lie algebra. Notice that 

{Si,Sj\ = eijkSk, (5.20) 

where k e {x, y, z} and eijk is the Levi-Civita symbol. Hence the vector space is 
closed under commutation, as specified in the definition of a Lie algebra given above. 
The Lie group is generated by Sx, Sy, Sz or strictly speaking by the corresponding 
Pauli matrices ax, cry, az, in the sense that each U{a, b) = exp (— ^^^) , where n is 
the axis of rotation in 3D space and 9 is the angle of rotation. In fact, for any finite 
n, SU(n) is a Lie group of rank n — 1 with — 1 independent generators. 



5.3 SU(4) symmetry in graphene 

Graphene is SU(4) symmetric in the sense that the ordinary (non-interacting) Hamil- 
tonian commutes with the elements of SU(4). In addition, the Hamiltonian of 
the Coulomb interactions is also SU(4) symmetric up to small symmetry-breaking 



terms [199 . This is due to the equivalence of the two valleys and the two spin states, 
including the continuous interchange between spins and valleys. We assume here, as 
above, that the ground state corresponds to the n = LL being either completely 
full or completely empty. The four basis states that form a fundamental representa- 
tion of SU(4) are | Ht), | Jj-t); I 1li>> I J]-i>- I shall label them in analogy with the four 
flavours of down, up, strange, and charm quarks, by {JJ-i, Hi, JJ-t) Ht} = {d, u, s, c}. 

The generators of SU(4) can be constructed using bilinear combinations of 
the fermionic creation and annihilation operators conserving the number of electrons 
and holes. We have 

Cij = ^ (^nmi^nmj — ^ d\^^dnmi, (5-21) 



nm 
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where i,j e {d,u,s,c}. The generators satisfy the commutation relations 



[Cij,Cki] = SjkCii — 6iiCkj, (5.22) 

so are closed under commutation. The operators of spin and pseudospin can be 
expressed in terms of the Cij. For example, Tz = ^ (Ccc + Cuu — Cgs — Cddj stnd 
T+ = Tj; + iTy = Cud + Ccs- The elements of SU(4) may then be expressed as 
U = exp (^i^jij ^ij^ij^ ^ with @ij the transformation "angles' 
commutes with the Hamiltonian, applying it to an eigenstate of the Hamiltonian, 
yields another eigenstate of the Hamiltonian with the same energy. This means that 
the eigenstates of the Hamiltonian exist in multiplets such that applications of the 
SU(4) operators to any state always result in another state in the same multiplet. 
In addition, all states within the same multiplet are degenerate. It is possible to 
calculate the multiplicities of our states (CEs) and also to learn something about 
their symmetry, upon interchanging the quantum numbers of pairs of particles, 
before recourse to numerical methods. This is done using Young diagrams. 



199 
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Since U 



5.3.1 Young diagrams 

Young diagrams are a useful technique for determining the SU(n) multiplet structure 
of a composite system of particles. In the following brief description, I treat them as 
a useful tool, stating several results without proof. For a more rigorous development 
of the theory, see Ref. [230| . 

Young diagrams consist of n or fewer left-justified rows of boxes, such that 
each row is as least as long as the row beneath it. Roughly speaking, each box 
represents a particle. Each SU(n) multiplet is uniquely labelled by a list of n — 1 
non-negative integers, (ai, 02, . . . , a„_i). The integer in the i^^ position corresponds 
to the number of boxes in the i*^ row minus the number of boxes in the {i + 1)*^ 
row i.e. the overhang of the i*^ row. Some examples of Young diagrams for different 
SU multiplets and their labels are given in Fig. |5.1[ By the definition of a multiplet 
label, adding columns of length n to the left of a Young diagram is redundant. Young 
diagrams can be used to calculate the multiplicity of a multiplet. The multiplicity, 
[A^(n)] of an SU(n) multiplet is 

n— In— 1 . _ . -. 

«=1 J^l 
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SU(3) SU(4) SU(5) 

(3,1) (0,1,1) (0,0,1,0) 



[24] [20] [10] 

Figure 5.1: Examples of Young diagrams for different SU multiplets together with 
their unique multiplet label and the corresponding multiplicity. 



Thus for the n = 4 case, 

rjYC4^1 - («i + 1) ("2 + 1) (Q3 + 1) (qi + a2 + 2) (02 + 03 + 2) (ai +02 + 03+8) 
^ ^ ~ 1 1 1 2 2 3 

(5.24) 

The boxes in Young diagrams can be labelled by positive integers to represent 
the quantum state of a particle, according to certain rules. These are i. in each row 
the labels must not decrease, as read from left to right, ii. in each column the 
labels must strictly increase as labelled from top to bottom. The Young diagram 
then represents a state, which is first symmetrised with respect to exchanging pairs 
of quantum numbers for boxes in the same row and then antisymmetrised with 
respect to exchanging pairs of quantum numbers for boxes in the same column. 
Of course, for bosons the state should always be completely symmetric and the 
Young diagram a single row and for fermions, the state should always be completely 
antisymmetric and the Young diagram a single column. However, this holds if all 
quantum numbers that can label a particle are used. Often this is not the case, for 
example, spin quantum numbers or in the case of graphene, spin and pseudospin 
quantum numbers, are treated separately to the orbital quantum numbers. Under 
these circumstances other symmetric patterns may emerge. 

Let's illustrate this on the familiar example of a two-particle system, where 
each particle has spin ^ and two possible spin projections, which I shall call 1 and 
2. The first step is to create all allowed Young diagrams consisting of two b oxes. 

As we saw in Section 



5.2 



and multiplicity [3]; the latter 



In fact there are only two possibilities, and 
these are SU(2) multiplets. The former has label (2 
has label (0) and multiplicity [1]. Following the labelling rules, the Young diagram, 
can be labelled as 1 I 1 I, 1 I 2 or 2 I 2 |. Since the state should be symmetric 



with respect to permuting the quantum numbers of the two particles, we deduce 
that the states are |11>, ^ (|12>+ |21» and |22>, where \ij} means particle one is 
described by the quantum number i and particle two is described by the quantum 
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number j. We recognise these as the symmetric spin triplet states. The only way 

T 



of labelling 



IS 



Since the two boxes are now in the same column, the state 



_2] 

must be antisymmetric with respect to interchanging the two particles' quantum 
numbers, which is satisfied by (|12) — |21)). This is the familiar spin singlet. 

There are two fundamental representations of SU(n): [n], with Young di- 
agram and label (1,0, ...0) and [n] with Young diagram composed of n — 1 



vertically stacked boxes and label (0 . . . , 0, 1). In general, whenever two multiplets 
have the same entries in their labels, but in reverse order, they are said to be con- 
jugate and one of their multiplicities is denoted with a bar. The latter represents 
an antiparticle; this makes sense if we identify n vertically stacked boxes with the 
vacuum. I shall also use the symbol, 1 to denote the vacuum. For the current 
purpose, the hole in a CE behaves as an antiparticle. In order to create SU(n) 
multiplets corresponding to e-h and e-e-h complexes, we need to be able to combine 
two Young diagrams correctly. There is a standard way for doing this, which I shall 
now explain. The general idea is to take the first Young diagram as being the upper 
left portion of the emerging combined diagrams and label the second diagram such 
that all boxes in the first row are labelled 1, all boxes in the second row 2 and so 
on. The boxes from the second diagram should then be added to the ends of rows 
of the first diagram or the bottom of the first diagram in all the different possible 
ways subject to the following rules: 

i. The resulting diagram takes the form of an "allowed" Young diagram as previously 
described with a maximum of n rows for SU(n) multiplets. 

ii. Each number should appear no more than once in any column. 

iii. The sequence of numbers formed by reading from right to left in the first row, 
then the second row etc. is admissible. I define a sequence as being admissible if at 
any point in the sequence for any integer j, iVj ^ Nj holds for all i ^ j, where Ni 
is the number of times the integer i has appeared in the sequence up to that point. 
This is best illustrated by an example. Two SU(3) octets may be combined in the 
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following way: 
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To obtain the final line of Eq. (5.25), I removed 



0[2] (5.25) 
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'01. 



from the left of Young diagrams 



where possible. Eq. (5.25) can also be written in terms of multiplet labels as 



(1, 1) ® (1, 1) = (2, 2) 8 (3, 0) 8 (0, 3) 8 (1, 1) 8 (1, 1) 8 (0, 0) (5.26) 
and multiplicities as 

[8] (X) [8] = [27] 8 [10] 8 [10] 8 [8] 8 [8] 8 [1] • (5.27) 
We are now in a position to work out the multiplet structure for CEs. 
5.3.2 SU(4) multiplet structure for excitons 

Our neutral CEs can be thought of as a particle-antiparticle (e-h) complex where 
particles can have one of four flavours, as described previously. They are thus 
in complete analogy with mesons, which are elementary particles composed of a 
quark- antiquark (qq) pair, provided the quarks are either first (u, d) or second (c, s) 
generation. We find the multiplet structure for excitons by combining the Young 
diagrams for an electron and hole as follows: 



1® 



1. 



(5.28) 



Eq. (5.28) is equivalent to (1,0,0) ® (0,0,1) = (1,0,1) 8 (0,0,0) or [4] ® [4] 



[15] 8 [1]- In terms of the classification of excitons described in Section 5.1, the 



vacuum singlet corresponds to the spin and pseudospin singlet in Eq. (5.7). It can be 
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Figure 5.2: SU(4) [15]-plet describing the excitonic states in terms of flavour isospin 
(S^), liypercharge (Y) and charm (c). 



expressed diagrammatically as 



and symbohcally as Qo = 5 {dd + uu + ss + ccj . 



From Eq. (5.28) we can deduce that all exciton energies will have degeneracy 1 or 
15. This is useful, since it reduces the numerical effort. For the case of a completely 
filled LL, we need only calculate the energies of states where there is no spin or 
pseudospin flip. We have in fact already seen these degeneracies in the excitonic 
dispersion relations calculated in Section |4.8.3[ Note that for u = 2 there are 12 
unmixed 0^1 transitions and 4 with no spin or pseudospin flips, which are mixed 
by the exchange Coulomb interaction. There are of course no pairs of transitions 
(one 0^1 and the other —1^0), which are mixed by the direct Coulomb 
interaction for this filling. The possible dispersion relations are then E^^\ which 
has degeneracy 12 + 3 = 15 and Ej^'^\ which is non-degenerate. 

A way to visualise the [15] multiplet is using the standard elementary par- 
ticle physics SU(4) representation |229| , as shown in Fig. 5.2 Namely, the three 
orthogonal axes show: the flavour isospin projection T,z = ^ (Cuu — Cdd^ , the hy- 



percharge Y 



Cuu + Cdd + Co 



\Css and the charm c = Ccc- For example, 



the exciton dc (exciton where the electron has flavour d and the hole flavour c) has 



2' 



y = and C 



-1. 
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Figure 5.3: (a) Quartet ([4]), (b) [20] and (c) [36] SU(4) multiplets describing the 
trion states. States are plotted in terms of flavour isospin (Sz), hypercharge 
(y) and charm (c). Shading denotes the different levels of charm. 

5.3.3 SU(4) multiplet structure for trions 

In contrast to the exciton case, the charged CEs or trions have no direct particle 
physics analogy, since bound structures of the form qqq or qqq have colour and so do 
not exist. The closest thing is a baryon, which is a bound structure of three quarks. 
However, we may still use the same techniques. For the trion, the multiplet 
structure is: 



(X) (X) 




(5.29) 



5.3 



or [4] (X) [4] ® [4] = [36] 8 [20] 8 [4] 8 [4]. The multiplet plots are shown in Fig. 
It is not clear what the physical meaning is behind the presence of two identical 



Young diagrams in Eq. (5.29). 
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The multiplet structure for the X+ trion is 



® 



01 



(5.30) 















or [4] ® [4] ® [4] = [36] [20] [4] [4]. Although the Young diagrams in Eqs. 
( |5.29 ) and (5.30) look different, the multiplicities are the same. The multiplet plots 



for the X'^ trion look similar to those for the shown in Fig. 5.3, except that 
all electrons are replaced by holes and all holes by electrons, which results in each 
coordinate of the three-particle complex changing sign. Indeed we shall see in the 
following Chapter that a symmetry exists between positively and negatively charged 
trions. 
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Chapter 6 

Charged collective excitations in 
pristine graphene 



The possibility of the formation of charged collective excitations (charged CEs) or 
trions in pristine graphene has already been mentioned in Section [4.2[ The present 
Chapter expands upon this by applying the techniques detailed in Chapters |4] and 
[Hand examines the results. 

6.1 Techniques for charged collective excitations 
6.1.1 Creation operators 

We saw in Section [44] that charged CEs can be labelled by the quantum numbers Jz 



and k. Eq. (4.6) can thus be rewritten to give the creation operator for a negatively 



charged CE or X particle as 



m,i,m2 



Notice that the sum over has been removed, since it is now fixed by mi,m2 and 
Jz according to Jz = \n-i \ — 'nil — |?^2| + "i2 + \n-i\ — ms — ^. Similarly the creation 
operator for a positively charged CE or particle is 



For the , Jz = l^^il — itii — \n2\ + m2 — \n^\ + ms + 
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Figure 6.1: States resonantly mixed by the Coulomb interaction for various filling 
factors, V. Filled circles (•) denote electrons, empty circles holes (o) and shaded 
regions, the Dirac sea. The ringed particle structures are bound trions. (a) v = 
—2 + £. The electron in the n = 1 LL is only mixed for angular momenta Jz ^ \- 
(b) u = 2 — e. The hole in the n = —1 LL is only mixed for angular momenta 
Jz ^ (c) u = 2 + e. The electron in the n = 4 LL is only mixed for angular 
momenta ^ |. 











C V, 














1 



Figure 6.2: Photon polarisation operator tt(uj). The box shows the resonant exci- 
tation of the e-h pair with the de-excitation of the electron to n = LL. Because 
of the multiple Coulomb interactions (wiggly vertical lines) the intermediate state 
becomes a bound trion Xqq_^, with two electrons in the n = LL and a hole in the 
n = —1 LL, leading to a singularity in 7r(a;). 



6.1.2 Mixing of charged collective excitations 

We focus on charged CEs where there is no spin or pseudospin flip in the initially 
excited e-h pair, since these are most likely to be optically active. As for the neutral 
CEs, we only take into account mixing between resonant excitations. Fig. |6.1| shows 
examples of which trions are mixed for different filling factors. Surprisingly, a single 
particle state may be mixed to the three particle complexes. This is an electron for 
filling factors of the form = fi + e and a hole for filling factors of the form u = fi — e. 
Looking at Fig. |6.1| [^a), one can think of this as being due to the recombination of 
the electron in the n = LL and the hole in the n = — 1 LL, which allows the second 
electron to be promoted from the n = to n = 1 LL, since the energies are the 



same. Such a process is illustrated in Fig. 6.2 For the case u = 2 + e shown in Fig. 
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|6.1| [c), an electron from the n = 4 LL may become mixed, since £4 — 61 = ei — eo- It 
is important to note that a charged CE wiU only have a non-zero amplitude for the 
single particle state, provided the generalised angular momentum projection for the 
single particle, jz, is the same as that for the three particle complexes, Jz- Suppose 
the single particle has orbital quantum numbers n, m. Then for the , an electron 
will only be mixed to the three particle states when jz = |n| — m — ^ = Jz- Hence it 
is only mixed for Jz ^ |n| — ^. Similarly for the , a hole will only have a non-zero 
amplitude when Jz ^ 5 — |n|. 

Once it is understood which states are mixed (significantly) by the Coulomb 
interaction, the energies and coefficients, (mi, 7712) or i?j^^^^_^^(mi, m2), 

may be found by diagonalising the matrix for the Hamiltonian, given by Eq. ( |4.7| ) 
with some additional terms to describe the scattering to a single particle state. In 
our calculations, we typically used ^ mi,m2 ^ 15 which involves diagonalising 
square matrices with between 600 and 900 rows. 



6.1.3 Self energy corrections 

Mostly, the self energy correction is easy to handle, since it is the same as that for 
neutral CEs. Let's first consider the case u = —2 + e shown in Fig. |6.1[ a). The 
excitation is a mixture of four states of the form -Rqts -Its offt I ~ ^^'^ case, 
Jz ^ 5, an additional single electron state, c| 1 | — 2). Note that the ground 

state is Cqa|^| — 2) for some general m, not simply | — 2). Then, considering the 



single particle state, the self energy contribution to {—2\c^^i_j^Hc'^^^i ^ | — 2) — 
(— 2|coff-t-m-ffcQ||,|^| — 2) is E'q^I — O.TSE'o- The self energy correction to 

(-2\Rors -Irs 0^t^<. -Irs O^fl " 2> " (''^l^mmH cl^^J " 2> is i^^rl = E%%,, aS 

for neutral CEs. 

There are a couple of cases, where calculating the self energy corrections 
is more involved. These are v = 2 + e for the and v = —2 — e for the 
X^ . As we shall see in the following section, these two cases are dual, so it 



excitation is a mixture of four states, ^i^s Ors ifi-f 1^)' case, Jz ^ - 



is enough to consider only the former, which is shown in Fig. 6.1 c). The X 

7 

2' 

an additional single electron state, c'^^7_j \ — 2). The ground state is c|||,|^|2). 
The self energy correction to (2\Rirs Ors ii\]HR\^^ ifrtl^^ ~ ('^\'^n\rnHc\^^^\2} is 
i?Orl' before. However, the self energy correction to {2\c^^7_j^Hc'^^^t_j |2) — 
<2|ci^l^ij4^^„|2> is 2"Lo[^SE (4, -n) - E^e (1, -n)]. Using n, * 810 at B = 20T 
as before, j:f=o[EsE (4, -n) - Ese (1, -n)] * 0.89J;o- 

As previously mentioned, the charged CEs appear as discrete states below 
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the neutral CE continuum discussed in Section 4.8 The binding energy of a charged 
CE is the difference between the energy of the lower continuum edge and the energy 
of the charged CE. Thus in the graphical presentation of numerical results for the 



energies (see Section 6.3), the energy of the lower continuum edge is set equal to 
zero. For the cases /x = 0,+l, it is possible for an exciton to be created with 
both electron and hole in the n = LL and for the second electron to be in the 



n = 1 LL. Using results from Section 4.8.3, this gives a lower continuum edge at 
the energy ei + Eq^^ — Eq. For the cases fi = +2, the lower continuum edge is at 
ei + -Eq" — 0.75Eo, since only excitons with one composite particle in the n = LL 
and one with |n| = 1 can be formed. 

6.1.4 Symmetry between X"*" and X 

It is only necessary to perform calculations for either positively or negatively charged 
CEs, since they are linked by a symmetry relation. The energies and expansion 
coefficients are identical for an state formed at filling factor v = fi — e with 
quantum numbers Jz, k and an state formed at filling factor u = —fi + e with 
quantum numbers — Jz, k. The duality is an immediate consequence of the inherent 
particle-hole symmetry in pristine graphene. An example of this can be seen in 



Figs. 6.1 a) and (b). From now on, the discussion will focus on the X particle. 



6.2 Optical selection rules for charged collective excita- 
tions 

Optical selection rules for charged CEs can be derived from the commutation rela- 
tions: 

[6'H+,k^] = 0, (6.3) 
[6n+Jz] = +6n+, (6.4) 



where 671+ is defined in Eq. (4.21 ). Suppose the ground state is |e) = CnrsmlfJ'} with 
jz = I'^l ~ ""^ ~ 5 ^i-iid the photocreated trion state, |ee/i), has geometric quantum 
numbers k, Jz- The intensity of the resulting absorption peak is then directly 
proportional to the square of the dipole matrix element, {eeh\5'H+\e}. Thus from 



Eq. (6.3), we may derive the optical selection rule, k — m = 0. This reflects the 
negligible linear momentum of a photon in the dipole approximation and the fact 
that k is associated with translations; it is the discrete analog of linear momentum for 



charged states in a magnetic field. From Eq. (6.4), Jz — jz = +1 for the a- circular 
polarisations. Combining the selection rules for k and Jz yields Jz+k = |n| — ^ + 1. In 
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Figure 6.3: Bound X states and their multiplicities as a function of Jz for (a) 
= —2 + e and (b) i' = 2 + e. Only the seed states (A; = 0) are shown. Energy zero 

0.75£;o- 



corresponds to the ORE at ei + Eq^^ 



Section 4.4 , we saw + A; is equal to the generahsed angular momentum projection 



of the seed state, Hence if a family of states is optically active in the 

polarisation, ji^-* = |n| + ^; if it is active in the polarisation, ji^^ = |n| — |. 

These selection rules are a necessary, but not sufficient condition for a family 
of states to be optically active. There is another selection rule associated with 
the electron flavour discussed in Chapter Sl It arises from the fact that dTi+ccQo 



(see Eq. (4.21)), so the photon is flavourless in the SU(4) sense. From this and 
considerations of the dipole transition matrix elements, (ee/i|(5o|e), it immediately 
follows that the states in the [20] and [36] multiplets are all dark, while the states 
in the [4]-multiplet may be bright, provided they have the proper orbital quantum 
numbers just described. 



6.3 Numerical results 

We performed calculations for trions with no spin or pseudospin flips for the filling 
factors = /X + e with fi = 0, +1, +2 and charged CEs were only found below the 
band for = +2. This can be attributed to the fact that the lower continuum edge 
for the fi = +2 cases is 0.25-E'o higher than for the remaining filling factors, as seen 



in Section 6.1.3 Hence less e-h attraction is required to push the energy of the trion 
below the band. For the SU(4) symmetric case, /U = +2, the trion states with spin 
or pseudospin flips have the same energies as those with no flips. This is not the 
case for fi = 0, +1, but it is unlikely that charged CEs will be found even if spin and 
pseudospin flips are taken into account, due to the continuum reaching low energies. 
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Figure 6.3 shows the energies of states relative to the lower continuum 
edge for (a) u = —2 + e and (b) u = 2 + e. The relevant mixed states for these 
filling factors are shown in Figs. |6.1[ a) and (c) respectively. The multiplicities have 
been marked on, indicating which SU(4) multiplet the states belong to. Calculations 
were made assuming no spin/pseudospin flips in the exciton and that the additional 
bound electron had flavour c. Hence the numerical diagonalisation found the state 
Qqc from the [4]-plet, two states combining cuu, cdd and ess from the [20]-plet 
and three states combining cuu, cdd, ess and ccc from the [36]-plet. The remaining 
states in a multiplet can be found by repeatedly acting on a single state by the 



generators. The structure of the Young diagrams (see Eq. (5.29)) should tell us 



about the symmetry of the states they represent, although this is less clear when 
a hole or "antiparticle" is involved, so that a single box no longer corresponds to 



a single particle. The diagrams for the [20]-plet, 



and the [36]-plet, 



both have a hole to the left and two electrons to the right. For the [20]-plet the 
electrons are in the same column and states within this multiplet are antisymmetric 
with respect to interchanging their flavours. For the [36]-plet the electrons are 
in the same row and states within this multiplet are symmetric with respect to 
interchanging their flavours. 

Out of all the states formed at /i = +2, only those at one energy, indicated by 
the red circle in Fig. |6.3[ a), are optically active. Notice that these states obey both 
the flavour selection rule (they are in the [4]-plet) and the geometrical selection 
^ for the cr"*" polarisation. The state in Fig 



7(0) 



6.3 



'b) with J; 



(0) 



rules, J; 

satisfies the geometrical selection rules, but not the flavour selection rule, so is not 
optically active. 



6.4 Experiment 



Charged CEs of the 2DEG created by confining electrons to a semiconductor quan- 
tum well have been observed in many optical experiments using absorption and lu- 

. BothX+ 



minescence techniques 
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and X 
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231 



232 



234 



235 



particles have been detected in the presence and absence of magnetic fields. Sig- 



natures of charged CEs have also been seen in semiconductor quantum dots 234 



The binding energy is of the order of meV. Theoretical predictions [220 
should be 



suggest it 



O.li^Oi which agrees with experiment, since Eq ~ lOmeV. The binding 
energies we predict for charged CES in graphene are also ~ O.IEq, so they are sim- 
ilar to those for the 2DEG in units of Eq. However, Eq depends inversely on the 
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Figure 6.4: Predicted absorption spectra in pristine graphene for the briglit state 
at energy — 0.09ii^o with = — 2 + e. The inset schematically shows the relevant 
optical transitions. Energies are given in units of Eq with zero energy corresponding 
to the CR energy, ei + - 0.75£;o. 



dielectric constant, e. For GaAs quantum wells, e ~ 12. It is difficult to assign a 
value to e for the case of graphene, since this will depend on its environment. At 
least for suspended graphene, e < 12, meaning that trion binding energies should 
be higher in suspended graphene than in the 2DEG in GaAs quantum wells. Thus 
trions should be easier to observe in graphene, at least in principle. 

A recent ARPES study of potassium doped graphene in the absence of a 
magnetic field, indicates that the simple picture of two conical bands touching at 



the Dirac point is not the full story 236 . In addition to these two charge bands, 
a pair of plasmaron bands are clearly visible in their high resolution data, so that 
the crossing at the Dirac point is in fact resolved into three crossing points. A 
plasmaron is formed when a hole becomes bound to a plasmon. Hence it is the 
zero-field analogue of our charged CEs. 

To the best of my knowledge, charged CEs in graphene in high magnetic 
fields have not yet been observed, which is perhaps surprising since they are now 
routinely seen in semiconductor quantum wells. Our study found one bright 
state at u = —2 + e with binding energy k 0.09£'o {Jz = \, active in cr"'", see 
Fig. |6.3[ a)). Note that the eeh complexes involved are created when an exciton 
with the hole in the n = —1 LL and the electron in the n = LL become bound to 
an additional electron. As such they are not optically active in the a"'" polarisation. 



However, as indicated in Fig. 6.2, they "dress" the optical response of the resonant 
process by which an electron is promoted from the n = to n = 1 LL, which is 
active in the o""*" polarisation. This should be observable in absorption experiments, 
provided the noise is low enough. It is expected to give rise to a separate absorption 
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peak below the CRE or may show up as a low-energy CRE shoulder, depending on 



the broadening. This is shown in Fig. 6.4 , where the spectra are broadened ~ 0.05£'o, 
which amounts to 2 meV at = 10 T. This broadening is slightly below that found 



in absorption spectra of the cleanest graphene samples currently available 195 



which may explain why signatures of charged CEs are not seen in this data. The 
intensity of both the peak and the main CRE peak scales according to e in 
the cr"'" polarisation. Notice that because of the e-h symmetry, there exists the X"*" 
counterpart at i/ = 2 — e, which is bright in cr^ polarisation. 

The dark bound and X"*" states may in principle be observed by tunneling 
spectroscopy, which is a sensitive tool for probing discrete states in the spectrum. 
This can be pursued, e.g., in tunneling experiments involving gate-tunable graphene 



quantum dots [237 . Also, the various symmetry-breaking terms such as external 
fields, disorder and lattice defects, ripples, and deformations p| |103 may partially 
lift the limitations following from the orbital and SU(4) selection rules. This can 
make some of the dark states "grey" and detectable in the absorption spectra as 



well as possibly in photoluminescence 238 
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Chapter 7 



Localised neutral collective 
excitations in dirty graphene 



The previous three chapters have studied CEs of pristine graphene. Although 
graphene is relatively clean, when compared to other materials, disorder is still 
always present. For the purpose of applications, there is a lot of interest in deter- 
mining the dominant scattering mechanisms, which limit the DC conductivity of 
graphene. One would expect the type of disorder in a graphene sample to depend 



on its method of fabrication 74 , 75 . Although mechanical exfoliation produces high 
quality samples, useful for investigating fundamental properties of the material, scal- 
able growth methods such as chemical vapour decomposition and epitaxial growth 
are the techniques likely to produce large enough samples for industrial applications. 
Whenever graphene is grown on a substrate, it will experience strain due to the mis- 
match between the crystal lattices. This may result in ripples or blisters forming in 
the sheet, or alternatively the breaking of cr-bonds and the formation of structural 
defects such as pentagons and heptagons. In samples of graphene grown on nickel. 



lines of structural defects composed of pentagons and heptagons were observed 239 
They act as dopants resulting in a line of charge. Other sources of disorder may 
be point-like impurities either in the substrate or on top of the graphene layer. In 
this Chapter we study the localisation of neutral CEs on a single such impurity, 
which models the physical situation of a low impurity density. Two different im- 
purity types are considered for comparison: i. a short range ((5-function) scatterer 
due to a vacancy or neutral foreign impurity at one of the graphene lattice sites, ii. 
a long-range charged impurity on or close to the graphene layer. The position of a 
Coulomb impurity within a graphene unit cell is irrelevant due to its long-ranged 
potential, which barely changes over such length scales. A similar study was carried 
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out for the 2DEG for the case of the Coulomb impurity [240 



7.1 Impurity Hamiltonian 

7.1.1 Single particle impurity Hamiltonian 

The impurity is assumed to be at the origin and the symmetric gauge is chosen. 
We use the single particle Hamiltonian derived in Ref. |217| to incorporate the 
impurity potential. The amplitudes of the electronic wave function at positions Ra 
in sublattice A and in sublattice B are written as 



V'^(Ra) = e^'^-'^^FfiKA) + e^^'-^^Ff (R^), 
V'b(Rb) = -a;e^K-^^Ff (Rb) + e^^'^^Ff' (R^), 



(7.1) 
(7.2) 



where lo = exp(27rz/3). Here the -Fs are envelope wave functions which change 
slowly on the scale of the lattice constant, a. Upon substituting these into the tight 
binding equation and considering nearest neighbour hopping only, the effective-mass 
equation is found to be: 



vpcr.tl 
VFa*.U 



+ V 



E 



where 



V 



( i'A(r) 
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t'A(r) 
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VB{r) 
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Ff 
Ff 



(7.3) 



(7.4) 



and E is the single particle energy. In Eq. (7.4), 



VA{r) = 25(r-RA)fA(r) 



and 



Ra 



i(K'-K).RA„/ 



VAi^), 



(7.5) 



(7.6) 



where f^(r) is the onsite potential energy and g{r) is a normalised real function, 
which only has an appreciable amplitude for |r| less than a few times a. Similar 
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expressions for i'_B(r) and u_B(r) are obtained by replacing ^4 ^ i? in Eqs. (7.5) and 

Let us now examine the form of the impurity matrix, V, for the Coufomb and 
5- function impurity potentials. The Coufomb impurity potential is Vc(r) = 
where Z is the charge of the impurity in units of e and eimp is the dielectric constant 
associated with the electron-impurity interaction. Screening of Coulomb impurities 
has been seen to be relevant in graphene 241 , although the situation is not yet 
fully understood, particularly for a strong magnetic field. Coulomb impurities with 
charge \Z\ = 1 belong to the subcritical regime and hence screening effects due to the 
substrate and the electron system in graphene can be modelled via an effective charge 
\Zcs\ < \Z\ [241] , which is equivalent to assuming an effective dielectric constant, 
einip- Since energies are always measured in units of Eq here, which contains e, it 
is only the ratio e/eimp, which is important. Its value for graphene is not known 
, so for most results presented below e/eimp 
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1 and Z = +1 is assumed. 
Based on the previous discussion, the Coulomb impurity is long-ranged, so 
that VA(r) = vb{y) = Vc(r). In addition, for any long-range impurity, VA(r) and 
{'^(r) are negligible, since they are reduced to a sum over oscillating functions (see 
Eq. (7.6)). Hence for the Coulomb impurity, V*^ = Vcir)!^. For a 5-function 
impurity on an A sublattice site at the origin, fA(r) = ^'^(r) = Vo<5(r) and ^^^(r) = 
t?B(r) = for some constant Vq = Wa^, where W is the impurity potential. This 



gives 



I 5(r) 

5{v) 

\ 



5{t) 



5{r) 










(7.7) 



Conversely for a 5- function impurity on a i? sublattice site at the origin, w_B(r) 
^^(r) = Vo5(r) and va{y) = va{y) = 0, yielding 



/o 







5{r) 


-uj5{v) 










\ 



5{v) ) 



(7. 



Notice that the Coulomb impurity does not scatter between the valleys or sub- 
lattices, whereas the 5- function impurity does scatter between the valleys, whilst 
preserving the sublattice. 
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7.1.2 Second quantised impurity Hamiltonian 

The full second quantised Hamiltonian describing the single particle energies, the 
e-h interactions and the interaction with the impurity is H + Hirap, where H is given 
by Eq. (4.7) on page 66 and 

M,M' rn,m' AfJ^' m,m' 

The impurity matrix elements are given by Vij^lf^ = Sgg' ^ldr'^^^j^i{r)V^^Ni^), 
where i 6 {A,B,C} denotes the impurity type and = {nmr}, J\f = {tits} as 
before. For the Coulomb impurity, 

^cj/m = ^AfAf'^mm'an{'SnV\n\-lm + y\n\m) , (7-10) 

where Vnm = Hd^'^'iPnm{^Wc{f)'^nmir) are the 2DEG matrix elements calculated 
in Appendix [bJ We ignore off-diagonal terms with n ^ n' oi m ^ m' , since these 
give the second order contributions to the energy, which are small compared to the 
first order contribution. 

The general form of a matrix element for a (5-function impurity is 

VAiB)'^::'^s,,v^rk^, (0) i^i^ (0) , (7.11) 

where le {|n|, |n| — 1} and ke {\n'\, \n'\ — 1}, with exact values set by r, r' and whether 
the impurity is on the A oi B sublattice. This imposes selection rules determining 
which states are connected by the delta impurity, since ijjnm (0) 7^ if and only if 
n = m. For completeness tpnn (r) = ^^/^-^ -/^n (^'^^ ^^P (i^)' where the L„ are 
Laguerre polynomials, so that Vnn (0) = Exact expressions for '^a{b)^^ 



uciic pui^iiuiiiiais, su iiiai lynny^) — ^ ■ i^^aui CApicasiuiis iui y A[B' 

are: 

2^ 



VaSI'"' = -^V^itir' = 5„2-|(2-^0"-^o.')j^^,5|„,|^,5|„|_^^y„!_, (7.13) 
V^il'"' = = 5.'2-^(2-^0"-V^')5..,<5|^„,|_i„,5|„|„Vb^, (7.14) 

= VB^i""' = 2-^(^-^--^°"')<5..'5|n'|^^'5|n|„.^0^. (7.15) 
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We treat both impurity types as a perturbation, using the first order term 
only. For the Coulomb impurity, the ratio of the second order energy correction, 
E^'^\ to the first order energy correction, E(^\ is ~ aS. For the case 

of the (^-function impurity, let us illustrate the details on a particular example, 
which works similarly to other scenarios. We shall consider an impurity on the B 
sublattice and calculate the ratio of the first and second order corrections to the 
single particle energy of the |0 HJ, 0) = cJ||,j^q|0) state. The first order correction is 
E(^^ = <0 Hi OliJimplO Hi 0> = The second order correction is 

^(2) = g I; l<n Hi m|gi^p|0 Hi 0>P + \(n U m|ffi^p|0 Hi 0)P 

n=l m=0 ^0 - en 

This represents virtual processes where an electron spontaneously hops to a higher 
LL and then decays back to its original state. Using (n Hi rn\Hiuip\0 Hi 0) = 
^Boll^ = ^\n\m^^^ and <n Hi m|^imp|0 Hi 0> = -a;(5|„|_i„^;^^^, we obtain 



for Eq. (7.16 



T/2 1 

= — ^ — y 1 . (7.17) 



This gives the ratio of the first and second order terms as 



EW 2V2TThvF£B ^iVn' ' 



It may seem disappointing that ~ Vi?, since we are in the high magnetic field 
regime, but the ratio is still small. In fact the ratio for the n}^ order correction 
behaves as ~ (tsJI^)" ' with ^ « 1. Using B = 20T, n, = 810 and 

Vo = leV X a , we find j^^^j k 0.03. This suggests that second and higher order 
terms in the perturbation expansion can be safely ignored. 

7.2 Mixing of excitations by the (5- function impurity 
Hamiltonian 

The Coulomb impurity matrix is diagonal, so it does not mix different neutral CEs, 
'3jviA/'2 Jz 1^)' Hence the mixing of excitations is determined by the e-h interaction. 



which was described in Section 4.6 The 5-function impurity behaves quite differ- 



ently in this respect, scattering between single particle states in different valleys. 



'^Note that V!"!, is the generahsed Harmonic number of order Uc of 1/2. 
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Figure 7.1: Set of mixed transitions for impurity on B sublattice, polarisation 
and v = —1. 



This results in mixing between excitons with and without a pseudospin flip. We 
seek optically active states, so the procedure is to begin with an optically active 
transition and then determine which other transitions it is mixed to. Recall for 
the Coulomb impurity, the mixed transitions were those four with no spin and no 
pseudospin flips and with Jz = +1. For the (5-function impurity, either six or eight 
transitions are mixed depending on the sublattice the impurity is located on, the 
filling factor and the light polarisation. Notably the 5-function impurity mixes states 
with different values of Jz, so Jz is not a good quantum number in this case. 

To make the approach clearer, we now consider a specific example, where the 
impurity is on a i? sublattice site (at the origin), the filling factor is = — 1 and 
the graphene sheet is illuminated by cr"*" circularly polarised light. As can be seen in 



Fig. 7.1 , we have taken the K valley (ft") to be lower than the K' valley (JJ-) in energy. 
This valley splitting leads to a possible inequivalence of the A and B sublattices, so 
that results may be different for an impurity located on the A sublattice and one 
located on the B sublattice. However, the same results are always obtained when 
moving the impurity from the AtoB sublattice, whilst simultaneously switching the 
energies of the K and K' valleys. Let |A/imi,A/'2 772-2) = '^Mxmx^Mim-^^ ~ ~^)- 
least energetic optically active inter LL transition takes the form |1 fff m, fl^t fn) 
for m = 0, 1, . . ., which has Jz = 1- The only transitions, which are in resonance 
are other transitions from LLs n = to n = 1 and those from LLs n = — 1 to 
n = 0. We do not correct the CRE for the valley splitting, setting hujy = 0, 
since fiLOy « huc- Thus transitions with pseudospin flips may be considered in 
resonance with those without. Spin flip transitions are not involved in the CE, since 
neither the e-h Coulomb interaction nor the interaction with the impurity flips the 
spin. The transition |1 fl"! m, fl^t "rn) is mixed via the Coulomb exchange interaction 
to the optically dark transitions |0 |lt — 1 J|t + 2), |0 Hi m, — 1 '(|i m + 2) and 
|0 JJ-i m, —1 JJ-i m + 2). Here the oscillator quantum numbers are determined by the 
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conservation of J,. 



According to the selection rules given in Eqs. (7.12)-(7.15), a particle can 
only be scattered to or from a state with r (r =J],) by an impurity on the B 
sublattice, when |n| = m (|n| — 1 = m). Hence the latter three transitions are not 
mixed with any others by the impurity interaction. The hole in |1 fff m, fl^t rn) 
only interacts with the impurity when m = 0, in which case its quantum state is 
unaltered. The electron only interacts with the impurity when m = 1 and the cor- 
responding excitation is |1 fft 1,0 1). This is scattered to either |1 fft 1,0 fft 1) 
(unchanged) or to |1 JJ-t 0, fft 1) by the electron-impurity (e-imp) interaction. Note 
that = 2 for |1 |Lt 0, Ht 1>, so excitations of the form |1 |Lt m-, Ht "i- + 1> must 
be included. Such transitions are mixed by the direct e-h Coulomb interaction to 
those of the form |0 Jjt ~1 111"'- + 3), which are not mixed with other types of 
transition by the impurity. We have now obtained the set of six excitations illus- 
trated in Fig.[7l], |1 Ht ?7i,0 Ht m>, |0 |Lt m, -1 JJ-f + 2>, |0 Hi m, -1 ffi m 2>, 
|0 |li m, -l\\,im + 2>, |1 lit "i,0 Ht "T- + 1>, |0 |lt -I ^] m + 3>, such that each 
member excitation is only mixed to other members of the set and which contains 
the only excitation optically active in the cr"*" polarisation for this filling factor. One 
could of course construct different sets of mutually mixed excitations, but these 
would be dark and the aim of this work is to focus on the optical problem. 

It is now straightforward to construct the Hamiltonian matrix and to diag- 
onalise it numerically. This requires cutting off the sum over m at a finite value. 
This is entirely justified physically, since we seek excitations, which are localised and 
only have large amplitudes for low m values. We include the first 30 m-values in 
our calculations, so that the total matrix size is 180 for the example detailed above. 
Changing the cut off value for m close to m = 30 does not change the energies 
appreciably, so that the errors are within symbol size for all results presented below. 



7.3 Results 

7.3.1 Coulomb impurity 

A certain duality exists in our results as a consequence of particle-hole symmetry in 
graphene. The energies and absorption strengths are the same under simultaneously 
changing the sign of the filling factor [v <-> —ly), the direction of light polarisation 
{a~^ <-> cT^), the sign of the angular momentum projection (J^ <-> —Jz) and the sign 
of the impurity's charge (donor, D"*" acceptor, A^). An example of this is shown 



in Fig. 7.2 The neutral CEs which are localised on the impurity appear as discrete 
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T non-degenerate 
□ doubly degenerate 



(al 



-5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9 10 
J 




Figure 7.2: neutral CEs localised on (a) a donor impurity a,t u = —1 and (b) an 
acceptor impurity at = 1. Energies are given relative to ei + -Eg^f — O.TSE'o in 
units of £^0 as a function of the generalised angular momentum projection, J^. The 
hatched area of width O-TSi^o represents the continuum of extended neutral CEs. 
Quasibound states within the continuum are not shown. Insets show four branches 
of resonantly mixed inter-LL transitions conserving spin and pseudospin. 



states outside the continuum of delocalised neutral CEs. All the presented results 
are for e/ejmp = 1 unless otherwise stated. 

Figure |7.2| ^a) shows for v = —1 four low-energy branches of neutral CEs 
localised on the -D"*" for > 1; two of these branches are degenerate. Their nature 
is explained as follows. For large positive J^-, the hole is on average much farther 
away from the impurity than the electron (see Eq. (4.13)). Therefore, the e-D'^ 
attraction dominates over the h-D~^ and e-h interactions. Hence we find branches 
with asymptotic Jz » 1 energies equal to ei + Eq^^ + Vcj/m = 0; 1; • • •)• 
an example, notice the three branches approaching energy — 0.25-E'o and the single 
branch approaching zero energy in Fig. |7.2[ a). These originate, respectively, from 
the three n = —1 ^ n = transitions (denoted hereafter as r_io) and from the 
single Tqi transition for u = —1, since Vco^so ~ ~Eo and Vcirso ~ — 0.75£'( 



0- 



Similar asymptotic behaviour can be seen for other filling factors in Fig. 7.3 



The neutral CEs above the band develop for Jz < 0, when the hole is closer to the 
than the electron (or the electron is closer to the than the hole). Such unusual 
excited states are bound in 2D because of the confining effect of the magnetic field. 
Due to the duality, results for the A'^ at v = —1,0 (not shown) can be obtained 
from those for the by changing Jz —> — Jz and v —u. 

For the case = 2, we have the full SU(4) symmetry, so I have marked 
on the true degeneracies of the states in Figs. |7.3[c) and (d), as determined in 



Section 5.3.2 For the other filling factors, the symmetry between different single 
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Figure 7.3: neutral CEs localised on the at (a) = 0, (b) = 1, (d) z^ = 2, and 
(c) localised on the A'^ aX v = 2. 
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particle spin/pseudospin flavours is broken, but only partially. Hence if the energies 
of excitons with spin and or pseudospin flips were calculated, I would expect some 
of them to be degenerate with the energies shown. However, this is left to further 
work and the degeneracies given are due to excitons with no spin or pseudospin flips 
only. 



Figure |7.4| shows the optical properties of the localised states for various 
different parameters. The duality between states in Figs. |7.4[ a) and (b) and also 
between states in Figs. |7.4[ c) and (d) is evident. For this reason, we restrict the 
discussion of results to the D~^. Note that for each polarisation, there are two 
branches of dark states indicated by diamonds. Optically active localised states are 
shown using circles, with size which scales as the modulus squared of the dipole 
matrix element, d = {MiM2Jz\'H+\u}. Two types of localised states can be optically 
observed: (i) truly bound states, which are split off the continuum and have nor- 
malisable wave functions with exponentially decaying tails, (ii) quasibound states 
within the continuum, which have high probability amplitudes on the impurity and 
long-range oscillating tails. The latter may exhibit asymmetric Fano-type optical 
signatures [242 ], which is beyond the scope of the present work. 

For both polarisations (Figs. |7.4| ^a) and (c)), the upper branch originates 
mostly from the Tqi transitions with some small (zero at = 2) admixture of the 
r_io. With increasing z^, the number of (strongly mixed by the e-h exchange) Tqi 
transitions increase (see Fig. 7.2 insets), which leads to the enhanced contribution 
of the repulsive e-h exchange interactions. This explains the blue shift of the up- 
per branch to higher energies with increasing v. Also, its optical strength, 
increases for the polarisation (Fig. 7.4 ^a)) and decreases for the polarisation 
(Fig. |7.4[ c)). This is explained by the fact that the Tqi transitions are optically ac- 
tive in the polarisation while T^io transitions are dark. Conversely, the strength 
of the upper branch in the polarisation originates solely from the T-iq. There 
are fewer of them with increasing v, and eventually the upper branch becomes 
completely dark in the at v = 2. Similarly, the lower-energy branch mainly 
originates from the r_io transitions with some small admixture of the Tqi. Its red 
shift to lower energies with increasing u is explained by the decreasing number of the 
T_io transitions leading to the decrease of the repulsive e-h exchange contribution. 



Figure. [73| addresses the effect of impurity screening, as discussed in Section 
7.1.1 The energies of states localised on a impurity for v = \ and the o""*" 
polarisation are plotted as a function of e/eimp. As e/eimp increases, the impurity 
screening is reduced and the states are pushed away from the band. 
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Figure 7.4: Evolution with filling factor v of energies and optical strengths of neutral 
CEs localised on (a) the with = 1 active in the cr~^ polarisation, (b) the 
with Jz = —1 active in the polarisation, (c) the with = — 1 active in 
the polarisation and (d) the A'^ with = 1 active in the polarisation. The 
optically active states are indicated by circles with sizes ~ Their de generacies 



are not given for clarity; this information is contained within Figs. 7.2 and 7.3 The 
diamonds represent optically dark states. The degeneracies implied by the close 
overlapping symbols are due to excitons with no spin or pseudospin flips only. The 
dotted lines are guides to the eye. Inset of (c): Dipole strength \d\'^ versus energy for 
= 0. The spectra were convoluted with a Gaussian of width O.OSEq. The arrow 
indicates an impurity-related feature below the CRE. 
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Figure 7.5: Energies and absorption strengths of neutral CEs localised on the 
for filling factor u = 1 and light polarisation (T+ as a function of e/eimp- 



7.3.2 5-function impurity 

As for the case of a Coulomb impurity, states localised on a 5-function impurity 
possess a high level of symmetry. Their energies and optical properties are the same 
under simultaneously changing the sign of the filling factor (i^ <r-> —ly), the direction 
of light polarisation {a~^ a^), the sign of the impurity potential (Vq <^ —Vq) and 
the sublattice on which the impurity is located {A B). In addition, changing 
only the impurity from one sublattice to the other, yields identical results when the 
filling factor is even, but qualitatively different results when it's odd. This is due 
to the valley splitting, which allows the sublattices to become inequivalent if the 
valleys are not equally filled. 



Figure [7^ shows the energies and absorption strengths of localised states for 
a selection of values of filling factor, light polarisation and sublattice plotted as a 
function of the impurity strength, Vq. In each plot only one of the parameters dif- 
fers from those in Fig. |7.6| [a) to see how this affects the states. As for the Coulomb 
impurity case, the localised neutral CEs occupy discrete energy levels outside the 
continuum of extended neutral CEs. In contrast to the Coulomb impurity, the ener- 
gies of states localised on the (5-function impurity are non-degenerate. In particular 
there is no evidence of the SU(4) [15]-plet seen in clean graphene and the Coulomb 
impurity case for v = 2. This is because the 5-function impurity breaks the sym- 
metry between the sublattices and hence the valleys. In contrast, the long-ranged 
Coulomb interaction varies slowly on the scale of the lattice spacing and preserves 
the sublattice symmetry. 
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Note that the continuum has width Eq in Figs. |7.6| [|a),(b),(d) which have 
u = —1 and width O.JBEq in Fig. |7.6[c) which has u = 0. The band width is 



determined by which transitions are mixed; this is shown in Figs. 7.1 and 7.7 Notice 
that in ah cases, there are the four transitions with no spin or pseudospin flips, which 
are mixed by the exchange e-h Coulomb interaction. If these were not mixed to any 
other transitions, which is the case for the Coulomb impurity, the continuum would 
have width 0.75£'o, as seen in Section 4.8.3 In some cases (specifically those shown 
in Figs. |7.7[ b),(c)) the 5-function impurity admixes transitions with a pseudospin 
flip, which are not otherwise mixed by the Coulomb e-h interaction to any other 
transitions; these do not alter the width of the band. For fllling factors v = —1 and 
V = 1 (latter not shown) however, the (5-function impurity also admixes pseudospin 
flip transitions, which are each mixed to a single other transition by the direct e- 
h Coulomb interaction. This broadens the band from 0.75ii^o to Eq (see Section 



4.8.3). 



Our results indicate critical Vq values for the formation of localised excitations, 
although within our approach, it can be difficult to resolve what may be a low 
energy bound state from the lower continuum edge. The critical value for which 
bound states appear varies slightly depending on the filling factor, light polarisation 
and sublattice of the impurity. For example, for a system with filling factor v = —1 
illuminated by cr"'" circularly polarised light with an impurity on an A sublattice 
(such a system will henceforth be denoted by A, fj"*", = — 1), there are no bound 
states for -30eV <W< 120eV (for B = lOOT). Thus it takes a very strong 5- 
function impurity potential to localise excitations. This is due to the fact that the 
(5-function impurity only couples to a small subset of the basis states of a neutral 
CE, which have m = 0, +1. Strictly speaking, such large energies are greater than 
the width of the vr band, where the electrons can be reasonably treated as massless 
Dirac fermions and so should not be treated as a small perturbation. 

For both light polarisations, impurity positions and for all four fllling factors, 
larger magnitudes of Vq are required to push states above the band than to pull them 
below. States below the band form when either the electron or hole is nearer to the 
impurity and attracted to it; in this case the additional e-h attraction also lowers the 
energy. States above the band form when one particle is held nearer to the impurity 
due to magnetic confinement, but is repelled by it, whilst the other particle is further 
away. In this case the e-h attraction works against this, trying to lower the energy, 
so that a larger impurity strength is required to overcome this. 

The bound states move further from the band as the impurity strength in- 
creases, in accordance with intuition. An interesting question is what happens to 
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Figure 7.6: Evolution of energies and optical strengths of neutral CEs localised on a 
(5-function impurity located on either an ^ or .B sublattice site as a function of the 
impurity strength, Vq, for various filling factors and light polarisations (see labels 
on plots). The parameters B = lOOT and e = 5 are used. The yellow shaded area 
represents the continuum of neutral CEs. The energies are renormalised, so that the 
CRE, which equals ei + E^;i - O.TS-Eo, is at zero energy. Optically dark states are 
represented by diamonds and optically active states (in the specified polarisation) 
by squares with sizes ~ \d\'^. 
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Figure 7.7: Diagrams indicating which neutral CEs are mixed when the (^-function 
impurity is on the A sublattice and the sample is illuminated by (a) polarised 
light with = — 1, (b) (T+ polarised light with z/ = and (c) polarised light with 
u = -1. 
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Figure 7.8: Evolution of energies and optical strengths of neutral CEs localised 
on a (^-function impurity located at a carbon site on (a) the A sublattice and (b) 
the B sublattice site as a function of the filling factor v. The yellow shaded area 
represents the continuum of neutral CEs. The energies are renormalised, so that the 
CRE, which equals ei + - ^.IbE^, is at zero energy. Optically dark states are 
represented by diamonds and states active in the a"'" polarisation by squares with 
sizes ~ 



the bound states when they approach the band. One possibility is they continue 
to exist within the band as quasibound states, with a high probability to exist on 
the impurity in addition to long-range oscillating tails. Although we detect pos- 
sible signatures of such states, our method is not accurate enough to claim their 
existence. The hypothesis may be supported however by the observation that in all 
cases the number of branches below the band for positive (negative) Vq equals that 
of branches above the band for negative (positive) Vq. 

In most cases, for large enough |Vb| values, bound states may be found both 
above and below the band for both Vq > and Vq < One exception (see Fig. 
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^a)) is the A, cr"*", u = —1 system. In this case, the impurity selection rules 
forbid the hole to interact with the impurity, so that all diagonal matrix elements 
have the same sign as Vq and states above the band are only seen for large enough 
positive Vq and states below the band are only for large enough negative Vq. This 
behaviour is of course mimicked by the symmetric system B, , v = 1 except that 
the electron is forbidden to interact with the impurity, so that states above the band 
only appear for large enough negative values of Vq and states below the band only 
for large enough positive Vq values. The other exception is the pair of symmetric 
systems A, c^, v = 2 and B,a^ ,v = 2, where in both cases only the hole interacts 
with the impurity so that states above the band only exist for large enough negative 
values of Vq and states below the band only for large enough positive Vq values. The 
number of branches due to the e-imp and the number of branches due to the /i-imp 
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interaction can be understood by studying the transitions involved and the impurity 
selection rules. 



Figure 7^ shows the energies and absorption strengths of states localised 
on a 5-function impurity with Vq = SOOeVo^, as function of filling factor v, when 
the sample is illuminated by polarised light. The cases when the impurity is 
located on the A and B sublattice are compared. One can see that the energies and 
absorption strengths are identical for filling factors = 0,2, but differ for v = +1, 
as previously mentioned. 

7.4 Discussion 



It is interesting to compare Figs. |7.4[ a) and 7.8 and see how the results differ for 
the Coulomb and 5-function impurities. For the 5-function impurity, there are more 
localised states and they do not seem to evolve smoothly as a function of filling factor 
in contrast to the case of the Coulomb impurity. This is because for the Coulomb 
impurity there are always four transitions with the same Jz for any filling factor, 
whereas for the short-range impurity, where the intervalley scattering is significant, 
changing the filling factor may change the number of transitions that are mixed and 
their Jz values, which influences the number of branches of bound states that will 
be seen. The energies may be degenerate for the Coulomb impurity, but not the 
5-function impurity. This can again be attributed to the greater complexity of the 
transitions involved for the J-function impurity. 

One may also compare the results for the 2DEG with those for graphene in 



the case of a Coulomb impurity 240 . They are qualitatively similar with branches 
of neutral CEs appearing both above and below the continuum with energies ~ Eq. 
However, there are more branches in graphene than in the 2DEG, due to the the 
additional LL pseudospin splitting. 

A possible way of detecting the localised neutral CEs predicted above is 
via polarisation resolved spectroscopy experiments, in much the same way as was 
described for charged CEs. The localised neutral CEs below the band may be 
visible as additional lower intensity peaks in the absorption spectrum below the 
main CRE peak. This is shown schematically in Fig. |7.4[ c), inset. Our single- 
impurity theory is applicable to samples with finite impurity density nimp < I/tt/^, 
i.e. when the mean separation between impurities exceeds the size of localised neutral 
CEs. The intensity of impurity peaks will then be / ~ ?T-imp|c^P- In current samples, 
nimp ~ lO^^cm^^ K l/vr/^ 92,243 . Hence it is reasonable to expect that with 



the rapid advances being made in sample fabrication, clean enough samples will be 
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available in the near future. Experimental verification of our results would enable 
the presence of short-range and long-range (Coulomb) scatterers to be distinguished, 
thus contributing towards the ongoing debate regarding the nature of disorder in 
graphene. 
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Chapter 8 

Conclusions and outlook 



This work has focused on several aspects of disorder and many-body interactions 
in quantum systems. The analytical treatment of cither of these issues is fraught 
with difficulty and much of the progress is made using numerical techniques. An 
even more challenging problem is the interplay between the two effects, especially 
when neither one is clearly dominant. Many theories focus on one or the other, 
but usually both will be present. In this thesis, three systems were examined: i. 
the Anderson model of BCC and FCC lattices in which disorder is modelled using 
random matrices and electron-electron interactions are ignored, ii. an exciton in a 
nanoscopic Aharonov Bohm ring, which is disorder-free, but where the electron and 
hole interact via a (5-function potential and iii. a monolayer of graphene in a strong 
perpendicular magnetic field, where excitations are collective, due to the screened 
Coulomb electron-electron interactions. There are a few other themes that connect 
the projects. For example, impurity (disorder) induced localisation of electrons is 
considered in i. and of collective excitations in iii. In projects i. and ii. quantum 
interference of the particle's wave function plays an important role. In this final 
Chapter, I shall summarise the main results arrived at in the previous Chapters, 
discuss the remaining open questions and suggest directions in which the work could 
be taken forward. 

In the Anderson transition study, we implemented the transfer-matrix method 
to calculate localisation lengths for electrons in the 3D BCC and FCC lattices with 
onsite disorder. Finite size scaling was performed on these results for a range of dif- 
ferent disorders, energies and system sizes in order to obtain critical disorders (at a 
constant energy), critical energies (at a constant disorder) and the critical exponent, 
which governs the behaviour of the correlation length as criticality is approached. At 
energy £' = 0, we found an increase in the critical disorder, as the coordination num- 
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ber of the lattice increases, in accordance with intuition: Wc = 16.54 [l8 31 [32 135| 
for the SC lattice (coordination number 6), 20.81 for the BCC lattice (coordination 
number 8) and 26.73 for the FCC lattice (coordination number 12). All the calcu- 
lated critical exponents for both lattice types, however, are v ~ 1.5, which is also in 
agreement with previously calculated results for the SC lattice. This is because it is 
a universal localisation property of a 3D system in the Gaussian orthogonal univer- 



sality class, and hence should not change by the scaling theory of localisation 13 



In addition to high precision calculations, we used some of our slightly lower preci- 
sion data to obtain the energy-disorder phase diagrams for BCC and FCC lattices. 
The diagram for the BCC lattice has a similar shape to that previously calculated 
for the SC lattice, but the diagram for the FCC lattice is assymetric about E = 0, 
due to its non-biparticity. 

It would be interesting to apply the techniques outlined in this project to 
other lattices. Graphene would perhaps be a good candidate. According to scaling 
theory, as a 2D system, the electronic wave function should be completely localised 
in the presence of any finite disorder, leading to insulating behaviour. Indeed this 
was found to be the case in numerical calculations of the localisations lengths for 
the honeycomb lattice |141| . However, the majority of experiments on graphene 
indicate that it is in fact a good conductor. In addition, other numerical studies 
claim to have found an MIT [244 , 245 . This prompts the question of how large the 
system size should be before localisation is observed and how this compares with 
the typical sample sizes used in experiments. 

In the Aharonov-Bohm study, we developed a simple ID model for an exciton 
in a nanoring with an in-plane electric field. We were able to express the problem 
in terms of an eigenequation, where the matrix depends on the excitonic energy. 
Correct energy and wave function amplitudes are provided, when the eigenvalue is 
close to 1. We observed clear AB oscillations with period equal to the magnetic flux 
quantum. The amplitude of these oscillations first increased as a function of electric 
field and then rapidly decreased once the exciton had been pulled apart. Most 
importantly we also observed oscillations in the oscillator strength as a function of 
magnetic fiux. Beyond the critical electric field value, these became inverted i.e. the 
maximum became a minimum. Interestingly the minimum reached zero, indicating 
a probability zero for the recombination of the electron and hole and the emission of 
a photon. Thus by tuning the electric and magnetic fields, we found a mechanism 
whereby the exciton could be manipulated between being light and dark. This could 
have potential applications in excitonic data storage and the trapping of light. Two 
recent experiments have observed excitonic AB oscillations in the presence of an 
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electric field [176 , 177] . Both an inversion in the oscillation pattern at high enough 
electric fields [176 and a dramatic reduction in the emission intensity for particular 
values of the external fields 177] were observed, in accordance with our predictions. 

An obvious extension to this work would be to create a more realistic ring 
model by introducing a width and/or height to the ring in order to see how this 
might affect the oscillations. Another interesting avenue of research would be to con- 
sider electrons (Cooper pairs) in superconducting quantum rings instead of excitons 
in semiconducting devices. In the thermodynamic limit, the free energy [55] and 
all thermodynamic quantities [57| should have periodicity <l*o/2. However, this limit 
is not reached in mesoscopic systems. The excitonic Bohr radius will be replaced 
by the coherence length (roughly, the distance between electrons in a Cooper pair) 
and the penetration depth of the magnetic field may also play a role. In s-wave 
superconducting rings much smaller than the coherence length, oscillations in the 
supercurrent with period <I>o are now well established 2461 . In d-wave superconduc- 



tors a <I>o current component has also been found [247]. It would be interesting to 
look at the problem for p-wave superconductors. These are particularly interesting 



in this context, since $o/4-quantised vortices exist there [248] . One topic of investi- 
gation could be how the periodicities vary with the ring size close to the coherence 
length. 

In the graphene project we made a comprehensive study of collective excita- 
tions of the monolayer system in the presence of a magnetic field. I first reviewed 
the work that has been done by others on neutral collective excitations of pristine 
graphene. These can be thought of as excitons, but the term "collective" is used to 
emphasise that several excitations are mixed i.e. there is a superposition, due to the 
two-body Coulomb interaction. We then considered the localisation of such neutral 
collective excitations on a single impurity. Both long range (screened Coulomb) and 
short range (5-function) impurities were treated. The localised excitations occur 
in discrete states outside the band of extended collective excitations. They form 
branches for the case of a Coulomb impurity, although for a 5-function impurity 
things are somewhat (and maybe counterintuitively) more complicated, as different 
excitations are mixed for different filling factors. The main difference between the 
two impurity types is that the Coulomb impurity cannot scatter between the two 
valleys, K and K', in graphene, whereas the 5- function impurity is able to cause 
this large change in an electron's momentum. Certainly the Coulomb impurity is 
much more effective at localising the collective excitations. We observed a symmetry 
between the localised states. Their energies and optical properties are the same un- 
der simultaneously changing the sign of the filling factor (z^ ~^); the direction of 
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light polarisation ((T+ ^ a^), the sign of the impurity potential (Vq ^ —Vq) and the 
sublattice on which the impurity is located for the case of the (^-function impurity 
{A <r-> B). We hope that our work, particularly if it is experimentally verified, may 
contribute towards the ongoing debate regarding the nature of disorder in graphene. 

We also looked at charged collective excitations, which may be thought of 
as electron-electron-hole (X^) or electron-hole-hole {X~^) complexes. Such discrete 
states may exist below the continuum for very nearly full {X~^) or very nearly 
empty (X^) Landau levels, making them energetically favourable; these two types of 
excitation are connected by a symmetry relation similar to that for localised neutral 



collective excitations. The very recent observation of plasmarons in graphene 236 



which are the zero magnetic field equivalent of our X~^ states, makes us optimistic 
about future experimental observation of the states we predict. We found that a few 
are optically active and as such may be visible in optical spectroscopy experiments 
as a low energy shoulder to the cyclotron resonance or additional peak below it. 

In our calculations we made use of both the geometrical and dynamical sym- 
metries of our system. There is always rotational symmetry about a perpendicular 
axis through the origin, where the impurity is located. This enables us to label 
the collective excitations by a generalised angular momentum, J^, which has an 
orbital and sublattice component. For the pristine system, translational symme- 
try meant that the charged collective excitations could be labelled by a /c quantum 
number associated with the operator of magnetic translations. By the dynamical 
symmetry, I mean the SU(4) symmetry associated with the presence of two spin 
and two valley pseudospin degrees of freedom, yielding four possible flavours in to- 
tal. For a fully or approximately fully filled Landau level, a symmetry arises due to 
the possible exchange of these four equivalent flavours. (Note that we do not have 
the (5-function impurity in mind here, since in this case the impurity introduces 
an inequivalency between the two pseudospin projections). With this in mind, we 
introduced an analogy between neutral collective excitations and mesons and also 
between charged collective excitations and baryons, where only the first and second 
generation quarks are considered. Using this, we could employ Young diagram tech- 
niques to understand and correctly predict the degeneracies of our states. We also 
produced multiplet plots similar to the style of the SU(4) multiplet polygons seen in 
particle physics, as a way of visualising our states and their dynamical symmetries. 

There are still several unanswered questions within the work presented here. 
Most importantly, we would like to fully understand the dualities that exist between 
our results for different filling factors, light polarisations and impurity potentials. 
There should of course be a similarity transformation connecting any two Hamilto- 
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nians describing symmetric cases. Another way of taking the graphene work further, 
would be to examine collective magneto-excitations of the bilayer system, which has 
already attracted a lot of attention 103 . Its neutral collective excitations (magne- 
toplasmons) have already been studied for a clean system [249 , 250| , but not, so far 
as I'm aware, for a system containing impurities. 
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Appendix A 

Connectivity matrices 



Given below are explicit expressions for the connectivity matrices for BCC and 
FCC lattices with M = 3. Recall that Cj is the connectivity matrix describing the 
connections of the i"^ slice to the (i — 1)*'^ slice. Element Cjk of the connectivity 
matrix equals 1 if site j in the i^^ slice is connected to site k in the {i — 1)**^ slice; 
otherwise Cj^ = 0. The boundary terms are indicated in italics. For the BCC lattice 
for odd layers, 
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For the FCC lattice for odd and even layers 



( 1 


1 

















1 


\ 





1 


1 

















1 


1 





1 


1 




















1 





1 


1 




















1 





1 


1 




















1 





1 


1 




















1 





1 


1 




















1 





1 


1 



















1 





1 } 



(A.3) 



In all cases, i is a positive integer. 
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Appendix B 



Coulomb matrix elements for 
the 2DEG 



As seen in Eqs. (4.10) and (7.10), graphene matrix elements of the Coulomb inter- 
action can be expressed in terms of those for the 2DEG. The method for calculating 
such 2DEG matrix elements is explained in Ref. [218 and outlined here for conve- 



nience. Let us first examine the impurity matrix elements 

Vnm = (nm\Vc\nm) = 
The general technique is to express Vc (r) in terms of its Fourier transform 



Vc(.) = ^ (B.2) 

€imp9 



and the operator exp (iq • r) in terms of displacement operators 



/■ \ , -i a' qa \ I iqb^ iq*b ■ 

exp(.q.r) = -P(-:^ + ^J-P(:^ + ^l (B-3) 



V2y VV2 

where q = {qx + iQy) ^B- Since the two different sets of Bose ladder operators com- 
mute, 

{nm\Dn (a) {/3) \nm) = <n|D„ (a) \nXm\D.m (/3) \m). (B.4) 
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One can show using the Baker-Campbell-HausdorfF formula and exphcit expressions 
for Laguerre polynomials that performing the real space integrals gives 

(n'\Dia)\n)=\ V r -1' . ,^ ^^'^^ 



Using the above one can show 

ZEr '"^ 



Vr,. 



—S x-V2e-L„(x)L„(x), (B.6) 
Jo 



which can be calculated explicitly for particular n, m values. In fact the only ele- 
ments used in our calculations are 

_ ZEo r{m + l/2) _ 4m-l 

VOm — — 1= j , Vim—-. 7;V0m- [^■') 

/vr m! Am — 2 



Now consider the two-body Coulomb interaction matrix elements 

(B.8) 

The techniques used are the same as those for the impurity matrix elements. Two 
important examples are the direct and exchange matrix elements 

rm\sm', _ I minilmius! 
Urmism,2 - — : — 7777 — : — 777 1^5. yj 



rOO 

Jo 



sm\rm2 




mini!min2! 



7r(mini -|- t)!(min2 -|- 1)\ 



(B.IO) 



dx x'-'+'-y^-^^L. 

JO 



s—r I 
r 



(a;)) -^^niirii (^)-^min2('^)' 



where t 



m'^l = |m2 — mini = min(mi, m'^) and min2 = min(m2,m2). 
Eq. (4.28) gave an expression for the correction, E'se {n-,n'), to the energy 
of an electron in the n^^ LL of graphene due to exchange with electrons in the n'^^ 



LL. By Eq. (4.10), E^e {n, n') can be expressed in terms of 2DEG matrix elements. 
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Thus we calculate 



= S J^q^{^(.)<n;|l)n (-^) |mXn',|i)„ (^) |n,>(B.ll) 



X 



{m\b^ ( --^j |m'><m'|D„ (^-^j |m> [/4vr2 



= Jdq't/((/)<n;|^„ (--^^ \ni){n'^\bn |n2>/47r2, 

where f/((?) = 2Txe^ jtq is the Fourier transform of the Coulomb potential, C/(r). Eq. 
(B.5) may now be utilised to obtain results for particular ni , n'^^ , n2 , n2 values. 
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Appendix C 

Spin and pseudospin eigenstates 



The excitonic spin and pseudospin eigenstates are 

|5 = 0, T = 0; 5, = 0, = 0> = ^ (| tt> + Ut ^t> + I Hi) + I U U)) 

(C.l) 

\S = 0,T = 1;S, = 0,T, = -1> = -L (I tt> + I , (C.2) 

\S = 0,T = l;S, = Q,T, = Q}=^i\nm-\ ^T> + I ti> " I U U}) , 

(C.3) 

\S = 0,T = l-S, = 0,T, = l) = ^i\ilU^)+\ Hi m , (C.4) 

\s = i,T = o-s, = -i,T, = o> = (I ny + 1 u m) , (c.5) 

|5 = 1, T = 0; 5, = 0, T, = 0> = ^ (I Ht frt> + I ^t> " I fTi 1li> " I U U)) , 

(C.6) 

|5 = i,r = 0; 5, = 1, = o> = (I ti> + I , (C.7) 

|5 = 1,T = 1;5, = -1,T, = -1> = I U m, (C.8) 
|5 = 1, r = 1; 5, = -1, = 0> = (I tt> - Ui ^t» , (C.9) 

\S = 1,T=1;S, = -i,r, = i> = I ^T>, (c.io) 

\S = 1,T=1;S, = 0,T, = -1> = -L (I ^1 _ I , (c.ll) 

|5 = 1, T = 1; 5, = 0, = 0> = ^ (I tt> - I ^T> " + I , 

(C.12) 
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\S = l,T = l;S, = 0,T, = l) = ^{\ilU^)-\ Hi U)) , (C.13) 

\S = 1,T=1;S, = i,r, = -1> = I Hi), (C.14) 

|5 = i,r = = i,r, = o> = (I - I n Hi}) , (C.15) 

\S = 1,T = 1;S, = 1,T, = 1} = \ n Uy (C.16) 
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